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CONVEXITY  OF  FUNCTIONALS  BY  TRANSPLANTATION 


ty 

C.  Polya  - M.  Schiffer 


In  the  present  paper,  we  wish  to  discuss  the  dependence  of  various 
functionals  on  their  domain  of  definition.  The  functionate  considered 
will  be  defined  by  certain  extremum  problems  or,  equivalently,  by  boundary 
value  problems  for  partial  differential  equations  of  elliptic  type.  In 
either  case,  it  is  difficult  tc  compute  the  functional  for  an  arbitrarily 
given  domain,  while  explicit  solutions  are  availably  if  at  all,  only  in  the 
case  that  the  domain  of  definition  is  of  simple  or  symmetric  structure. 

There  arises  the  question  of  utilizing  the  knowledge  of  the  functional 
for  a few  special  domains  in  order  to  obtain  information  about  the  same 
functional  in  the  general  case.  Two  different  methods  have  been  applied 
successfully  for  this  purpose;  namely,  the  method  of  transplanting  extremum 
functions  and  the  method  of  -.filiation. 

The  method  of  transplantation  can  be  applied  if  the  functional  in 
question  is  characterized  as  the  extremum  value  of  another  functional  over 
a certain  function  class  with  respect  to  the  domain  of  definition  D.  Suppose 
that  we  know  the  function  in  D which  yields  the  correct  extremum  value  of 
the  functional;  if  D'  is  another  domain  referred  to  D by  a one-to-one 
transformation  of  a sufficient  degree  of  continuity,  we  may  consider  the 
extremum  function  of  D when  properly  transplanted  as  a new  function  in  D' . 
Here  it  will  not,  in  general,  be  the  extremum  function  for  the  problem 
considered  with  respect  tc  D';  but  it  may  still  be  a competing  function  for 
the  extremum  problem  and,  in  this  event,  the  value  it  imparts  to  the  D/-func- 
tional  may  be  used  as  an  estimate  for  the  actual  extremum.  Finally,  we  may 


obtain  in  this  way  a comparison  between  the  extrema  with  respect  to  D and 
to  D'  so  that,  if  the  b- extremum  is  knowq  a bound  for  the  D'-extremum  has 
been  found.  The  method- of  transplanting  tne  extremum  function  has  been 
used  systematically  in  [15  ^ and  has  provided  numerous  inequalities  and 
theorems . 

The  method  of  variation  starts  with  the  assumption  that  a functional 
is  known  for  a given  domain  D and  provides  a formula  for  the  infinitesimal 
change  of  the  functional  in  case  the  domain  D of  definition  is  changed 
infinitesimally.  The  variational  formula  may  serve  to  characterize  those 
domains  for  which  the  functional  at  t s in  s c xin  omum  values  under  various 
side  conditions.  It  happens  frequently  that  the  formal  structure  of  the 
variational  formula  indicates  monotonicity  or  other  growth  characteristics 
of  the  functional  studied.  In  [61  the  fact  was  stressed  that  various 
convexity  properties  of  functionals  could  he  established  by  considering 
the  second  variation  of  the  functional  with  the  domain,  and  several 
applications  of  convexity  results  were  given. 

*In  this  paper  a certain  synthesis  between  the  above  two  methods  is 
intended.  We  shall  use  the  method  of  transplantation  within  a family  of 
domains  D(t)  which  arise  from  one  given  domain  D»D(l)  by  a one-parameter 
transformation.  In  this  way,  we  are  able  to  obtain  variational  formulas 
and  convexity  statements  for  the  functionals  considered.  The  .information 
obtained  is  more  extensive  thar  that  found  by  one  discrete  transplantation 
of  the  extremum  function  while  the  results  available  by  infinitesimal 
variation  are  improved  to  statements  in  the  large.  We  consider  in  this 
paper  only  a few  simple  types  of  transformation  such  as  conformal  transfor 
mations  in  the  plane  case  or  one-sided  stretchings  in  the  two-  and  three- 
dimensional  case.  But  it  is  easy  to  apply  the  same  method  to  other 
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transformations,  as  well,  and  various  analogous  results  can  then  be  obtained. 

7fe  indicate  briefly  the  partial  one-sided  stretching  of  a domain  which 
consists  in  a stretching  of  the  part  of  the  domain  that  lies  to  one  side 
of  a given  line  or  plane  while  the  rest  of  the  domain  is  left  unchanged. 

By  using  proper  partial  stretchings,  many  general  variational  formulas  can 
be  derived  by  the  transplantation  argument,  a fact  which  shows  the  significance 
of  this  simple  type  of  transformation. 

In  Chapter  I,  we  develop  our  method  for  functionals  such  as  torsional 
rigidity,  virtual  mass,  outer  conformal  radius  and  electrostatic  capacity. 

Most  of  these  quantities  are  defined  by  two  complementary  extremum  problems 
leading  to  upper  and  lower  bounds  for  the  functionals.  In  these  cases,  our 
method  leads  to  corresponding  complementary  convexity  statements  which  are 
of  value  in  the  numerical  treatment  of  the  quantity  considered. 

In  Chapter  II,  we  discuss  a useful  theorem  of  Poincarl  which  permi+r., 
an  easy  simultaneous  estimation  of  the  N first  eigenvalues  of  a general 
type  of  eigenvalue  problems.  This  lays  the  groundwork  for  the  application 
of  our  method  to  eigenvalue  problems . 

In  Chapter  III,  the  convexity  of  various  combinatiohs  of  eigenvalues 
is  studied  in  the  case  in  which  the  domain  of  definition  is  deformed  either 
by  stretching  or  by  conformal  transformation.  In  the  problems  treated, 
the  convexity  statement  is  of  particular  interest,  since  the  eigenvalues 
considered  are  defined  as  minima  of  certain  quadratic  functionals  and  since 
from  the  definition  itself  only  upper  estimates  for  the  eigenvalues  are 
available  in  an  elementary  fashion.  Here  convexity  results  lead  sometimes 
to  lower  bounds  and  are  an  important  numerical  tool,  ffe  refer  to  [10]  where 
the  convexity  theorem  of  this  paper  has  been  applied  in  order  to  find  estimates 
for  the  membrane  eigenvalues  of  isosceles  triangles. 


■ 
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In  Chapter  IV.  we 
the  initial  dona  in  D ( 1 % 
ca^f j the  derivative  of 
be  fo'.uid  .at  the  point  t 
examples . 


indicate  the  use  which  car.  be  made  of  the  fact  that 
has  symmetry  properties.  It  is  shown  that  in  many 
the  functional  with  resuc't  to  the  parameter  t can 
-1.  Various  inequalities  are  derived  as  illustrative 


In  Chapter  V.  we  show  that  the  invnriar.ce  of  the  class  of  harmonic 
functions  in  l>  under  c-nformal  mapping  can  be  used  in  order  to  derive 
convexity  statements  for  some  interesting  fur.cti  nals  connected  with  the 
Green's  function  for  Laplace's  equation.  In  fact,  the  invariance  of  this 
Green's  function  under  conformal  mapping  provides  the  simplest  and  most 
frequently  used  instance  of  transplantation . 

Finally,  in  the  Appendix,  a numerical  application  is  riven  for  the 
torsional  rigidity  of  isosceles  triangles  and  rectangles.  This  example 
indicates  how  far  one  car.  go  miner ica lly  bv  a proper  utilization  of  a few 
comparison  and  convexity  theorems. 


H 


Charter  ' 


Transplantation  of  the  Extremal  Function 


1.  Torsional  Rir.iditv. 

1.1  We  consider  a s :m: ly- connected . bounded  domain  D in  the  (x.y)-plane 
and  define  a one-sided  stretch! r.g  of  D in  the  x-direction  as  follows.  We 


subject  the  ur.ole  ; lane  to  the  linear  transformation 
( 1 . 1 ) x " - tx  y ' -*  y , t > 0 , 

and  define  D(t)  ns  the  im .age  of  E under  this  stretching;  we  call  t the 
parameter  of  the  stretching  ar.d  allow  for  it  arbitrary  positive  values. 

The  original  domain  b may  also  be  denoted  by  D(l). 

We  wish  to  study  the  dependence  of  various  functionals  of  the  domains 
D(t'  upon  the  parameter  t and  to  demonstrate  a rather  general  method  which 
is  applicable  if  the  functional  in  question  is  defined  by  an  extremum 
problem.  We  explain  the  rr.etr.od  for  the  case  of  the  torsional  rigidity 
P(t)  of  the  domain  D(t'. 

The  functional  F is  defined  by  means  of  the  stress  function  f(x,y) 

'’or  the  domain  b.  f his  function  satisfies  the  differential  equation 


(1.2)  £f  « -2  . A - * - 

c3x2  dyc 

in  b and  vanishes  on  the  boundary  C of  D.  In  terras  of  f(x,y).we  have  [15': 

(2  J J f dx  dy)^ 


(1.3) 


j/  (V  f )*'dxdy 


Using  integration  by  parts  and  (1.2).  we  may  bring  ( 1 . 3 ) also  into 
alternative  forms 


(1.4) 


P ^ 

P - 2 ij  fdx 

D D 

The  advantage  of  the  definition  (i.3)  over  the  formulas 


(1.4) 


the 


comes 


from  the  following  theorem;  Let  F(x,y)  be  any  continuously  differentiable 


function  in  D which  vanishes  on  the  boundary  C;  then 

(2  JJVdx  1y)2 

(1.5)  

Jj  ;7F"dxdy 

D 


Thus,  the  stress  function  f(x,y)  solves  an  extremum  problem  for  the 
ratio  on  the  right  side  of  (l.r>)  and  the  torsional  rigid ity  P is  the 
maximum  value  of  t:;is  ratio. 


Let  now  f(x,v)  t®  the  correct  stress  function  of  the  domain  D(t  ). 

o 

Then  f(t  t *x,y)  is  defined  in  D(t).  vanishes  on  its  boundary  C ( t ) and  is  an 
admissible  function  in  the  variational  problem  (1.5/  that  defines  P(t). 


Therefore , 

(1.6'  Pit)  ^ 


\2  ff  f(t  t."  "x,y)dx  dy)2 

D(t)  G 


Jj  Lt^t  "f.  It  t ~x,y)2*f  (t  t *x,y)2lc:xdy 
D(t)  ° - ° 2 ° 


where  and  f.  denote  the  derivatives  of  f with  respect  to  its  first  and 
second  argument,  respectively,  ffe  change  now  our  variables  of  integration 
by  referring  back  to  the  domain  D(t^);  thus,  (i.6)  becomes  after  some  easy 


transformations  Jf  , < „fdx  dy  . ft  f (x,y)2dxdy 

° nc.  ' i * nu  1 2 


(1.7) 


ml6  b 


nit  : 

_a 


D(t  ) 


(2  jj  f(x,y)dxdy)‘ 


D(t  ) 


o 

-2 

Here,  we  put  T - t . Observe  that  the  right  side  of  (1.7)  has  the  form 

aX  +b  and  that  for  T - XT  - t we  have  equality  in  ( 1 . 7)  by  virtue  of 

o o ^ 

the  definition  (1.3). 

Consider  the  curve  with  the  ordinate  tP(l)  and  the  abscissa  T - t 
As  (l.v)  shows  (since  t is  arbitrary),  there  passes  through  each  point  of 
this  curve  a straight  line,  no  point  of  which  lies  below  the  curve.  A curve 
with  thi3  property  is  obviously  a convex  curve  and  the  straight  lines 
considered  are  its  supporting  lir.es. 


We  fin:  from  (1.7)  also  the  3lope  of  the  supporting  line  at  the  po^nt  Tq 

J / f?  dx  dy 
D(t  ) 1 

(1.8)  a - Z~>U  2 5 ^0 

° p(t  r 

o 

Thus,  the  curve  considered  in  nor.-decrea3ing.  We  collect  our  results  in 

Theorem  1.  Let  P-F(t)  be  the  torsional  rigidity  of  the  simply- 

connected  bounded  plane  domain  D-D(t)  which  arises  from  D(l)  by  or.a-sided 

stretching  in  a fixed  direction  with  the  parameter  t.  Then  tP  considered 

-2 

as  a function  of  the  variable  T • t , is  an  increasing  function  and  convex 
from  above. 

1.2.  We  wish  to  add  a few  additional  remarks.  Instead  of  the  expression 
tP  1 we  might  also  consider  the  functional  AP  ^ where  A(t)  is  the  area  of 

the  domain  D(t).  In  fact,  A(t)  is  proportional  to  t and  the  curve  AP  ^ 
with  the  abscissa  C ' t 2 has  the  same  properties  as  the  curve  tP  1. 

If  D(l)  is  a circle  with  unit  radius,  D(t)  is  an  ellipse  with  semi-axes 
1 and  t.  Obviously, 

(1.9)  f(x.y)  • . T-t'2  . 

is  the  stress  function  of  D(t)  and  we  calculate  easily  from  (1.4) 

3 

(1.10)  p(t)  - , tP(t)"1  - ^(H  *1) 

1 * t 

Thus,  in  this  particular  case  the  convex  curve  becomes  simply  a straight 
line  which  shows  that  the  inequality 

(1.11)  ( tP~ 1 ) — 0 
d r 2 

is  the  best  possible  one. 

1.3.  Let  us  return  to  the  general  domain  D(l)$  it  contains  a certain 
circle  of  radius  r and  is  contained  in  some  circle  of  radius  R.  We  know, 
therefore,  that  the  domain  D(t)  contains  an  ellipse  with  semi-axes  r and  rt 


and  is  contained  ir.  another  oliijvse  with  semi-axes  R and  Rt . The  torsional 
rigidity  P is  a mono  tom:  set  function  as  follows  immediately  from  the 
extremum  problem  (1  .t)  which  defines  it.  Hence,  using  (1.10),  we  have 
the  estimates 

(1.12)  (~~R4'  1 v X * 1 ' — tP(  t ) ~ 1 — ( TTr4, ) ^ ( X * 1 ) 


Thus.  tP(t)  ^ as  function  of  X is  enclosed  between  two  straight  lines 
both  of  positive  slope  and  having  a positive  intercept  on  the  vertical 
axis.  In  particular,  we  3ee  that  tP(t)  is  strictly  increasing  with  X 
and  that  the  equality  sign  m vi.3)  is  excluded.  The  slope  a of  the  curve 
does  not  increase  since  the  curve  is  convex  from  above;  but  we  have  now  shown 
that  even  the  limit  of  this  slope  for  X — *oo  must  still  be  positive. 

1.4.  Using  theorems  on  the  continuity  of  the  Green's  function  in 
dependence  on  its  domain  of  definition,  we  can  prove  that  P(t)  is  differentiable 
in  t.  Prom  the  fact  that  tP(t)  * is  an  increasing  convex  function  of  X, 
however,  it  follows  directly  and  in  a nucn  more  elementary  way  that  it  has 
a derivative  everyw.oere  except  for  a null  set  in  X,  Wherever  this  derivative 
exists,  it  is  given  by  (1.8).  From  this  formula,  we  car.  easil7  compute  the 
derivative  of  P(t)  with  respec  t to  t;  ’using  1 1 . *1 ) . we  obtain 


(1.13) 


tfU)  = 1 

dt  t 


ft,  idx  iy 


9 


D(t) 

where  f(x,y)  is  the  stress  function  of  D(t).  Thus,  P(t)  is  a strictly 
increasing  function  of  t. 

We  can  combine  formulas  (1.4)  and  (1.8)  in  order  to  derive  the 
inequality 

(1.14)  ^ [tpur1^  x'^pur1  - ± ttpur1! 


The  right  side  term  represents  the  slope  of  the  straight  line  that  joins 
the  point  considered  on  the  curve  to  the  origin.  This  fact  could  have  been 


9 


deduced  immediately  from  the  convexity  cf  the  curve  without  the  use  of  an 
explicit  expression  for  the  derivative  of  t.P(t)  \ 


1.5  The  torsional  rigidity  P has  been  expressed  by  (1.5)  aa  the 
maximum  value  of  a certain  functional  m a specified  f'unction  class  F(x,y). 
This  fact  leads  to  easy  lower  bounds  for  P and  was  utilized  systematically 
in  the  preceding  sections  ir.  i:de:  to  obtain  convexity  results  for  the 
torsional  rigidity  We  wish  to  formulate  now  a minimum  problem  connected 
with  the  torsional  rigidity  which  will  provide  upper  bounds  for  P and  yield 
complementary  convexity  tr.ecrems  for  this  functional. 

Let  D be  a given  simply-connected  bounded  domain  in  the  (x,y)-plane 
and  q^(x,y)  (1=  1,2/  be  a continuously  differentiable  vector  field  in 
it  such  that  its  divergence  satisfies 


(1.15! 


d B 

d x ql 


Let  f(x,y)  be  the  stress  function  of  D defined  in  Section  1.1.  We  clearly 
have  by  (1.4),  (1.15/  and  the  'ar.isnmg  of  f on  C: 


(1.16) 

Hence 

(1.17) 


ax  ay 


f dx  dy  = P 


(q^  * q2)dxdy 


I (fc  * f^)dxdy 

! * > 


J ^fx”  ql)2  * (fy”  ^2)2'dxdy 


and,  finally,  in  view  of  (1  ..) 
(1.18) 


(q,  * )dx  dy  ^ P 

± t 


where  equality  can  hold  only  for  q*Vf.  Thus,  the  torsional  rigidity  appears 
as  a minimum  for  a functional  defined  for  all  continuously  differentiable 
vector  fields  in  D with  the  constant  divergence  -2.  No  boundary  conditions 
an  imposed  upon  the  vector  field,  ■‘he  extremum  problems  (i  5)  and  (1.18) 
stand  to  each  other  in  the  same  relation  as  do  Dirichlet's  and  Thomson's 
principles  in  potential  theory.  [Cf.  Section  4.4.' 


Let  now  D^t)  be  again  the  family  of  domains  obtained  from  D*D(l)  by 

one-sided  stretching.  If  f(x,v)  is  the  stress  function  for  a domain  D(t  ) 

o 

it  defines  for  this  domain  the  extremum  vector  field  3-  Vf.  We  can 
transplant  this  vector  field  into  D(t.)  and  utilize  it  as  a test  vector 
field  in  the  minimum  problem  ( 1 . 13 ) for  this  new  domain  We  have  to 
transform,  however,  the  vector  field  in  a proper  way  in  order  to  safeguard 
the  side  condition  ( 1 . 15 ) Let  us  define,  therefore,  in  D(t)  the  vector 

field 

t t to 

(1.19)  f-  qn  ("f  x,y)  , q2(-jT  x,y) 

In  view  of  (1.18),  we  then  obtain  the  estimate 

(1.20)  P(t)  ^ J C(^*)2  q,  (-—  x,y)2  ♦ q2(~  x,y)2)dxdy 

D(t)  ° 

and,  referring  back  to  the  domain  D(tQ),  we  find  in  the  usual  way: 

(1.21)  ^ ^ jJ  Jf  i^axdy 

*“0  D(t  ) 0 D(t  ) 

o o 

Thi3  inequality  leads  to 

Theorem  la:  Under  the  assumptions  of  Theorem  1 the  functional 

wl  2 

t xp(t)  is  an  increasing  function  of  the  variable  s-t  and  is  convex 

from  above. 

The  complementary  character  of  Theorems  1 and  la  is  seen  best  if  we 
express  the  convexity  statements  as  differential  inequalities  for  P(t) . 


We  may  express  Theorem  1 in  the  form 


(1.22) 


while  Theorem  la  leads  to 


(1.23) 


3 


El 

t 


11 


He  see  that  the  theorems  provide  upper  and  lower  estimates  for  P ( t ) , 
respectively. 

The  method  applied  in  this  section  may  be  called  the  transplantation 
of  an  extremal  vector  field.  It  is  of  use  in  the  case  of  a Thomson 
principle  for  a given  functional  while  the  method  of  transplanting  an 
extremal  (scalar)  function  is  appropriate  in  the  case  of  a Dirichlet 
principle.  Friedrichs  (5"1  has  studied  Legendre  type  transformations  which 
can  bring  a wide  class  of  minimum  problems  with  prescribed  boundary 
conditions  into  ar.  equivalent  form  involving  a maximum  problem  with 
natural  boundary  conditions . By  such  transformations  one  can  frequently 
obtain  for  the  same  functional  a Thomson  as  well  as  a Dirichlet  type 
extremum  characterization.  The  duality  between  the  two  extremum  problems 
is  reflected  in  the  duality  between  the  two  transplantation  methods. 

In  transplanting  a vector  field  we  have  to  change  the  components  of 
the  field  as  well  as  the  independent  variables  of  the  field  as  is  apparent 
in  (1.19).  The  transplantation  of  such  a field  under  general  linear  trans- 
formations of  the  independent  variables  will  be  discussed  in  Chapter  TV. 
Various  further  applications  of  the  method  will  be  given  in  the  following 
sections  of  thi3  chapter 

2.1.  Let  B be  a finite  body  ir.  the  (x,y ,z)-space  bounded  by  a smooth 
surface  S.  Suppose  that  B is  immersed  in  an  incompressible  fluid  which 
flows  irrotationally  past  the  body  under  the  influence  of  a dipole  at 
Infinity  of  strength  1 and  directed  along  the  x-axis.  The  velocity 
potential  of  such  a flow  has  the  form 

(2.1)  ^(x,y,z)  = x f Cp(x,y,z) 


regular  namcnic  outride  B We  have  on  5 the  boundary 


where  Cp(x,v,z'  is 
condition 

u-)  ’ 


"'zr*-  = - cos  (n,x) 
dr. 


which  determines  <p(x.y,z)  in  a unique  way. 


(2.3) 


Let  D be  the  exterior  of  B;  we  define  [15, l?"1 

rrp 

*x  ■ |//1  vqpi‘dx  dydz 


as  the  virtual  mass  of  B m the  x-dlrection.  W is  a functional  of  3 

x 

which  plays  a role  in  many  hydrodynamics!  applications. 

Consider  a function  f(x,y,z)  which  vanishes  at  infinity,  is  continuously 
differentiable  in  D and  has  a finite  Lirichlet  integral  there.  Clearly 

(2.4)  ^ Jvf'^Cpdxdynz  * - f d 6~  * jT  f cos  (n,x)d  <T 

n 


Thus,  by  Schwarz1  inequality 


(2.5)  [jlJ(Vf!;dxd7d,.»x  3 >f[  f cos  (n,x)d(T  ) 


and  the  equality  sign  can  held  only  for  Vf  = c V^P,  where  c is  any  constant. 
We  have,  therefore,  the  folic* ing  definition  of  the  virtual  mass  by  means  of 
an  extremum  problem:  pp  ? 


f cos  (n,x)d  CT ) 


(2.6) 


W = max 
x 


u 


r ~ 2 

I : 1 Vi  ) cix  dy  az 

V 


2.2  If  we  introduce  tne  transformation  x'  - tx,  y'»y,  z/-z,  we 
obtain  again  a family  of  bodies  Bit'  with  exteriors  D(t),  boundary  surfaces 
3(t),  and  virtual  masses  in  the  xdirection  W (t) . In  order  to  study  the 
dependence  of  W (t)  cn  t,  we  use  once  more  the  extremum  definition  (2.6) 
of  the  virtual  mass  and  the  transplantation  of  the  correct  extremum  function 
of  a domain  D(t  ) into  a domain  Bit)  as  a competing  function. 
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Let  cp(x.y.z)  be  tne  correct  extremal  function  for  some  fixed  D(to)  and 
consider  f(x,y,z)  - Cp(to*.  *x,y,z)  in  the  domain  D(t).  Clearly,  this 
function  has  a finite  Lirichlet  integral  and  may  serve  to  estimate  W^(t) 
by  means  of  the  inequality  (2.6).  Wo  have 

' fi  ^P^'o*'  "x,}'.z)cos  (n,x)d<T 
U.l)  Wx(t'1  5-  2 -2  _ 2 . _2  . " " 


J f t*  t"?  9^  +(p^dxdy  ds 


We  change  our  variables  in  order  to  use  D(t  ) a.;  domain  of  integration  instead 

o 

of  D(t)j  this  means  simply  a replacement  of  t t ^x  by  x.  Observe  now  that 

o 

cos  (n,x)dCT  ■ dydz  ir>  not  affected  by  this  change  of  variable.  Hence,  we 


(2.8)  [t*  ( t ) ~ 

x 


‘cjff  ‘7>ldxdydl‘  jjj  (Pjldxdyd! 

!>(t  } D(t  ) 


V if  <p  cos  (n,x)d  <T  )2 
S(t  ) 

O 

using  again  the  variable  T.  » t . Thus,  we  have  an  estimate 


(2.9) 


ftW  (t0]'x  ^ M> 
x 


— y . 

and,  for  T7 » a t , we  have  by  (2.6) 


(2.9') 


[tt  (t  ) T’-*'  at  *b 

oxo  c 


_ i _2 

This  result  shows  that  tr.e  curve  [tW^(t)]  with  'Z  -t  as  independent 
variable  has  at  every  point  a supporting  line  and  is  convex  from  above. 

Thus,  it  has  a tangent  everywhere,  except  for  a null  set;  wherever  the 
tangent  exists,  its  slope  is  given  by 


(2.10) 


W (t) 

x 


2 j jj  qpi^dydz^  0 


This  follows  from  (2.8)  and  the  equation 
(2.11)  Jf  Cp  cos  (n,x)  d (T  = 


S(t  ) 
o 


which  i3 , in  turn,  an  immediate  consequence  of  (<-.2)  and  (2.3). 

We  conclude  from  (2.10)  that  tW  (t)  is  a non-decreasing  function  of  t. 

x 

Since  the  vnlume  V(t)  of  the  body  B(t)  is  given  by  V(t)  ■»  tV(l)  we  may  also 

say  that  the  combination  VI  is  non-decreasing  in  t. 

x 

We  summarize  our  results  in 

Theorem  2.  Let  W^(t)  be  the  virtual  mass  in  the  x-direction  of  a 
body  B(t)  which  arises  from  B(  1 ) by  one-sided  stretching  in  the  x-direction 
with  the  parameter  t.  Then  [tl^(t) J ^ is,  considered  as  a function  of  the 
variable  't-  * t , an  increasing  function  and  convex  from  above. 

2.3 . The  curve  [tW^(t)l  ^ with  abscissa  “V  behaves  asymptotically 
quite  differently  from  the  curve  tP(t)  ^ considered  in  the  preceding 
paragraph.  In  fact,  it  is  easily  seen  from  (2.10)  that  the  slope  of  the 
curve  converges  to  zero  if  T — ^ao,  that  is,  t— >0.  This  is  due  to  the 
fact  that  even  a plane  plate  has  a finite  virtual  mass  if  it  stands  under 
right  angles  to  the  flow  direction;  and  for  t — >0  the  function  W^(t)  will 
converge  to  the  virtual  mass  of  the  plane  plate  obtained  by  projecting  the 
body  B(l)  into  the  (y,z)-plane. 

The  plane  plates,  perpendicular  to  the  x-direction,  give  an  elementary 
family  of  bodies  B(t)  for  which  the  explicit  dependence  on  X.  of  (tW^(t)^ 
can  actually  be  computed,  iince  3(t)~  B(l),  we  see  that  in  this  particular 
case 

(2.12)  (tW^t)"!  ^ , k independent  of  X , 

which  is  a monotonic  curve  convex  from  above. 

2 .A . We  can  derive  further  results  for  the  virtual  mass  by  using  the 
method  of  transplanting  an  extremal  vector  field.  For  this  purpose,  we  have 
to  characterize  W^  by  a Thomson  type  extremum  principle.  We  define  in  the 
domain  D a solenoidal  vector  field  q-q^^.q^,  that  is,  we  require  the 
differential  condition 
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(2.13)  7q  - 0 

We  suppose  further  that  |q|*0("^)  as  r — ^ oo  and  that 


U.u) 


If 


q dx  dy  dz  <.  ac 


ana  that  the  boundary  condition 
(2.15)  q-fl  " -cos  (n,x) 
be  satisfied  on  d . 

The  velocity  potential  2p(x, y,z)  defined  in  (2.1)  leads  to  the  vector 
field  V Cp  of  the  required  type.  If  q is  another  admissible  field,  we  have 


(2.16) 


iff 


q dx  dy  dz 


Jj  ( V93  )^dxdy  dz  * 2 J(q  - VCp  ) VCf>dxdy  dz 


Using  the  harmonic  character  of  Cp, 
we  find 


' pp 

♦ (V cp-  q 

JdJ 

D 

the  conditions 


)^dxdydz 

(2.2),  (2.13),  and  (2.15), 


(2.17) 


(q  - VCp  )-Vcpdxdy  dz  - 0 


thus  (2.16)  leads  to  the  inequality 


(2.18) 


r r r 2 

q"dx  dy  dz  ^ 


JJJ 


J^dx  dy  dz 


and  equality  can  hold  only  for  q"  V<p.  By  the  definition  (2.3)  of  the 


virtual  mass,  we  have  therefore 

2 


(2.19) 


Wx  - mm 


r 

k 

sj 


dx  Jv  dz 


where  ail  solenoidal  vector  fields  with  the  boundary  condition  (2.15)  and 
the  proper  behavior  at  infinity  are  admitted.  This  minimum  condition 
characterizes  the  velocity  potential  and  yields  at  the  same  time  estimates 
for  the  virtual  ma3s  W^;  it  is  called  'he  Thomson  principle  for  the  virtual 
mass  problem.  It  leads  to  easy  upjer  bounds  for  W , in  contradistinction 
to  the  extremum  problem  (2,6)  nhich  permits  estimates  for  W^  from  below. 
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2^5  We  consider  now  the  one-sided  stretching  of  the  body  B-B(l)  in 
the  x-direction.  Let  cp  be  the  correct  velocity  potential  in  the  exterior 
D(tQ)  of  the  body  B(tQ).  We  transplant  the  vector  field  q E V Cp  into 

the  domain  D(t)  by  putting 

/ ,t  to  t \ 

(2.20)  qt  x>y'^  * Q2  x’y>z)  ’ q3^"?  x’y>2^  J 

Clearly,  is  well-defined  In  D(t)  and  solenoidal  there.  Its  norm  is 
incegrable  over  D(t)  and  it  is  easily  verified  that  on  the  boundary  S(t) 
we  have 


(2.21) 


qt-a 


.L 


cos 


3 (n,x) 


t 

Thus,  we  have  to  use  the  vector  field  q^  as  test  field  in  (2.19) 


We  find  in  the  usual  way 


f (*) 


(2.22) 


xC  JL. 

~ t 


. qidldyd2‘  $ JR 

D(t  ) t D(t  ) 

c o 


(q^  * q^)dx  dy  dz 


This  inequality  shows  that  t ^W  (t)  is  a convex  (from  above)  function  of 


_ 2 

^ - t ; the  slope  of  the  supporting  line  at  every  point  t is  given  by 


(2.23)  a - t 


ff/V 


tj  (q^  * q7  )dx  dy  dz  - t j ML^fy  * ^ z Idxdyd 


D(t) 


D(t) 


We  proved  thus 

Theorem  2a.  Under  the  assumptions  of  Theorem  2,  the  functional 

t hi  (t)  is  a non-decreas i nfy  function  of  £ _t  ^ which  is  convex  from 
x 

above . 

We  may  formulate  Theorems  2 and  2a  in  the  form  of  the  differential 
inequalities  ^ 

(2.24)  f- 

t X 

and 

(2.25) 

V 
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We  see  again  that  the  two  results  lead  to  upper  and  loner  estimates  for 
We  can  combine  both  inequalities  and  obtain 


(2.26) 


(w^M2  ^ i/t2 


3.  The  Outer  Radius. 

3.2  Let  Z be  a smooth  closed  curve  in  the  complex  z-plane;  let  D be  the 
exterior  of  C There  exists  n unique  univalent  analytic  function  in  D which 
maps  D onto  th°  exterior  of  the  unit  circle  and  which  has,  near  infinity, 


the  series  development 


(3.1) 


1 al 

f ( z ) - "^  z * a ♦ — ♦ . . 


, r 7>  0 


The  coefficient  r(C) * r,  a functional  of  the  curve  C (or  its  exterior  D), 
is  called  the  outer  radius  of  C. 


(3.2)  g(z)  - log  | f(z)  I - log  | zl  + log  ♦ O(-j-^-r) 

is  the  Green's  function  of  D with  the  logarithmic  pole  at  infinity,  we 
may  define  also  the  outer  radius  potential-theoretically  by 


(3  .2  / ) log  ■A-  * iia  (g(z ) - log  ! z | ) 

r lzj-*CD 

If  z-z(s)  is  the  parametric  representation  of  the  curve  C in  terms 
of  its  length  parameter  s,  we  define  for  z(s)tC 


(3.3) 


_L_  dgiz) 
2TT  jn 


where  jj  is  the  interior  normal  with  respect  to  D.  We  clearly  have 


(3.4) 


!s)ds  « 1 


It  is  also  easy  to  derive  from  (3.3)  and  Green's  identity  the  formula 


^ 7 i Vis  )ds  - log  ^ - g(£  ) 
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This  result  holds  still  for  ^,£C  and  we  find 


(3.6^  log 


1. 


log 


izr  z2' 


V(s^ ) V(a^  )ds^ds2 


C t 

with  z^z^),  z -z(s2). 

Let  next  p(s)  be  an  arbitrary  continuous  function  with 

(3 .”)  J'  p(s )ds  3 0 

C 

Consider  the  integral 

A 

(3.8)  "£}('*' , y;  ) - ^ log  <^(s )ds  , £ - C ♦ i q , 

c ’ * 

which  represents  a harmonic  function  hg( ^ , I'j  ) in  D and  another  harmonic 
function  h^(  V , ) in  the  interior  D of  C.  The  two  harmonic  functions 
have  the  same  limit  values  at  ^ach  point  ^ € C , but  their  normal  derivatives 
jump  by  the  amount  2*TT  p(s  ) if  we  cross  the  curve  C at  the  point  z(s). 
Green's  identity,  we  have  the  relation 


(3.9) 


(Vhe)2d^dr|  *jj  (Vhi)2d^drj) 

n C 


ah  ah, 

- 1 h t ~rr  - ^~'d. 


e d n 3n 


('  r 


Z',TJ  J l0g  Tcfcij  PV  Pls2,dsld32  ? 0 


c c 


'1  2 


Finally,  we  observe  that  by  (3.5)  and  (3.7) 


(3.10) 


log 


; ^ — ~ Vis.  ) P(s0)ds,d3-  - 0 

z^-  z^l  1 \ 2 12 


Let  now  p(s)  be  continuous  on  C and  satisfy 


(3.11) 


Jf 


ds  ■ 1 


Then,  using  p(s)*'V(s)-  p(s),  we  easily  see  that  in  view  of  (3.6),  (3.9), 
and  (3.1C): 


(3.12)  log  -Jr  ^ 


log 


C C 


zrV 


p>( s ^ ) p(3-, )ds„ds. 


,2/ual'ia2 
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Thus.  the  outer  radius  can  be  defined  by  an  extremum.  problem  with  respect 
to  continuous  functions  on  C viith  integral  1. 

The  extremum  problem  defining  log  , already  considered  by  Gauss, 
has  a simple  electrostatic  interpretation.  In  fact,  consider  a charge 
distribution  ^jl(s)  with  logarithmic  potential,  total  amount  1 and  spread 
over  the  curve  C.  The  integral  in  (3.12)  repre3<nts  the  total  energy  of 
this  charge  distribution.  The  equilibrium  charge  distribution  "V(s)  is 
distinguished  by  the  fact  that  it  minimizes  this  energy,  ^he  minimum  energy 
log  may  be  called  the  capacity  of  the  conductor  C. 

3.2  We  consider  now  a family  of  curves  C(t)  which  are  obtained  from 
one  fixed  curve  C(l)  by  one-3ided  stretching  with  parameter  t.  Let  C(tQ) 
be  one  fixed  curve  and  'i^s)  its  corresponding  charge  density.  If  C(t)  is 
any  other  curve  of  the  family  we  define  on  it  a density  by 
(3.13)  yi^sOds'  -*  V(s)ds  , 

where  s'  is  the  corresponding  length  parameter  on  C(t).  In  view  of  (3-12), 
we  have 


(3.14) 


C(t)  C(t) 


log 


yi(s^)  ^(s'Ms^ds^  . 


Referring  back  to  the  curve  C(to),  we  find  easily 


(3.15)  106  ~ / f 


r(t)‘ 

We  may  interpret 


log 


°(to)  C(t0)  ^ (xx- x2)2  * (yj- j7Y 

o 


■V(s^)'V(s2)ds1ds2  . 


expt  f f log  — ~T^~2 

c(t  ) eft  ) 1 *to<yi- 


O O 


- G[ 


t2  *o 


t2(xrx2)2*to(yry2)2 


(3.16) 


J/(s1)  J^(s2)ds1ds2l 
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as  a geometric  mean  with  the  nor.- negative  weight  function  ~p( s);  thus,  we 


car.  write- 


(3.1?) 


2 2 

t t 


>2  -GI  2 


y 1 

2,  l2  , . 2 , -.2 

^ (xj - * ) t (y  - y2) 


rte  apply  this  result  for  two  different  values  of  t,  say  for  t and  T. 

n 

Using  the  well-known  Holder-Minkowski  inequality  [?f  p.  21,  (2.7.1) 1 


(3.18) 


0^  a 1 * G T b ^ trrG[a  - b 1 


we  find 


.2  t2  t2  T~  t~ 

(3  19)  — — ♦ — r;r — ♦ -2 

(+)2  It)2  " \2,  \2  t.2,  x2  _2,  2,  y 

r(t)  r(T ) t (x1-x2)  t0(y1_y2'  T U^^)  *t  (y^y^ 


T2  t2 


v2  * ~2 , 


Now,  we  utilize  the  fact  that  if  A 0'*0,  B 2>0,  and  w > 0 

(3.20)  = “ ~~^AB'  o 

dw*  A"*h  (Aw*B)3 

This  shows  that  w/(Aw*5)  is  convex  from  above  and  satisfies  the  subadditivity 


condition 


(3.21) 


1 r_!i 

2 Aw:  ♦ Aw2  ♦ B A.|  (Wi  ♦ „2)  * B 


Let  us  suppose  that  t,  t , and  T are  choser.  so  that 


(3.22) 


t2  - \ (t*  ♦ T2) 
o <■ 


Then  we  find  from  (3. 19)  and  (3.2l): 


(3.23) 


1 t-^  • £ ot 

2 r(t)2  r (T ) 


2 2 
(x1-  x2)  ♦ (y1-  y2) 


3ut  2>( s)  was  just  the  equilibrium  charge  distribution  of  the  curve  C(t  ) 

o 

and  hence  by  (3.6)  the  right  side  of  (3.23)  is  just  t“V(r(t  )*).  ^hus, 

oo 


we  finally  proved: 


(3.24) 


o 2 ^ (t2  ♦ T2) 

[-t—  + _xL-;  ^ 

r (t)2  r(T)2  r(J  (t2*T2)) 


^ ^ 2 
that  is,  t /(r(t)^)  is  convex  from  above  as  function  of  t . 
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can  transform  this  result  by  observing  the  identity 

(3.25)  ^ f(x)  - [xf(x)1  , 9 - l/x 

dp  dx  J 

which  holds  for  every  twice  differentiable  function  f(x).  This  formula 
shows  that  from  the  convexity  of  xf(x'  in  x > 0 follows  the  convexity  of 

fix)  in  dependence  of  - l/x  > G.  ^hus,  (3.24)  leads  to  the  result-  that 

2 2 
r ( t)  ‘ is  convex  from  above  as  f '.unction  of  'L  ■ t 

_2 

Finally,  it  is  easily  seen  from  (3.14)  that  r(t)  does  not  decrease 

with  increasing  T. . T^.us  , we  may  state 

Theorem  3-  Let  C(t;  be  the  smooth  closed  plane  curve  which  arises 

from  C(l‘  by  one-sided  stretching  in  a fixed  direction  with  the  parameter  t; 

let  r(t)  be  its  outer  radius.  Then  r(t)  * is  non-decreasing  and  convex 

_2 

from  above  as  a function  of  the  variable  't  - t 

3.3.  It  is  easy  to  extend  Theorem  3 to  arbitrary  continue  in  the  plane 
by  approximating  them  by  smooth  closed  curves  for  which  the  result  has  been 
established.  There  exist  two  continue  for  which  r(t)  ^ is  a linear  function 
of  'C  . In  fact,  consider  any  line  segment  of  length  £ ; it  is  easily  seen 
that  its  outer  radius  is  r 3 ^/u.  If  the  segment  lies  in  the  direction 
of  the  stretching,  we  have 

(3.26)  r ( t )"*”  - *^4  C 

r 

while  if  it  lies  in  normal  direc  tion  r('t)r<"  will  be  unaffected  by  the 

stretching.  Thus,  in  the  inequality 

,2 

(3.27)  r(t)  " 0 

d r" 

equality  does  hold  for  these  two  types  of  segments. 

If  C(l)  is  a circle  of  radius  R the  curves  C(t)  will  be  ellipses 
with  semi- axes  R and  Rt  It,  is  easy  to  compute 
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(3.28)  r(t)‘'-'^b(i*t)^-‘^  — P 

R (aTT  * 1) 

_2 

This  formula  enables  us  to  determine  the  asymptotic  behavior  of  r(t) 
in  the  case  that  C"  -">0  or  oo  for  many  curves  C(t).  In  fact,  let  us 
suppose  that  C(l)  contains  a circle  of  radius  Rq  in  its  interior  and  is 
enclosed  by  another  circle  R^.  The  outer  radius  r is  a monotonic  3et 
function  as  can  be  seen  easily  from  (3.2  ) and  the  maximum  principle. 
Thus,  by  means  of  (3.28),  we  have  the  estimates 


(3-29)  (1  * tT2  r(t)"2  \ (1  * t)'2 

RT  R 

1 o 

3 . L . We  defined  the  capacity  log  l/r  of  the  curve  C as  the  minimum 
energy  of  a charge  distribution  of  total  amount  1 spread  over  C.  Suppose 
now  that  we  bring  n charged  particles  on  the  conductor  C,  each  having  the 
charge  l/n;  what  is  the  minimum  energy  of  this  set?  Clearly,  we  have  to 
solve  the  minimum  problem: 

(3.30)  “^7  log  77  -min  , i,k-l n 

2n"  i/k  l?i'Zkl 


th 

if  z^£C  is  the  location  of  the  i ' particle  in  the  equilibrium  position. 

Closely  related  to  this  problem  is  the  question  for  the  n*^  order 

diameter  of  a curve  C which  i3  defined  by  [3"1 

(n) 

(3 .31)  (d  (C) ] ^ - max 

n 


7T 

i<k 


z.  - z.  ; 
1 k*  ’ 


i,k-  l,...,n;zi<eC,  zk£C 


It  can  easily  be  shown  that  the  d (C)  form  a montonically  decreasing 

sequence  and  that  they  converge  to  the  outer  radius  r of  C. 

Ie  can  deal  with  the  d in  the  same  way  as  we  dealt  with  the  outer 

n 

radius.  Transplanting  the  extremum  points  z^  under  a stretching,  we  can 

_2  -2 
prove  that  d is  a non-decreasing  convex  function  of  f * t 


I 
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We  remark,  finally,  that  the  can  be  defined  for  an  arbitrary  plane 

point  set  and  that  lin  d does  always  exi3t.  This  limit  i3  called  the 

n— >cn 

transfinite  diameter  d of  the  sei.;[3'1-  It  can  be  shown  from  our  preceding 
result  that  d is  a non-decreasing  function  of  'T  ■ t and  convex  from 
above . 


4 • Electrostatic  Ca ? acity . 

4.1-  Let  5 be  a smooth  closed  surface  in  the  (x,y , z )-space  and  let  D 
be  its  exterior.  We  define  in  D a harmonic  function  V(x,y,z)  which  satisfies 
on  3 the  boundary  condition 
(4.l)  V = 1 on  i 

V is  called  the  conductor  potential  of  the  surface  S;  it  has  at  infinity 
the  form 


(4.10 


V ( x , y , 2 ) * ~ ♦ 0(^1 ) 


2 2 2 

x ♦ y ♦ z 


and  the  factor  K is  called  the  capacity  of  3.  K may  be  interpreted  as 
the  electric  charge  which  induces  on  the  conductor  3 the  potential  1. 
By  Green's  theorem,  we  hav« 

(4.2) 


{ VV)^ dx  dy  dz 


no  a y 

V ~~  d C“  ■ 4TTK 


an 


If  f(x,y,z)  is  an  arbitrary  function  which  is  continuously  differentiable 
in  D,  has  a finite  Dirichlet  integral  over  D and  vanishes  on  3 and  at 
infinity,  we  have 

(4.3)  ^^|'v7V-  Vf  dx  dydz  =*  0 


Thus,  clearly  by  the  usual  arguments 

J(Vu)^dx  dy  dz 

) 

where  all  continuously  differentiable  functions  U(x,y,z)  with  finite 
Dirichlet  integral  are  admitted  which  have  on  S the  boundary  value  1 and 


(4.4) 


min  — 
47T  J J 
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which  vanish  at  infinity.  This  minimum  principle  characterizes  at  the  sane 
time  the  conductor  potential  V and  the  capacity  K.  It  is  called  the  Dlrichlet 
principle  for  the  conductor  problem. 

If  S(l)  is  a fixed  surface,  we  may  obtain  the  family  3(t)  by  one-sided 
stretching,  say  in  the  x- direction,  and  denote  by  K(t)  the  corresponding 
value  of  the  capacity.  From,  the  definition  of  K by  means  of  the  extremum 
problem  (4-4).  we  can  derive  the  following  results  by  the  method  of  trans- 
planting the  extremal  function: 


(4.5) 


U&l 

t 


,2. 


T~^  " ( j 1 [~r  V2  ♦ V2  * V2  Idx  dy  dz 
4rrt  J j 2 : y z J 

0 D(t  ) 1 


where  V(x,y,z)  is  the  conductor  potential  of  the  surface  S(tQ).  Thus, 
we  have 

Theorem  4-  Let  3(t)  be  a closed  smooth  surface  which  arises  from 
3(l)  by  one-sided  stretching  in  a fixed  direction  with  the  parameter  t; 
let  K(t)  be  its  capacity.  Then  K(t)t  ^ is  a non-decreasing  function  of 
■ t and  convex  from  above. 

The  slope  of  the  supporting  line  at  every  point  t of  the  curve  is 
given  by 


(4.6) 


4~ 


V2 


dx  dy  dz 


D(t) 


From  this  formula  we  can  easily  compute 
U.7)  £ K(t)  - ^ 

D(t) 


can  define  on  the  surface  S the  charge  density 
’ " 4^tK  ~?~n 

which  is  by  the  maximum  principle  non-negative  on  S and  satisfies  the 
condition 

(4.9)  If  V(P)d<T  - 1 

S 


a 
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If  F and  Q are  points  of  D we  shall  denote  their  distance  by  r(P»Q)- 


(4.10) 


V(Q)  - ff  v(P)r(P,Q)"1drrp 


and  by  virtue  of  (>4.9): 


(4.11) 


K "jf  J VlP)  V(Q)r(p,Q)'ld  CjA 


It  is  also  easy  to  verify  Saus3:  theorem: 


(4.12) 


min  \^\  J u(P)  |-UQ)r(P,Q)  1d  (7^>d 


where  all  functions  ^x(P)  on  S are  admitted  which  satisfy 


(4.13) 


f v(?)d(T  - 1 


He  obtain  from  the  extremum  definition  (4.12)  for  the  capacity  an 

estimate  for  the  quantity  K(t).  Transplanting  the  correct  charge  distribution 

V(P)  of  S(t  ) onto  3(t)  by  the  law 
o * 


(4.14) 


u(P)dr7  - V(P)d<T 


where  diT'  and  dtT  are  corresponding  surface  elements,  we  obtain: 

(4*15)  KUl-ff  fj  )2*  (y-  7)2*(z-£  )2rl/2^(PMQWpd<rQ 


S(t  ) d(t  ) "o 
o o 

We  remark  next  that  if  a,  b,  and  t are  positive 


(4.16) 


>t  ♦ b 


< 0 


Hence,  we  obtain  from  (4.15): 


(y-  r?  )2  ♦ (z  - 'Z  )2 


-17>  iao-Krtr*(t-to)  II  Jf  ~ 

o c i c/.  i rlP, 


5 (t  ) S(t  ) 
o o 


MP)r(Q)dcTpd^ 


This  estimate  proves  that  fcK(t)  i3  a non-decreasing  function  of  t 
which  is  convex  from  above.  This  can  also  be  formulated  as  the  statement 
that  K(t)  ^ is  convex  from  above  in  dependence  on  t (See  (3.25).)  If  S 
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is  n plane  surface  perpendicular  to  the  x-direction,  we  have  3(t)  S(l) 
and  in  this  case  tK(t)  becomes  a linear  function  of  t while  K(t)  ^ is 


independent  of  t 


_ i 


Thu.,  our  result  is,  in  a certain  sense,  a sharp  one. 


By  (4.17),  we  find  as  the  slope  of  the  supporting  line  to  the  curve 


tK(t)  versus  t 
(4.18)  a 


) 

c 

.S  v 
Oit 

-1  +> 

1 

1 

o ' o 

We  can  compute  from  this  result 


(v  - ri  ) * (z  - )‘ 


•(P,Q)' 


■v(P)  MQ)d(TpdcrQ 


(4.19) 


ft  K(t) 


0 ^ 


o(t  ) S(t  ) r 
o o 


v(?)  V(u) 
(P,Q)3 


(x  - ^ ) d^d^ 


This  formula  shows  that  the  capacity  grows  with  a one-sided  stretching. 

Comparing  formulas  (4.?)  with  (4-19)  we  obtain  an  interesting  integral 
identity.  It  leads,  in  particular,  to  the  inequality 


(4.20) 


1 


V dx dy  dz  ^ 
x J 


l V*-  ♦ V2  Idx  dy  dz 
y z 


Equality  holds  here  if  the  surface  3 lies  entirely  in  the  (y,z)-plane,  since 
in  thi3  case  we  have  obvious  ly  K ' (t)  • 0. 

In  order  to  compare  the  ccnvexity  results  of  the  last  two  sections, 
we  formulate  them  as  differential  inequalities  for  K(t).  The  convexity  of 
K(t)t  ^ in  dependence  on  t * leads  to  the  inequality: 

(4.21)  t2  & . t # £ K 


dt 


dt 


while  the  convexity  of  tK(t)  in  dependence  on  t yields: 


(4.22) 


(t^*2  ^)K  2tK'2 

dt2  dt 


d*K 


We  see  that  the  two  inequalities  lead  to  estimates  of  14_2  from  both  sides 

dt 

and  thus  complement  each  other.  This  was  to  be  expected  since  they  were 


derived  from  the  complementary  extremum  problems  (4-4)  and  (4.12)  which 
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pcrr.it  upjier  and  lower  estimates  for  the  capacity. 


4»J.  The  extremum  problem  (4.4)  for  the  capacity  K can  be  generalized 
by  admitting  continuous  functions  U[x,y,z)  in  D which  are  continuously 
differentiable  except  for  a smooth  surface  in  D across  which  their 
derivatives  are  permitted  to  jump.  In  fact,  the  equation  (4-3)  hold3 
for  functions  f(x,y,z)  of  this  class  and  this  implies  (4.4) 

•^he  fact  that  we  can  relax  the  continuity  requirements  for  the  competing 
functions  allows  us  to  introduce  a type  of  one-sided  stretching  slightly 
more  general  than  that  used  till  now.  We  consider  the  transformation  of 
the  (x,y ,z)-space 


tx  for  x £ 0 


(4.3) 


7'  - y 


x 


for  x 0 

which  is  one-to-one  and  continuous  but  which  has  discontinuous  derivatives 

across  the  plane  x-0  for  t/l.  We  call  this  transformation  a partial 

one-sided  stretching  in  the  x-direction  with  parameter  t. 

Let  B-B(l)  be  a given  body  with  surface  S-S(l),  exterior  D»D(l), 

and  capacity  K-K(l).  By  B(t).  5(t),  D(t),  and  K(t),  we  denote  the 

corresponding  concepts  after  a partial  one-sided  stretching  in  the  x-direction 

with  parameter  t.  Let  V(x,y,z)  denote  the  correct  conductor  potential  for 

the  body  B(t  );  then  the  function 
° ^ t 


(4.24) 


V("^  x,y,z)  for  x^O 
V(x,y,z)  for  x 4:  0 


will  be  wel’.-def ined  in  L(0  and  an  admissible  test  function  for  the  minimum 
problem  (4.4).  Hence 
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U.  25)  i~TK(t)  4=  J j'  [(~^)  ^x^~t  x,y>z)  + ^y*  V^"!dxdy  d: 

2Ht)  ,x  ^0 


JJJ  ;^*'2‘^'dxdydr 

D(t),x£0 


We  refer  back  to  D(t  ) and  denote  the  part  of  D(t  ) n the  half  space 

o o 

- 

x ^ 3 by  D and  the  other  part  by  D . We  find 
J o * J 0 

(4.26)  4TT K(t)—  JJ[ Idx  dy  dz  *JJJ ‘~t*  * t~ 

D“  D* 

o o 

with  V-V(x,y,z).  We  rearrange  the  terms  and  observe  that 


dz 


(4.27)  4^K(t  ) - 

o 


f fT 

jjj, 

D(t0) 


[ 7^  ♦ ]dx  dy  dz 

x y z J 


thus,  we  obtain 


t-  t 


(4.28)  K(t)-K(toK  JFfJIJ^^-^Tdxdy 


(t-  to)"  rr 

4TTt  t 


dz 


dx  dy  dz 
x J 


We  U3e  the  well-known  fact  that  K(t)  has  arbitrarily  many  derivatives  with 
respect  to  t.  Hence,  (4.28)  implies  immediately 

u-29)  ■ 


• ♦ 


Thi3  formula  holds  even  when  the  whole  body  B lies  in  the  half-space  x <.  0; 
clearly,  in  thiis  case,  the  bodies  B(t)  will  coincide  with  B and  hence 
Ky(t)  = 3.  This  leads  to  the  integral  identity 
(4.30) 


f fftV2  ♦ V2  - V2  Idx  dy  dz  - 0 

JJJ  y z x 


if  H is  any  half-space  in  D bounded  by  a plane  x-  const.  This  identity  can 
also  be  verified  directly  by  integration  by  parts. 
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We  obtain  further  from  (4.23)  the  estimate  for  the  second  derivative 
of  K(t): 

(4.31)  K"(t  )£  — ^ I f 

0 2-n"t"  Jj. 

° D 

o 

Finally,  combining  (4.29)  and  (4.31),  we  find 

(4.32)  t*V'(t  ) * t K'(t  )^  -1-  ffftV2*  V2*  V2  dy  dz  ^ K(t  ) . 

o o o o 4TrJJJ  x y 2 3 o 

♦ 

n 

0 

Thus,  we  arrive  at  the  differential  inequality 
(4.320  t2K  "(t)  ‘ tK'(t)  ^K(t) 

which  generalises  (4.2l)  to  the  case  of  partial  stretching.  We  find  thus: 
Theorem  4a:  The  functional  t ^K(t)  is  convex  from  above  in  dependence 

of  the  variable  T = t 2 in  the  case  of  a partial  one-sided  stretching  with 
the  parameter  t. 


4.4.  We  formulate  now  the  Thomson  principle  for  the  capacity  of  a 


conductor.  We  observe  that  V V represents  a vector  field  in  D 

with  the  following  properties: 


(4.31) 

(4.34) 

(4.35) 


V-  a 3 0 

rC(  2 

jjjq  dx  dydz  <4  cc 
D 


-2, 


q | ■ 0(r  ) as  r 


ao. 


II 2' 


fid  <r  - 1 


If  q(x,y,z)  is  an  arbitrary  vector  field  with  continuously  differentiable 
components  in  D and  satisfying  the  above  three  conditions,  we  have 


(4.36) 


rrr 

jq^dxdydz-  ]|jq^dxdydz*  2jjjq^‘ {q- q^)dxdy  dz  + jJJ(q- q^^dxdy  dz 


Y 


We  observe  next  that  by  (4.33)  and  (4.35)  and  the  fact  V - 1 on  S 

qdX  dydZ  " 4^"  if  "s’  V(9’n'd0~  = 


Id 


and,  in  particular, 


(4.3”') 


J JJq0  dxdyd" " ZtTk 

D 

Thus,  (4.36)  leads  to 


(4.38) 


fjjq^dxdyd;  ? ^ 

D 

and  can  be  defined  as  the  minimum  of  the  left-side  integral  for  all 

adaissible  vector  fields  q.  This  is  the  required  Thomson  principle. 

It  is  important  to  observe  that  the  vector  field  q nay  even  be  discon- 
tinuous along  a smooth  surface  ^ in  D,  provided  only  that  (q-fl)  be  the  sane 
on  both  sides  of  In  fact,  this  latter  condition  is  sufficient  in  order 

to  establish  (4.3°)  which  is  equivalent  to  Thomson's  principle. 

Consider  now  again  the  bodies  B(t)  obtained  fron  B»B(l)  by  the 
partial  one-sided  stretching  (4.23).  Let  q*  VV  be  the  extremal  vector 

field  for  the  exterior  D(t  ) of  the  body  B(t  ).  We  define  in  D(t)  the 

o o 


transplanted  vector  field 

+ 

o 


z 

(4.39)  qt  =(q1(Y  x,y,z)  ♦ x>y>z^ 

if  x $0,  while  q^  * q if  x ^ 0.  The  vector  field  q^  is  well  defined  in 
D(t) ; it  is  d isccntinuous  across  the  plane  x*0  but  its  normal  component 
passes  continuously  acrcs3 . It  satisfies  the  conditions  (4.33)-  (4.35) 
with  respect  to  the  domain  D(t)  and  is,  therefore,  a test  field  in  the 
inequality  (4.38). 

We  arrive  at  the  inequality 

(4-40)  Zihfe)  -fff’-t  ql*f  (q2*  l3)ldltd5' d* 


t t t 

_2  _ „ _Q  /_£! 


* |/jtql+  qp*  q^dxdyd: 


31  - 


Rearranging  as  in  tho  preceding  section,  we  find 
(4,4l)  4TTK(t)  ' )~ 


t - t 



JjJ[q2*  q3'  ql'|dxdy  dz 


(t-  t ) 


rcc  2 

^cbcdydz 

o J J 


We  deduce  from  (4.4l) 
^*42^  dt  "4TTK(t) 


X 

t 


[q^  * - q2"ldxdy  dz 


D(t) 


which  agrees  with  (4.29).  Furthermore,  ne  derive  from  the  same  inequality 

.2 


(4.43) 


r 1__-'  A. 

dt2  4rr*^  t* 


. <hdxd>,dz-  tii  uwfe1 


Combining  (4.42)  and  (4.43)  we  find 


'4-44^  [4TTK(t) 


i"  ^ 


[ q~  * q2  ’’dx  0 


}y  dz 


D(t) 


and 

(4.45) 


f 1 V ♦ 1 r 1 , I ' y -i. 

4TTK(t)  t L4hTK(t)  • ~ 2 


tqj  ♦ <\2  + q^xdy  dz 


D(t) 


This  inequality  may  be  interpreted  as  follows: 

Theorem  4b:  Under  partial  one-sided  stretching  the  functional  tK(t)  ^ 

i3  convex  from  above  as  a function  of  the  variable  s ■■  t*~ . 

This  result  holds,  of  course,  a fortiori  in  the  case  of  one-3ided 
stretching  of  the  whole  space;  it  is,  however,  a weaker  result  than  the 
convexity  of  tK(t)  1 as  a function  of  t which  wa3  proved  in  Section  4.2. 
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5 • Angular  Stretching . 

5,._I.  It  is  possible  to  study  the  change  of  the  preceding  and  similar 
functionals  for  various  other  types  of  deformations  in  the  sane  way  as 
we  did  in  the  case  of  one-sid'd  stretching.  The  extremum  property  of  the 
functional  combined  with  the  transplantation  of  the  extremal  function  leads 
egain  to  statements  on  monotonic  or  convex  dependence  on  the  parameter  of 
the  deformation. 

We  mention,  as  an  example,  the  deformation  by  angular  stretching. 

Let  D be  a domain  in  the  (x,y)-plane  and  let  0 be  a point  outside  D. 

We  choose  0 as  the  origin  of  a polar  coordinate  system  r,  and  assume 
that  D appears  from  0 under  an  angle  8-c  2'TT.  We  subject  the  whole  plane 
to  the  deformation 

(5.1)  r'  - r , <p'  - tCp  , 0 ^ t < ^ 

and  define  D(t)  as  trie  image  of  D -under  this  transformation.  Clearly, 

D(t)  lies  still  3imply  over  the  plane  and  the  mapping  D— >D(t)  is  one-to-one. 
We  shall  say  that  D(t)  has  been  obtained  from  D»D(l)  by  angular  stretching 
with  the  parameter  t. 

The  torsional  rigidity  P(t)  of  the  simply-connected  domain  D(t)  can 
be  estimated  by  means  of  the  inequality  (1.5)  and  by  using  the  transplanted 
stress  function  from  the  domain  D(tQ).  It  is  useful  to  bring  the  formula 
(1.5)  into  polar  coordinates  and  to  consider  the  functions  F"?(r,<p}. 

We  have 


(2 


(5.2) 


P - max 


jj  F(r,  <f)r 


dr  d cfY 


jf  ^2  ^rdrdcp 

D r 

Using  the  same  method  as  in  Section  1.1,  we  find  the  inequality 


4 
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Z t 


l5-3>  PU) 


_ 2 

where  Z - t and  f(r,^)  is  the  correct  stress  function  for  the  domain  D(tQ), 
Wo  3ee  from  (5.3)  that  the  curve  of  tP(t)  ^ plotted  versus  the  abscissa 
Z ■ t has  at  every  point  a supporting  line  with  the  slope 


0 J u -<■  % 

D ( t ) 

a 


-^^rdrdcp*  jj  f^rdr 


dcf 


D(t„) 


(2  l 1 f (r , O ) r dr  d<E )‘ 

D(t'  ) 


(5.4) 


- t3P(t)‘2  jj  rrdrdCf> 
D(t)  r 


0 


W may  thus  state 

Theorem  5.  Let  F(t)  be  the  torsional  rigidity  of  the  3 imply- connected 
boundod  domain  D(t)  obtained  from  D-D(l)  by  angular  stretching  with  the 

parameter  t.  Then  tP(t)~3  is  increasing  and  convex  from  above  as  a function 

^ -2 

of  th6  variable  u - t 


5.2.  We  obtain  from  (5.4^  and  (1.4 ) the  estimate 

(5.5)  ^ [tP(t)_1l^l -^[tP(t)_1l 
This  inequality  can  be  written  as 

(5.6)  log  :t3P(t)'1l  0 

which  shows  that  tJP(t)  decreases  with  increasing  'Z  or  that 

(5-7)  P(t)t  3 decreases  with  increasing  t 

while 

(5.8)  P(t)t  1 increases  with  t. 

These  two  estimates  give  upper  and  lower  bounds  for  P(t)  if  P(l)  is  known. 

We  can  derive  from  (5.4)  also  the  derivative  of  P(t)  with  respect  to  t 
and  obtain 


a formula  very  similar  to  (1.13). 
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Chapter  II 

On  a Theorem  of  Poincare 


I • iiaifiseai  _&a4_Pr<2<2l . 

l.I.  We  consider  real  valued  functions  u,  v,  u -,u  , ...  <p.  %■■■■ 

that  belong  to  a certain  class  21.  We  assume  that  if  u and  v belong  to 
2*.  and  c is  a real  constant,  also  u*v  and  cu  belong  to  21.  We  consider 
§ symmetric  bilinear  functional  A[u,v!  over  21.  That  is,  Alu.v"1  is  defined 
for  any  functions  u and  v of  the  class  ZL  and  has  the  properties 

(1.1)  Atu^  + u^vl  - ACu^vl  ♦ A[u2>vl  , 

(1.2)  A[cu,v'1  - cA[u,v]  , 

(1.3)  A[u,v]  - A[v,b] 


Let  B[u,vl  denote  another  symmetric  bilinear  functional  over  XT* 
(Therefore,  (l.l),  (1.2),  and  (1.3)  remain  valid  if  we  substitute  B for  A.) 
We  assume  that  B is  definite  positive:  that  is, 

(1.4)  B[u,u1>0 

unless  u vanishes  identically.  With  the  functionals  A[u,vl  and  Bru,vl  we 
form  the  Rayleigh  rails 


(1.5)  Rtu1  - 

Rfu1  is  defined  for  a 
vanishing  function. 

(1.6)  RCcu1 


aLu.uI 

B [ u , u 1 

11  functions  u belonging  to  21,  except  for  the  identically 
If  c is  any  constant,  c/0, 


- Riul 


according  to  ( 1 . 2 ) , (1.3). 

We  assume  that  we  deal  with  one  of  the  many  examples,  familiar  in 
mathematical  physics,  in  which  it  is  possible  to  define  the  eigenvalues 
"Vy  V2,  . . . and  the  eigenfunctions  Cj>^,  cp^, . . . belonging  to  R[ul 


as  follows: 
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The  firot,  eigenvalue  is  the  minimum  of  Riu)  and  is  a function 
for  which  this  minimum  is  attained,  so  that 

(1.7)  » R[  qp  x 3 ^ R[  u"1 

for  any  function  u of  ^ (u  =0  excepted,  of  course).  In  view  of  (1.6)  we 
can  choose  cp^  so  that 

(1.8)  B(Cp1,CP11  = 1 . 


Now  assume  that  n -^2  and  ^ , . . . , ^ are  already  defined  so  that 


(1.9) 

for  i ,k  - 1,2 
conditions 

(1.10) 


. . ,n-l.  Let  us  seek  the  minimum  of  R[u"!  under  the  side- 

B^Cp^u1  - - ...  - - 0 


Let  be  the  value  of  this  minimum  and  & one  of  the  functions  for  which 
n ] n 

is  attained,  so  that  the  relations  (1.10)  imply 

(1.11)  -y  - R[Q  ^ R[ u l 

n in 

We  may  choose  <£>  30  that 
'n 

(1.12)  B(cpn,cpnl.l 

Then,  obviously,  the  relations  (1.9)  are  satisfied  for  i,k » 1,2, . . . ,n-l,n. 


1*2.1  There  is  a simple  and  efficient  method  to  estimate  the  eigenvalues 
■^I’-V-,...  that  3eems  to  have  been  first  observed  by  Poincare  [11,  p.  259-  261 
Let  g^,g2,...,gn  be  any  m lirearly  independent  functions  of  the  class 
C.  and  let  the  real  variables  x^,x2,...,xn  range  through  all  systems  of 
values  except  0,0,..., 0.  Then 


“Vi'Vz*  -‘V.1 


zrz:x1xv*(g1,gki 


(1.13) 
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by  virtue  of  (l.l),  (1.2),  (i.3),  and  (1.5).  By  virtue  of  (1.4),  the 
denominator  is  a definite  positive  quadratic  form  in  x^ ,x^ , . . . jX^. 

We  assume  that  th?  quadratic  foxas  and  Zl2Tbikx1X)(> 

ihs  latter  el  which  is  auppaaed  tQ  fee  definite  positive,  have  the  property 
that  lor  all  values  cl  x^,x2, . . . ,xm  (o,o,..., o,  qX  course,  e^cspt&d) 

2Z2Ia,ux  x 

( 1 • 14)  R(x  g * x g * ...  * x r 1 — 

il  c c mm  «r» r». 


saik’ti’tk 

is  let  V',  -i'~ V ' denote  the  roots  el  the 

1 c m 

(1.15) 


:isti£  equation 


a...  - IS b ...  a.  - 12b 

11  11  lm  lm 


a . - -^b  ....  a - TSh 
ml  ml  mm  mm 


0 


yi  2 3 — Vm 


increasingly  ordered  55  tfeat 

(1.16) 

Then 

(1.17)  i>  <=  -V'  * ^ V ■*/' 

1 1 <.  mm 

We  follow  the  usual  notation  for  quadratic  forms  in  assuming  that 

a..  * a,  . , b - b . 
ik  ki  ik  ki 

In  order  to  prove  the  theorem  stated,  we  make  U3e  of  a well-known 

theorem  of  Algebra  [ 1 , Chapter  13,  Sect.  59,  Theorem  III.  ITe  introduce 

m linearly  independent  linear  homogeneous  functions  y-,yof...,y  of 

1 c m 

X— , x. , . . . , x so  that 
1 c m 

(1  18)  . v(yi ‘ 

KL.LoJ  _ 2 2 2 

^ikVk  W -*y» 

We  set 
(1.19) 


u - xlgl*x2g2*  ...  + xmgm 


1 


:vn-  Lder  un  integer  r. . 1 — r. 
v’.iO'  and,  in  addition,  to 
(l  20'  v - v 

' ' Vl  Jn*2 

Thus  we  set  up  a system  of 


- m,  subject  x^x., . 
the  m-n  conditions 
. . . " y_  " 0 


,x  to  the  n-1  conditions 
m 


(n-l)  * (m-n)  - m-1 

homogeneous  linear  equations  for  the  m unknowns  x^  ,x^ , . . . ,xm*  Thi3  system 
possesses  a solution  different  from  0,0,...  ,0  which  we  call  x^ ,x^ , . . . ,x  . 

As  these  numbers  do  not  all  vanish  and  g.  ,g„  , . . . ,g  are  linearly  independent, 
the  function  u,  Riven  by  (1.19),  does  not  vanish  identically  and,  since  it 
satisfies  (1.10),  satisfies  also  (l.il).  We  infer  from  (l.ll),  (l.L4), 


(1.18),  (1.20),  and  (l  16)  that 

(1.21)  V 4 R[u ] 

n 


Hi 


V. 


2*2* 


V'/ 

an 


7?  + 


2 

^2 


• + y. 


and  so  we  have  attained  the  aesired  conclusion  "V  for  n ■ 1 ,2  , . . . ,ra. 

2 . Remarks  and  Applications . 

2.1.  It  is  important  to  observe  that,  in  view  of  ( 1 .13 ) * we  satisfy 
the  condition  (1.14.)  if  we  choose 
(2.1)  aik  - A[gi>gk1  » bik  * BUi»gk) 

We  may  make  the  following  deduction  from  the  fact  that  provides 
the  minimum  of  Rfc^1,  under  the  side  conditions  (1.10).  Let  4^  be  an  arbitrary 
function  of  the  class  which  is  orthogonal  to  ^ , . . . , Cp  ^ • Then,  we  have 
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by  (1  l).  (1  .2),  and  (l  3),  -he  equation 


Ai  CP  , Q 1 * 2tA[  q , 'f  ] + t A[  V , 'f  1 
(2.2)  R [CP  nfi  ■ -a  n ^ ~ 

R B['S  , CP  1*  2t3[  CD  ,<P1*  CBtHVV! 

n n in 


The  condition  that  this  ratio  be  a minimum  for  t-0  leads  by  elementary 
calculus  to  the  requirement 
(2.3)  A[epn,4> 1 - ^n9[Cfn’4M 

In  particular,  we  have  for  m r. 


3/  if  n ■ m 

‘ Vt<V?.1-”  ,,  , 

0 if  n r m 

Since  this  equation  is  symmetric  in  n and  m it  must  hold  for  arbitrary  n and  m. 

We  observe  now  that  we  can  obtain  equality  simultaneously  in  all  m 

inequalities  (1.17).  For  (1-13/  yields,  by  virtue  of  12.3'), 

2 . _ 2 


(2.3")  R[x1Cp1 


*...+  V x 

♦ - x (■  r 1 • — *— • Ai ■ — ■*»  -» 

••  xn  in J ^ o o 


2 2 

y ♦ r ♦ 

X1  *2 


Thus,  if  we  choose  gi  ■ cp4  and  (2.i),  the  determinantal  equation  ( 1 . 15 ) 
turns  out  to  be  simply 


(V  - V ) ( V - -j/  ) . . (V  - y ) - C 
1 ra 

and  so  becomes  V and  all  the  m inequalities  (1.16)  simultaneously 
become  equations  . 


2.2.  We  consider  as  given  the  m linearly  independent  functions  §^>§2* • • • >ga* 
Their  linear  combinations  (1.19)  form  an  m-dimensional  linear  snbspace  S 

CD 

of  Z.  Let  us  vary  x.,x  , ...,x  , that  is,  let  u vary  in  this  subspace  S , 
and  let  us  seek  the  maximum  of  R[ul.  If  we  assume  (2 . 1 ) , we  easily  conclude 
from  (1.13),  (2.1),  (1.18),  and  (1.16)  that 

(2  .4)  max  R[u^  - u 

, o m 
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Take  now  the  case  considered  m the  foregoing  Section  2.1  when  for  a 
suitable  choice  of  the  subspace  S (namely  g.  * , g- ■ g * Cf  ) 

the  inequalities  (l  17)  become  equations.  Considering  the  inequality  (l.2l) 
and  (2.4),  we  see  that 

(2.5)  V “ Min  max  Rlu] 

” S uLSm 
a ra 

We  used  capital  M in  Min  and  lower  case  m in  max  to  emphasize  the  difference 
between  a variational  extremum  problem  and  an  elementary  extremum  problem 
that  is  concerned  only  with  a finite  set  of  variables.  The  minimax 
definition  (2.5)  of  the  eigenvalue  ~U  is  very  different  from  the 
well  known  definition  due  to  Courant  [2,  p.  113  \ which  involves  two  successive 
variational  extremum  problems.  In  the  case  of  finite  quadratic  forms  (when 
is  a finite  dimensional  vector  space),  (2.5)  is  due  to  E.  Fischer  (4  7- 
We  may  take  the  liberty  to  call  an  m-diraens ional  linear  sutspace 
of  7Z.  an  m-dimen3ional  "convoy"  . As  the  ships  in  a convoy  must  travel 
together,  the  speed  of  a convoy  depends  on  the  slowest  ship.  Even  if  we 
wish  to  make  R[u"'  as  small  as  possible,  (2.4/  represents  the  best  that  S 

m 

can  do  for  R[ul  provided  that  is  considered  as  a convoy,  as  a set  with 
solidary  elements.  If  we  seek  the  best  m-diraens ional  convoy,  (2.5) 
represents  the  solution  of  our  problem.  And  so  we  may  describe  -v  , the 
eigenvalue  of  R[u\  as  the  "m-dimensional  convoy- minimum  of  Rfu’-". 


2.3.  The  present  section  deals  with  an  elementary  remark  on  the 
relation  (l.l?)  that  we  shall  need  in  the  next  section. 

We  consider  finite  sequences,  ®ch  consisting  of  m real  numbers,  such 
as  c*  • • • , ck^  or  (5^,  • • • » (3m*  Say  that  the  sequence 

C*i  » £*i  , • • • , c*i  is  a rearrangement  of  0^ , c*2  > • • • jCk  m if  i^,^,  •••  »im 
are  different  positive  integers  none  of  which  exceeds  m.  Let  , ck^ , . • • , c* 
denote  the  rearrangement  of  . . *,0^  in  ascending  order  so  that 
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(4  6) 


A ^ (/ 

x -t  m 


We  say  that  i\,o<2.  ..cX^  is  valorized  by  (2^, . . - , (S, 


if 


(2.7 


</,  4 %*  • -,£*  ^ j5_ 

1 ^ i m m 


We  express  the  fact  that  ex' ^ , . . . , ;x^  is  majorized  by  *n  writing 

( *-  • 8 ) ^ 1 ’ 2 * ’ ^ ^ ^ \ ' ^2  * ’ 1 ' * 

Thus,  as  the  and  ~V  ' are  arranged  in  ascending  order,  "24  3 "Z^, 

A “ 22  and  ?;e  may  express  (1  17)  in  the  more  concise  form 

(2.9)  V V • • v « 3A' . V,' V' 

l i m 1 m 

(a)  tes'xme  that  ihere  aiLS  rearrangements  cx^  » c*  and 

1 ^"2  m 

’ Pk  ’ * ’ ‘ ’ ^k  -■  — * fiX  ihi  3 , respectively,  gush  £bai 

(2.i°)  *V  - ?k  ’CX.^  (^k’ ‘ ' ‘,C*i  - ($k 

TMn  (2.?)  Iflllfiwfi,  and  M iiie  o<  are.  saj£xi££d  by  liia  (5. 


Since  a simultaneous  permutation  of  the  o<  and  the  % leaves  obviously 
intact  the  essential  hypothesis  of  this  lemma  (a),  we  can  assume  in  proving 
it  without  loss  of  generality  tnat  the  inequalities  (2.10)  are  reduced  to 
the  form 

(2.11)  ~ ^,*2  “ '<*m  — 


m 


where  ^ , . . . , ^ is  a rearrangement  of  the  ^ . If  no  pair  of  subscripts 
i and  k exists  such  that 
(2-12!  1<k-  ■ 

(°.ll)  coincides  already  with  the  desired  (2.7).  If,  however,  there  is 
such  a pair  i,k  as  described  by  (2.12),  we  infer  from  (2.1l),  in  view  of 
the  ordering  of  the  cX  (see  (2.6))  that 

<*i  " *k  - ^ Pi  ^ 
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and  so 

*1 ' bk- 

Therefore,  interchanging  and  we  bring  (2.1l)  one  step  nearer  to 

the  desired  final  situation  (2.7)  which  we  can  thus  attain  in  a finite 
number  of  similar  steps.  This  proves  (a). 


(b)  Suppose  thai  ihe  sequence  (A  >•••><*„  ia  aaJsxized  by  iM  aaqusi 

1 m 


and  there  are  bwq  elements . o<i  In  the  flxat  sequence  and  ,6k  in 
the  second , such  that  cx ^ " (3^  • If  is  delete  o< i and  |?k , thfl  remaining  o( 
(n-1  in  number)  are  still  aajcriied  by  th§  remaining  j 6 (also  n-1  in  number ) ■ 
If  one  of  the  defining  inequalities  (2.7)  happens  to  be  4=  ^3^, 
the  assertion  is  immediately  obvious.  Otherwise,  there  are  subscripts  i' 


and  k' , i'  / k, . k'  / i such  that 

(2.13)  ’ °V  ~ ?k 

occur  in  (2.7).  It  follows,  however,  that 


U.u)  'V^'3k  - Pi' 

Now,  by  virtue  of  lemma  (a),  we  may  rearrange  the  , and,  using  (2. 14) 
we  may  replace  (2.13)  by 

°V  - I3!'  ’ *i  - <?k 

This  brings  U3  back  to  the  immediately  obvious  case  and  proves  (b).  There 
is  no  need  to  formulate  the  obvious  extension  of  (b)  to  the  deleting  of 
several  pairs,  each  of  which  consists  of  an  and  a equal  to  that 
particular  (X.  We  shall  need,  however,  this  extension  in  the  next  section. 


2,4.  Let  us  assume  that  those  eigenfunctions  of  R[ul  that  belong  to 
eigenvalues  equal  to  the  same  number  V form  a vector  space  of  finite 
dimension.  (This  is  so  in  the  most  usual  problems  of  Mathematical  Physics.) 


If  this  dimension  is  (! , we  say  sometimes  that  ~P  is  an  eigenvalue  of 
multiplicity  £ , or  (better)  that  eigenvalues  coincide  w ith  . (In 
fact,  v is  the  common  solution  of  i different  extremal  problems,  see 
section  1.1  or  (2.5).) 

We  introduce,  however,  a further,  more  special,  assumption.  Let  p 
denote  the  set  of  variables  (the  point)  on  which  the  functions  of  the  class 
y*.  depend  We  assume  that  there  is  an  operation  J of  period  2 (J^  is  the 
identity,  J is  an  "involution")  that  transforms  the  point  p into  p,  so  that 
Jp  =*  p , Jp  -*  J(Jp)  * p 

We  assume  that  the  involution  J transforms  every  function  u(p)  of  the  class 
into  a function  u(p)  of  the  same  class,  and  that 

(2.15)  A[u(p)  ,v(p)  1 - Afu(p)  .v(p)  I , 3[u(p)  ,v(p) 1 =»  B[u(p),v(p)l 

(Example . Let  p stand  for  the  point  (x,y,z),  p for  (-x,-y,-z), 

D for  a bounded  domain  symmetric  with  respect  to  the  origin,  and  let  the 
class  C consist  of  the  functions  having  two  continuous  derivatives  within 
D,  continuous  within,  and  on  the  boundary,  of  D,  and  vanishing  on  the 
boundary.  We  set 

A[u,v1-  I(uv*uv*uv)dxdydz  , 
xx  yyzz 

B(u,v 1 - J ' | uvdx  dydz  ; 

both  integrals  are  extended  over  D.  The  conditions  (2.15)  are  satisfied. 

We  could  define  p differently,  as  (-x,y,z),  or  as  (-x,-y,z),  provided 
that  D has  the  corresponding  symmetry,  with  respect  to  the  y ,2-plane,  or 
the  z-axis,  respectively.) 

We  call  the  function  u(p)  even  or  odd.  according  as  it  satisfies  the 
first  or  the  second  equation 

(2.16)  u(p)  - u(p)  , u(p)  - -u(p) 

We  develop  a few  consequences  of  the  above  assumptions  and  definitions . 


U'i 


(fi)  II  u oqq  and  v la  even 
A [ u , v =•  B'u.v  ' = 0 

In  fact,  in  using  the  definition  of  even  and  odd  functions  and  (1.2), 
we  find 

A( u(p) ,v(£)  - Ar-u(p  ) ,v(p)  1 - -Afu(p) , v(p)  1 

In  view  of  the  first  equation  (-  15)  tne  assertion  concerning  A follows, 
and  that  concerning  3 follows  similarly. 

(b)  if  q>(p)  is  an  else&f unction  sX  Rtu1  ibai  belongs  in  ills 
eigenvalue  v,  aloe  C?(p)  is  an  eAgenfuaeiiSD  belonging  is  gar-g 

eigenvalue  v. 

In  fact,  Cp(p).  as  nr.  eigenfunction,  must  be  obtained  at  a certain 
step  of  the  recursive  process  described  in  Section  1.1.  If  this  step  is 
the  nth,  <qp(  p ) - qpn  ( p ) , ^ Yet,  in  view  of  (2.15).  we  can  sub- 

stitute in  the  defining  process 

^Ap) , , -,cpn(p) 

for 

<f,(P),<?2(p),..  ,qpn(P) 

without  altering  the  cor  responding  eigenvalues 

v ,-v_ , • • ,-v 

12  n 

This  proves  (b) . 

(c)  if  is  an  eigenvalue  oX  multiplicity  X,  the  eigenfunctions 

belonging  ifi  have  a bag  is  that  consists  OX  j even  and  k qM  functions, 

j*k-X,  j >0,  k >0 

Let  p)  be  an  eigenfunction  that  belongs  to  V.  By  (b),  Cp(p)  is 
also  an  eigenfunction  that  belongs  to  V,  Therefore,  Cp(p)  + Cf(p)  is 
either  identically  0,  or,  except  for  a constant  factor,  an  even  eigenfunction 
belonging  to  v„  Similarly,  qp(p)-  Cp( p)  is  either  identically  0,  or  an 
odd  eigenfunction  belonging  to  "ZA  Yet 
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U 1?)  q>(p)  - ^:q5(c)  ♦ Q ( p ) 1 4 o> (p)  * ^ (p) 1 

and  both  terns  cannot  vanish  cn  the  right-hand  3ide.  Thus,  there  exist 
either  even  eigenfunctions , or  odd  eigenfunctions,  or  both  kinds  exist. 

The  even  eigenfunctions  for m a linear  aubspace  and  possess  a basis 
consisting  of  j functions,  and  the  o’d  eigenfunctions  possess  a basis 
consisting  of  k functions.  (If  there  are  no  eigenfunctions  of  one  or  the 
other  kind,  j or  k may  be  0,  respectively.)  At  any  rate,  as  (2.17)  shows, 
any  eigenfunction  belonging  tc  U is  a linear  combination  of  these  j+k 
basis  functions.  These  j4k  functions  are  obviously  linearly  independent, 
since  only  the  identically  vanishing  function  can  be  both  odd  and  even, 
and  so  j4k-  H . This  proves  (c). 

We  can  describe  the  situation  viewed  in  (c)  and  its  proof,  by  saying 
that  there  are  j even  eigenvalues  and  k odd  eigenvalues  coinciding  with 
we  call,  of  course,  an  eigenvalue  odd  or  even  according  as  the  eigenfunction 
associated  with  it  is  odd  or  even.  Under  the  assumptions  of  the  present 
section,  namely  (2.15),  we  £_an  label  each  single  eigenvalue  as  even  or  odd. 

(d)  Assume  that  the  £unciians_  g^ , g2 , . . . , g^  cons idered  in  section  1.2 
are  even  and  that  aik  and  ai_fi  Elven  fcy  (2.]).  T^en  ihs  roots 

, . . . , v'  q£  ihs.  characteristic  squaiian  (1.15)  oajsxiie  not  only 
the  first  m eigenvalues , but  ihe  first  m EVEN  eigenvalues . 

It  is  also  true:  "If  th^  functions  g^.g^, . . . are  chosen  as  odd, 

i^’ , "h/ , . . . , "i/  majorize  the  first  odd  eigenvalues"  and  the  method  of  proof 
that  will  be  used  applies  equally  to  both  cases. 

Let  h denote  the  number  of  those  odd  eigenvalues  that  are  not  greater 
th 

than  the  m even  eigenvalue,  and  let  gm+i’£n4-2’  * • * »6m4.n  denote  the  odd 

eigenfunctions  belonging  to  these  n odd  eigenvalues.  Apply  the  theorem  of 

Section  1.2  to  the  m4n  functions  g, ,...,g  ,g  g . (instead  of  the 

1 m m4!  m4n 
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m functions  ) Applying  (a)  to  (.2.1) , we  see  that  there  are  two 

m by  n rectangles  in  the  characteristic  determinant  (see  ( 1 . 1 5 ) ) one  in  the 
north-east  corner,  the  other  in  the  south-west  corner,  all  elements  in 
which  vanish.  Thus,  the  determinant  is  a product  of  two  determinants,  one 
of  order  m,  the  other  of  order  n.  The  latter  has,  in  view  of  (1.9)  and  (2.2), 
non-vanishing  elements  only  in  its  main  diagonal  and  yields  as  roots  precisely 
the  n odd  eigenvalues  associated  with  the  theorem  of 

Section  1.2,  the  m+n  roots  of  the  determinant  majorize  the  m*n  first 
eigenvalues  of  R[u"U  The  n odd  eigenvalues  occur  both  in  the  majorizing, 
and  in  the  majorized,  set.  If  we  delete  these  n odd  eigenvalues,  :>nly  the 
m roots  corresponding  to  g^,g^, — ,gm  remain  in  the  majorizing  set  and 
only  the  m even  eigenvalues  remain  in  the  majorized  set.  Using  (b)  of 
Section  2.3,  we  attain  the  assertion  (d). 

It  does  not  diminish  the  interest  of  the  foregoing  proof  that  some 
more  or  less  definite  statement  approaching  (d)  has  often  been  tacitly 
assumed  in  the  computation  of  eigenvalues. 


- 46  - 


Chapter  III 
Eigenvalue  Problems 


1 . The  Convexity  Theorem . 

U..  In  this  chapter,  we  shall  combine  Poincare's  theorem  with  the 
method  of  transplanting  extremal  functions.  We  will  obtain  various  results 
on  the  dependence  of  eigenvalues  on  their  domain  of  definition.  We  formulate 
an  eigenvalue  {roblerr.  which  is  general  enough  to  cover  various  importer t 
particular  cases,  but  we  do  net  intend  to  complicate  our  presentation  by 
aiming  at  tne  greatest  possible  generality. 

Let  D be  a fixjd  finite  domain  in  the  (x,y)- plane  which  is  bounded  by 
a smooth  curve  C;  let  s be  the  arc  length  on  C,  counted  from  some  given 
point  on  C.  We  introduce  two  continuous  non-negative  functions  on  C 
(1.1)  h(s;  t:)  £ 0 , k(s;  ~ ) ^0 

which  depend  on  an  additional  real  parameter  "C  whose  role  will  become 
clear  in  the  sequel.  Let  further 


(1.2)  p ( x , y ; c ) > 0 , qvx,y;  H ) > 0 , r(x,y;  IT)  > 0 

be  three  continuous  positive  functions  in  D which  also  depend  on  the  same 
parameter  t . Finally,  let  21  the  class  of  all  functions  <^>(x,y)  which 
are  continuously  differentiable  in  D and  are  continuous  in  D*C.  We  define 
in  21  th®  functional 


(1.3)  M*!1  w 


£ \jhdx2  * k dy2^  [p<^2  * q^^dxdy 


II r 


dx  dy 


which  still  depends  on  the  parameter  't  . In  the  notation  of  Section  (II,  l.l), 
we  have 

V 


(1.4)  A[u,v]  - 

and 


2 2^ 
uv  A h dx  ♦ k dy  * 


[puxvx  + q ^y1^^ 


4”  - 


(1.5) 


B r u , w 1 


r uv  dx  dy 


We  define  a set  of  eigenfunctions  UyCx^y.'C)  end  eigenvalues  Xy(^  ) 

for  the  functional  R.-I $ 1 by  the  sequence  of  minimum  problems  described 

in  Section  (II,  l.l).  The  existence  of  extremal  functions  is  guaranteed 

by  the  general  calculus  of  variations  and  it  is  al3o  well  known  that  the 
th 

V eigenfunction  u satisfies  the  differential  equation 


(1.6) 


_2L  , 

dr  (p  ?r  } 9y  (c* 


ay  ) ♦ Av  ruy  - 0 


with  the  natural  boundary  condition  on  C 


. , av  . , I 

>vyh 


2 ^ 

(1.7)  p ~~  cos  (n,x)  * q cos  (n.  y)  - uvVh  cos  (n,y.>  *■  k cos'*  (n,x)  , 


where  n is  the  interior  normal. 

In  view  of  the  conditions  (II,  1.9)  imposed  on  the  extremum  function, 
we  have  the  orthonornal izotion  conditions 
(’••8)  B-VyUyl  - 

where  we  denote  as  usual  by  the  Kronecker  symbol 


(1.9) 


0 if 


tv 


yv 


1 if  I J • V 

I 

On  the  other  hand,  we  derive  from  the  general  statement  (II,  2.3/)  and  (l.S) 
the  set  of  equations 
(1.10)  Aru^,uy]  - 

We  introduced  the  eigenfunctions  u^  and  their  eigenvalues  by  a 

sequence  of  extremum  problems  and  showed  that  they  satisfied  the  linear 

partial  differential  equation  (1.6)  with  the  linear  homogeneous  boundary 

condition  (1.7).  It  can  be  shown,  conversely,  that  the  system  (1.6),  (l„7) 

has  only  a countable  set  of  eigenvalues  Av  which  tend  to  infinity  with  "V 

and  that  this  set  coincides  with  the  set  of  extremum  values  of  the 

functional  R,[<£1.  In  most  applications,  one  is  interested  in  solutions  of 
Z * 
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the  differential  oyster?,  and  the  extremum  problems  considered  are  introduced 
for  methodological  reasons. 


1.2.  We  wish  now  to  make  assumptions  regarding  the  dependence  of 
the  functions  (l.l)  and  (i.2)  on  the  parameter  . We  put 

(1 .11)  xv  ’ T) 

and  are  going  to  draw  conclusions  on  the  dependence  of  the  eigenvalues  upon 
the  parameter.  Let  us  consider  at  first: 


iiypQttieaia  I.  h,  k,  p,  and  q ar_s  continuously  differentiable  and 
convex  from  abqye  in  dSFandence  ^ f . r is  Independent  qX  'C. 

We  want  to  prove  the 

N 

Conclusion-  Ibe  5_ua  Av(t:  ) ii  fqr  arbitrary  iota  gar  N convex 

V-  l 

fcofl  abovfl  os  a function  qX  't  ■ 

In  fact,  let  U3  write  for  short  h-h(s;'C),  hQ-h(s;t:o),  etc.  where 
is  an  arbitrary  but  fixed  value  of  the  parameter.  By  our  hypothesis, 
we  have  at  each  point  of  C: 


(1.12) 


h^h  MC-T)h.  , fu  - tV  h(s;  D 

o O i 1 C 7 >- 


r-  r 


and  three  analogous  inequalities  for  k.  p,  and  q.  Observe  that  by  the 
inequality  between  geometric  and  arithmetic  means,  we  have 

(1.13) 


/ o o o Z~l  h*h  k ♦ k _ 

( hQdx  * kQdy  ) (h  dx  ♦ k dy  ) 4z.  ~ ^ — dx  * ~a~z — dy 


Z-Z 


T-Z 


(h  + — r— 2 h.  )dx2  ♦ (k  + — ~ 2 k1  )dy2 
o 2 1 o ^ 1 J 

We  define  now  the  functional  R^.  [ which  is  obtained  from 

’ o 

if  we  replace 

(1.14)  p by  p0+(x_c0)pi  * q by  qQ+^'1'o^qi 

and 

J ; p 2 2 -1/2  T-t  ? C-r 

(1.15)  /|  hdx  + k dy  by  (hQdx  * kQdy  ) [ (hQ  + )dx  + (kQ  * — 5 k^)dy~  1 
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In  view  of  the  inequality  (l  13)  and  the  inequalities  for  p and  q which  are 
analogous  to  (1.1.2).  we  can  assert  that  for  each  c|>  of  the  cla33  we  have 

(1.16)  !U<£1 &R'  „ f $1  . 

’“o 

We  may  put 

V,  ■_ 

(l.l?)  » ' r iT. i _ — 1— Q.. 


r/  [ * 1 - — 2-a- 

1,f'o 


Let  ^ »<$2 » • • ’ the  N first  eigenfunctions  belonging  to  R^.  [<$1. 

o 

We  form 

(1.18)  u » x1£1*  * • * xn$N 

and  have  by  (i„16),  (1.8),  and  ( 1 . 10 ) the  Inequality 


(1-19) 

with 

(1,20) 


f £*lkv* 

tu1  ^ ^ ^ 

iv 

i 1 

hk  * ‘r.r1  ^-^k' 


We  denote  by  .,Aj,  the  eigenvalues  of  the  symmetric  matrix 

arranged  in  increasing  order  of  magnitude.  According  to  Poincare's  theorem, 

we  have 

(1.21)  <=  h'j,,  V-1,2,...,N  . 

since  the  sum  of  the  A,,  equals  the  trace  o'  the  matrix  ((w^))*  we 


derive  from  the  definition  of  R^.  ^ 


(1,22)  YZ  - 

32-1 


x-t  ? f-i: 

(Rq-t0  hi)dx  l (ko'~rfl  ki)dy  2 

2 . , , 21  -■  - ^ 


h dx  + k dy 
o o 


12-  1 


* jj{‘V  <r-'to)c>l1  Z^(^)2]dldy  • 
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Thin  inequality  has  the  form 

. N . N 

(:  .23) 


*' - 1 V ■ 1 

The  geometric  interpretation  of  this  result  is  obvious:  If  we  plot 

N 

the  function  y -»  ^ Ky(V  ) against  the  variable  TT  , the  curve  lies  always 
below  the  straight  line  y ■*  yQ  ♦ ("t  - XQ)ajj  , passing  through  the  curve 
point  ( TT Q , yQ ) . Thus,  the  curve  has  at  every  point  a supporting  line 
and  is  convex  from  above.  This  proves  our  assertion. 

We  find  from  (1.22)  the  slope  of  the  supporting  line 


(1.24) 


1 r hidx2  4 kidy2 


sz:$ 


2 V-  1 


2 

| 

-V 


Aj hQdx2  * kQdy 

Let  us  suppose  that  for  some  given  value  of  'Z  the  N first  eigenvalues 
Kyiv)  are  all  different.  In  this  case,  it  can  be  shown  that  the  \y{V) 
have  derivatives  with  respect  to  T if  the  coefficient  functions  h,  k,  p, 
and  q ere  continuously  differentiable  functions  of  'C  . Clearly,  we  find 


(1.25) 


d 1 [ hidx2 * kidy2  *2  fr,  ,d&‘2 

dT  '2.1  | o ? 


J i , j 2 , . 2 

C ^fh^dx  ♦ k_dy 


th 


o o 

In  order  to  discuss  the  case  of  a multiple  eigenvalue  let  us  suppose 
that  Ajj  has  the  two  eigenfunctions  and  There  i3  an  ambiguity 

which  eigenfunction  should  be  taken  in  (I.24)  and  we  may  use  as  N 
eigenfunction  a3  well  the  linear  combination 

(1.26)  , *2+( S2  - 1 

This  leads  to 

(1.27)  aN  - aN_1*  <*2K  + 2*^1  * (32  M 


where  K,  L,  and  M are  easily  computed  from  0.24)  and  (1.26).  Except 
for  the  singular  case  L>0,  K-M  we  have  therefore  an  infinity  of  supporting 
lines;  the  minimum  and  maximum  slope  of  this  family  can  be  easily  determined 
from  (1.2*7).  Thus,  the  curve  y(TT)  will,  in  general,  have  a corner  point 
if  the  la3t  eigenvalue  ir.  the  sum  is  degenerate  and  possesses  an  eigen- 
function + T • 

An  extension  of  these  considerations  to  the  case  of  a greater  degeneracy 
of  the  eigenvalue  is  obvious.  We  will,  in  general,  find  a corner  for 
the  above  curve,  due  to  the  fact  that  the  freedom  of  choice  for  leads 
to  supporting  lines  with  different  slopes. 


1^4.  In  this  section, we  shall  make 

iiypQtteala  II-  h,  k,  p,  and  q aia  iadspaMent  of  'r  and  r is  contin- 
uously differentiable  and  convex  from  below  in  T . 

N 

Conclusion . The  sum  )>  ( 'll ) is  convex  from  below  m Z for 

arbitrary  choice  q£  ids  iniagex  N. 


The  proof  is  the  same  as  in  the  previous  case.  We  have  by  assumption 


(1.28) 


r > r * ( Z - T )r. 
o 1 


. 3 iLi. 

1 dZ 


if  r - r (x,y;  T- ) , ro*r(x,y;'  . We  replace  in  R^t^t  the  function  r in 

the  denominator  by  r ♦ ("C-  "T  )r.  and  obtain  so  a functional  R-.  _ (3)  1 

o oi  ZtV0  x 

for  which  holds  by  (l  28): 




(1.29) 


jj  CrQ  ♦ (T-  'tQ)r1  1 $ dxdy 


for  arbitrary  $ of  class  21 • 


Let  now  ^^x,y)  (■v' ■ 1,2,  , . . ,N)  be  the  N first  eigem  unctions  of 


(o) 


R [ ^ 1 and  denote  their  corresponding  eigenvalues  by  V - We  put 
To  V 
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(1.30) 


X.  X0 

CK  + — ;2 

fh[o) 


, $,*••• 


9„ 


By  (l.8),  (1.10),  and  (1.29),  we  have  the  Inequality 


(1-31) 

with 

(1.32) 


Hf[u]  ^ 


ik 


x2  ♦ x2  ♦ f x3 

!l.-  _ ? llllZH 

z£bik*i\ 


\ \ tr  * (1--  r )r.  1 
o o 1 


iiA. 


o)  v(o) 


=;  dx  dy 


We  denote  by  A^,...,/A^  the  eigenvalues  of  the  right-hand  side  quadratic 
ratio  in  ( 1 . 3 1 ) and  have  by  Poincare’s  theorems 


(1.33) 


^■1,2 N 


The  Xj,  are  defined  by  the  secular  equation  (II,  1,  15)  with  a^"  <5^ 


and  the  b^  defined  in  (1,32).  This  equation  may  be  brought  into  the  form 


(1.34) 


bix-  (A'r! 5ik!  ■ 0 


Thus  ,we  have 

N , N 

(1.35)  n -jr  - n \y- 

V-l  hv  v - 1 
and  by  (1.33) 


r - <fc2 

[r  ♦ (T.  - Z )r.  1 7ft  ( y dx dy  , 

Jo  • 1 M 


(1-36)  ft  Tfa  * ft  Vv* (c‘ ft  xfe dId7  ■ 


Thus,  we  have  shown  that  the  curve  y,. ( 1 ) - ^ * A^tT)  A ha 3 at  every  point 

v-1 

a supporting  line  which  lies  entirely  under  it  and  this  proves  our  statement. 
The  slope  of  the  supporting  line  is 


cl 


(1.37) 


ri.ft  V^T^7 


- ] 'V' 

U 

This  formula  permits  us  again  to  calculate  the  derivatives  X^{Z  ) if  they 
exist  and  to  discuss  the  corners  of  ..he  curve  y^(T  ) if  degeneracy  of  eigen- 


value occurs. 


It  should  be  observed  that  we  have  to  assume  In  this  section  that  the 
functional  A[<£,c|v  positive- definite  in  the  class  21.  In  fact,  we 
operated  with  the  reciprocals  of  the  eigenvalues  and  the  inequality  (1.36) 
becomes  meaningloss  if 


2 . Application  to  Stretching  q£  Domains . 


2,1.  In  the  preceding  sections,  we  considered  a variational  problem 
in  which  the  coefficients  dependec  upon  a parameter  'C . We  use  the  method 
of  trans plant Lng  the  extremal  function  insofar  as  we  estimated  the  minimum 
for  a given  value  'Z  by  using  the  correct  extremal  function  for  T-  as  a 
comparison  function.  In  the  following  sections,  we  want  to  consider  fixed 
variational  problem:;  with  respe-t  to  domains  which  depend  on  a parameter  , 
By  referring  the  variable  domains  back  to  one  fixed  domain  we  will  reduce 
the  problem  to  one  of  the  previous  form,  namely  a variational  problem  for  a 
fixed  domain  with  coefficients  depending  upon  a parameter  T . Thus,  the 
results  of  the  previous  sections  will  become  applicable  to  the  problem  of 
varying  domains. 

We  shall  study  the  minimum  problem  for  the  functional 


(2.1) 


The  corresponding  eigenfunctions  uv(x,y)  will  satisfy  the  Euler- La  grange 
equations 

(2  .2)  Au^*  uv  » 0 


with  the  natural  boundary  conditions  on  C 


Clearly,  the  eigenvalues  X v depend  on  the  parameter  k From  Poincare's 
theorem,  we  derive  the  result: 

The  eigenvalues  \,,(k)  increase  monotonically  with  k.  From  the  theorem 

of  Section  1.2,  we  find  further: 

N 

Theorem  1.  The  sum  ^ \y(k)  is  for  arbitrary  integer  N convex  from 

2V-i 

above  as  a function  of  k . 

We  may  apply  (l  24)  in  order  to  find  the  derivative  of  the  above  sum 
whenever  it  exists;  we  have 


(2.4) 


ds 


if  uy  is  the  Vth  eigenfunction  to  the  value  k of  the  parameter. 

If  k -~>od,  the  eigenvalues  X ^(k)  tend  to  limit  values  Xy  which 
correspond  to  the  eigenvalue  problem 
(2.5)  Au  * Avu  - 0 , u - 0 on  C , 

the  problem  of  the  vibrating  membrane  spanned  over  D and  clamped  in  at  C. 


2.2.  We  consider  the  family  D(t)  of  domains  which  are  obtained  from 
a given  domain  D«D(l)  by  the  one-sided  stretching  in  the  x-direction 
(2,6)  x ' - tx  y ' ■ y 

The  extremum  problem  for  the  functional  R[c£]  over  a domain  D(t)  may  be 
expressed  as  the  analogous  problem  for  the  functional 

k £ <|2  ?jdx?  * t 2dy2  + [ t 2 <£2  ♦ (J2  Idx  dy 

(£>2  dx  dy 

and  with  respect  to  the  original  domain  D-D(l). 

-2 

We  put  t ■ X and  observe  that  the  hypothesis  I of  the  Section  1„2 
is  fulfilled.  We  may  thus  draw  the  following  conclusion: 
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Theorem  2.  If  the  domains  D(t)  are  obtained  from  D-D(l)  b 7 ©ne- 
tt 

sided  stretching  with  the  parameter  t,  the  suras  ’ Ay(k)  are  convex  from 

-2  *“1 

above  as  iunctions  of  the  parameter  'U  m t for  arbitrary  fixed  values  N 
and  k. 


From  (1.24)  we  calculate 


(2.8) 


j.  ^ 

dr 


v« 


A v(k) 


dxtfy 


This  formula  shows  that  the  suras  increase  monotonically  with  the  variable 

We  consider,  in  particular,  a rectangle  with  sides  of  length  a,  b and 
parallel  to  the  x-  and  y-axis  respectively.  The  eigenvalue  problem  (2.2), 
(2.3)  with  k-0  ha 3 the  eigenfunctions 


(2.9) 

with  the  eigenvalues 


iTT 


slU 


u^  - cos  x cos  T“  y 

' a o 


(2.10) 


Av  ■ ^(4-  * *9>  . ^>0,  -so 


a b 

Under  a stretching  in  the  x-direction  with  parameter  t,  we  thus  obtain 


(2.11)  xi  v-  ^ ~T  n y*2 

V*  la  b 

where  the  summation  is  to  be  extended  over  all  pairs  ( J2 , m)  which  lead  to 
the  N first  eigenvalues.  As  long  as  no  \ is  degenerate  the  coefficients 
in  the  right  side  linear  function  remain  unchanged  and  the  left  side  sura 
is  a linear  function  of  T . This  result  shows  that  the  convexity  statement 
of  the  preceding  theorem  i?  the  best  possible  statement  with  respect  to  the 
second  derivative  of  the  sura,  considered.  If  Ajj  is  degenerate  for  X - 
and  we  have  the  two  representations 


K - n <Tc  2 

a 


♦ ) - TT2( 


L 


12  + 5 

2 .2 


(2.12) 


Kith  we  have  to  use  the  pair  (^,m^)  in  the  formula  (2.1l) 

as  long  as  t 6:  CQ  hut  have  to  replace  it  by  for  rC~?' 

Thus,  the  slope  of  the  line  11 ) suffers  a discontinuity  at  the  point 
of  degeneracy  T and  it  is  clear  that  its  value  decreases.  This  agrees 
with  the  convexity  theorem,  which  covers,  however,  much  more  complicated 
situations  than  the  one  discussed  here. 

Similarly,  the  eigenvalue  problem  (2.5)  of  the  membrane  leads  in  the 
case  of  the  rectangle  to  the  eigenvalues 

(2.13)  Ay-  TT2(-^—  ♦ °r)  , JL  > 0 , m>  0 

a b 

and  an  analogous  situation  arises. 

2.3.  The  value  of  the  convexity  theorem  of  the  preceding  section  comes 
from  the  fact  that  it  permits  to  estimate  the  eigenvalues  of  a whole  class 
of  domains  D(t)  if  the  corresponding  eigenvalues  are  known  for  a few  domains 
of  the  class  only.  We  want  to  illustrate  this  fact  by  a special  case  which 
will  be  iiscussed  in  detail. 

Let  D be  the  square  with  diagonals  of  length  IT;  we  assume  that  the 
coordinate  axes  coincide  with  the  diagonals  of  D and  we  consider  the  family 
of  rhonbi  D(t)  obtained  from  D by  one-sided  stretching  in  the  x-direction 
with  parameter  t.  The  eigenvalue  problem  to  be  studied  for  this  clahs  of 
domains  is 

(2.L4)  Au+  Au  - 0 , - 0 

a n 

The  lowest  eigenvalue  of  this  problem  is  always  A ■ 0 corresponding  to  a 
constant  eigenfunction.  We  want  to  study  the  value  of  - A2  Tor  this 
problem  and  for  the  rhorabi  D(t). 

By  virtue  of  the  convexity  theorem  of  Section  2.2  with  k-0,  N«2,  the 
function  u depends  convexly  upon  X - t . We  shall  compute  now  u(T  ) for 
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the  value  1.  » 1 , that  la  for  the  square,  end  determine  two  supporting  lines 
through  this  p^int  to  the  curve  |>(T).  This  will  lead  to  an  upper  estimate 
for  all  points  of  the  curve. 

dince  the  sides  of  the  square  D have  the  length  a - b - 2,  we  obtain 

from  (2.10) 


(2.15)  |*(1)  - 2 

This  eigenvalue  can  be  obtained  by  choosing  JL  m 1,  m«0,  or  Jl  m 0,  m*  1 
and  is,  therefore,  of  degeneracy  2.  We  choose  as  eigenfunctions 

(2.16)  f - cos  x sin  y g - sin  x cos  y 

which  obviously  satisfy  (2.14)  with  • 2 and  are  orthonormalized . We  can 
easily  compute 


(2.17) 


JJVdKU  - 


[ TT  - 81 


■ 


dx  dy 


[TT2*  81 


fxgxdxay  - 0 


o<2  ♦ 62  - i 


We  may  now  use  any  combination 

(2.18)  u - c*  f ♦ |Sg  i 

as  a normalized  eigenfunction  in  obtaining  a supporting  line  at  the  point 
'Z  « 1,  p - 2 by  means  of  (l  24).  We  find  for  the  slope 

(2.19)  a - ^ ("|^)  dx  dy 


and  by  (2.17),  (2.18) 


(2.20) 


TT 


r 

\ 


*2(  IT2-  8)  ♦ ^2 (tt2  + 8)| 


The  extreme  supporting  lines  through  the  point  considered  have  the  slopes 

(2.20  ) 


a - 1 ♦ -8- 
1 T2 


a - 1 - — 
a2  2 
TT 


Thus,  we  have  the  inequalities 


- 58  - 


(2.21) 


>(C)  ♦ (l- A:)  , if  t^i 


TP 


7T 


p(T  ) ^(i- -^)1  Mi* , if  tM 


TP  IT2 

for  the  second  eigenvalue  of  a rhombus  whose  diagonals  have  the  length 
-l/2 

T.  A'*u  and  rr. 

le  shall  now  determine  the  value  of  p(0)  and  the  direction  of  the 
supporting  line  through  this  point.  The  convexity  of  the  curve  p(TT ) 
will  lead  us  then  to  another  estimate  valid  for  the  entire  curve.  Since 
the  first  eigenfunction  of  the  system  (2.14)  is  a constant  and  each 
higher  eigei  function  must  be  orthogonal  to  it,  clearly  if  ^(xjy)  is  the 
second  eigenfunction,  we  must  have 


(2.22) 


2 dx  dy  - 0 


We  may  characterize  the  second  eigenvalue  by  the  extremum  problem 

jy<v«) 


i)2  dx  dy 


(2.23) 


r 


min 


r 


dx  dy 


among  all  not  identically  vanishing,  continuously  differentiable  functions 
JLn  D+C  which  satisfy 

(2.24) 


fr 

1 u dx  dy  * 0 


Consider  now  the  rhombus  whose  diagonals  lie  on  the  x-  and  y-  axis 
and  have  the  length  tP  and  T,  respectively.  Consider  an  arbitrary 
function  u(x)  which  is  odd  and  continuously  differentiable  in  it3  single 
variable.  Because  of  the  symmetry  of  the  rhombus  in  the  y-axi3,  the  conditon 

(2.24)  is  guaranteed.  Using  u(x)  as  a comparison  function  in  (2.23),  we  find 


u'(x)  dx  dy 


dx  dy 


(2.25) 
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Using  the  geometry  of  the  rhombus,  we  can  simplify  and  obtain 


'IT 
— t 


(2.26) 


J u'(x)2(^--i)dn 

ft 

I 


u(x)2 (-?•  - •^■)dx 

«-  t. 


We  introduce  the  new  variable 

(2-27)  S'  1 ' m 


TT 


which  runs  from  0 to  1 as  x runs  fror  ~t  to  zero,  Clearly,  we  can 
choose  any  continuously  differentiable  function  v(C')  defined  in  the 


interval  <^0,l')  and  vanishing  at  ^ ■ 1 and  define 
(2.28) 


u(x)  = v(l-^j^)  for  t 


and  continue  u(x)  as  an  odd  function  into  the  second  half  of  the  rhombus. 

Each  v(9  ) leads  to  an  estimate: 

1 

/Vv'(<p2  df 


(2.29) 


P6"  i 


/fvd^ 


_2. 2 
TT  t 


We  obtain  the  best  estimate  if  we  choose  v(<2  ) in  such  a way  that  the 
ratio  on  the  right  side  of  (2.29)  becomes  a minimum.  This  is  an  elementary 
problem  of  the  calculus  of  variations;  the  Euler- Langrange  equation  for 
v(^)  is 

(2.30)  Y*  v * v ' * ^ v - 0 

with  the  boundary  conditions 

(2.31)  v 2 (0)  - 0 v ( 1 ) - 0 

The  first  condition  is  the  natural  boundary  condition  for  ^ ■ 0 while  the 
second  condition  was  imposed  by  the  requirement  that  u(x)  be  an  odd  function. 
Jv  is  the  minimum  value  of  the  ratio  considered. 
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1 


He  derive  from  (.2.30)  that  the  best  possible  v(^')  has  the  form 

(2.32)  v(^)  - Jq( A^) 

where  J (x)  is  the  Bessel  function  of  order  zero.  In  order  to  fulfill  (2.3l), 

r\ 

we  have  to  choose  -A.  as  the  first  root  of  JQ(x).  He  have 

(2.33)  A -2.4048 

and  (2.29)  may  be  put  into  the  form 

, ,2  _ 

(2.34)  \x  ^ L , T - t 

TT 

This  estimate  determines  a straight  line  through  the  origin  with  slope 
a "2.3426  under  which  the  whole  curve  pi(T  ) has  to  lie.  We  utilized 
till  now  only  the  minimum  definition  (2.23)  of  the  second  eigenvalue.  If 
we  add  now  the  convexity  of  the  curve  the  following  additional  estimate 
arises.  The  curve  u(  'L  ) satisfies 

(2.35)  u(/n^’2T  for  O^T^l 

since  the  line  2T  is  a secant  of  the  curve. 

The  straight  lines  (2.2l),  (2.34),  and  (2.35)  bound  in  the  ( pA  )-plane 
a narrow  triangle  inside  of  which  the  curve  pi( T ) has  to  lie  as  long  as 
0 4T^1.  It  is,  on  the  other  hand,  sufficient  to  know  the  curve  ^(T  ) 
for  values  0 4-T^i.  In  fact,  if  we  perform  a similarity  on  the  domain  D 
with  eigenvalues 

(2.36)  x'-c*x  y ' - c*  y 

-2 

we  obtain  a domain  D'  with  eigenvalues  o<  Xy  as  is  easily  seen  from  the 
system  (2.14)  by  transplanting  the  eigenfunctions.  In  particular,  the 
rhombus  with  diagonals  tTT,  IT  goes  over  into  another  rhombus  with 
diagonals  IT,  ^TT  if  we  choose  c*  - t ^ . Hence,  we  have  the  functional 
equation 

(2.37)  p-(^)  - fx( TT  ) 

which  defines  ^u(T)  for  values  T 1 if  pi(T  ) is  already  known  for 
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Combining  (2.34-  and  (2.37)  we  find  the  estimate 

A2 

(2.38)  U(  1')  ^ - 2.3438 

TT 

I 

valid  for  all  values  of  T" . ^i^iiariy^  we  can  state  that 

(2.39)  p ( C ) ^2  for  1 

These  two  estimates  together  with  (2.2l)  delimit  the  curve  |u(  TT ) for 
T?  1.  [Cf.  Fig.  I.  1 


2*4-  »e  may  use  the  artifice  of  the  preceding  section  in  order  to 
obtain  a general  statement  concerning  the  asymptotic  behavior  of  the 
eigenvalues  for  the  problem  (2.1/*)  as  the  parameter  t of  the  stretching 
tends  to  zero  or  infinity.  In  fact,  let  D-D(l)  be  given  and  let  f(x) 
be  continuously  differentiable  in  this  domain  and  3ati3fy  the  condition 


(2.40) 


f (x ) dx  dy 


0 . 


If  D(t)  is  the  domain  which  is  obtained  from  D by  the  stretching  (2.6) 
we  define  in  it  the  function 


(2.41)  F(x)  - r(-J-) 

This  function  satisfies  the  condition  (2.40)  with  respect  to  its  domain 
of  definition. 

Hence,  by  the  extremum  definition  (2.23),  we  may  assert  the  inequality 

1 2 


(2.4 2) 


u(t)  ^ 


f ' (x)  dx  dy 
^dx  dy  t 


.JL 

2 


for  the  first  non-zero  eigenvalue  of  the  problem  (2.14).  Introducing 

— 2 

again  the  variable  Z • t , we  have 


(2.4 2')  p(t)4  K Z 

where  K is  an  appropriate  constant.  This  result  shows  that  tends  to  zero 
with  T and  is  bounded  from  above  by  a linear  function  of  T if  T"  tends  to 


infinity. 
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It  is  quite  interesting  to  compare  the  lowest  eigenvalue  A of  the 
problem  (2.5)  in  its  asymptotic  character  with  the  eigenvalue  pi.  For 
this  purpose,  we  consider  two  concentric  rectangles  with  sides  a,  b,  and 
A,  B,  respectively,  such  that  the  first  lies  in  the  domain  D while  the 
second  rectangle  contains  this  domain.  From  the  minimum  definition  of  A 
follows  immediately  that  this  eigenvalue  decreases  if  its  domain  of  defini- 
tion increases.  Hence,  we  find  the  inequality  between  the  lowest  eigen- 
values of  the  three  nested  domains,  considered: 

(2.43)  rr2(^  * ^TT2(-^  * 

A 3 a b 

If  we  subject  now  the  entire  plane  to  the  stretching  (2.6),  the 
relation  between  the  corresponding  image  domains  remains  unchanged 
such  that  (2.43)  leads  to  the  new  inequality  (with  T »t  ) 

(2.44)  * Ar)-=  A ^tt2(^  ~ 

A 3 a b 

We  see  that  A(t)  does  not  tend  to  zero  with  T . If  we  plot  A against  T , 
the  curve  will  lie  between  the  two  straight  lines  given  by  (2.44). 

It  is  not  difficult  to  extend  these  considerations  also  to  the  case 
of  higher  eigenvalues;  the  natural  tool  is  provided  by  Poincare's  theorem 
of  Section  II . 1.2 . 

2.5.  Let  us  consider  a domain  D which  appears  from  the  origin  0 


- 63  - 


We  subject  now  the  whole  plane  to  an  angular  stretching 
(2.46)  r ' " r , * t cp  , 0 ^ t < 2 tt  b ^ 

as  defined  in  Section  I,  51.  We  obtain  a sequence  of  domains  D(t)  and 
corresponding  functionals  Rf $ 3 with  respect  to  them.  We  may  express 
these  functionals  also  with  respect  to  the  original  domain  D-D(l)  by 
referring  the  independent  variables  back  by  means  of  (2.46).  We  obtain 
a sequence  of  functioned  in  the  fixed  domain  D: 

l-2 


(2.47)  Rt^] 


J $2/^t'2dr2*  r2dq>2*  jj  [ <$2  ♦ $ lr  dr  d^ 

r>  r ¥ 


JI*2 


r dr  d 


The  eigenvalues  of  are  the  same  as  the  eigenvalues  of  the  functional 

R[c£>1  with  respect  to  the  domain  D(t). 

-2 

If  we  put  again  T.  ->  t , we  observe  that  Hypothesis  I of  Section  1.2 
is  fulfilled.  Hence,  we  proved: 

Theorem  3.  Under  angular  stretching  with  the  parameter  t,  the  sums 


are  convex  from  above  as  functions  of  7T 


for  arbitrary  N and  k. 


We  illustrate  this  result  by  the  following  application.  Let  D be  the 
circular  sector 


(2.48)  O^r-^1  , Off  cf4z.  — , a > 1/2 

The  lowest  eigenvalue  of  the  membrane  problem  (2.5)  for  this  domain  belongs 
to  the  eigenfunction 

(2.49)  u^(r^cp)  - *Jfl(c*r)  sin  a pp 

where  J is  the  Bessel  function  with  index  a and  cX  is  its  first  root.  The 
a 

eigenvalue  corresponding  to  the  function  (2.49)  is 

(2.50)  - <*2 
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The  first  root  p<  of  J is  an  analytic  function  of  the  index 

a 


(2.51) 


OC  = A ' a ' 


imce  the  domain  D belonging  to  t.  is  obtained  from  D by  angular  stretching 
a 1 

with  the  parameter  a \ we  derive  from  Theorem  3: 

2 2 

The  function  A "3 1 is  convex  from  above  as  a function  of  a . 

Thus,  our  convexity  theorem  for  the  eigenvalues  of  the  membrane  equation 
leads  to  a statement  on  the  behavior  of  the  roots  of  Bessel  functions  with 
continuously  varying  index*.  More  results  of  this  form  can  be  easily  obtained 
if  we  consider  various  velues  of  N and  k in  the  above  theorem. 

It  is  clear  that  similar  considerations  can  be  carried  out  for 
eigenvalue  problems  in  more  than  two  independent  variables.  We  leave  the 
formulation  of  the  corresponding  convexity  theorems  to  the  reader. 

3 . Applications  Ii0_conlor3ial_¥flppiog  • 

3.1.  We  consider  again  a smoothly  bounded  finite  domain  D in  the 
(x.y)-plane  and  a function  f(z),  z--x*iy,  which  is  regular  analytic  in 
the  closure  of  3.  Consider  the  conformal  mapping 
(3.1)  z * z ♦ tf (z ) 

which  depends  on  a real  parameter  t.  If  t is  small  enough  this  mapping 
will  be  univalent  and  will  carry  the  domain  D-D(O)  into  a domain  D(t). 

We  shall  consider  the  membrane  eigenvalue  problem  (2.5)  with  respect 
to  the  different  domains  D(t)  and  will  have  to  deal  with  the  functional 

JJ  ( V $ )*  dx  dy 


(3.2) 


Rf  cj>l 


If  & d*  dy 

considered  for  functions  <^(x,y)  of  the  class  ^ in  D(t)  which  vanish  on 
the  boundary  C(t).  We  may  refer  back  the  independent  variables  x,  y in 
D(t)  into  D(0)  by  using  the  transformation  formula  ( 3 . 1 ) . ^ince  the  Dirichlet 
integral  in  the  numerator  of  Rf^>l  is  invariant  under  conformal  transformation, 
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? arrive  at  the  sequence  of  functionals 

dx  dy 


Jf  (^$)2 


(3.3) 


Rt:$i 


r ^ dx  dy 


|1 ♦ tf  '(z)l 


i- 


These  functions  are  to  be  considered  now  for  the  fixed  domain  D»D(0). 

The  eigenvalue  problems  connected  with  this  functional  lead  to  the 
differential  system 

(3.4)  Au*  A j 1 * tf'(z)  I2  u • 0 , u • 0 on  C 

Thus,  instead  of  working  with  the  fixed  differential  system  (2.5)  for 
the  variable  domain  D(t),  we  are  using  the  variable  system  (3-4)  for  the 
fixed  domain  D.  One  advantage  of  this  method  of  investigation  is  that  we 
do  not  need  to  require  that  the  transformation  (3-l)  be  one-to-one  in  the 
complex  plane  but  can  admit  all  analytic  mapping  (3.1 ).  We  are,  of  course, 
essentially  interested  only  in  domains  D(t)  which  lie  simply  over  the  complex 
plane;  but  it  is  often  useful  to  connect  such  a domain  D(t)  with  an  initial 
domain  D’D(O)  by  a continuous  deformation  (3.l)  which  is  not  univalent  for 
all  intermediate  values  of  t. 

Since 

(3.5)  r - 1*  2t  Re  {f'(z)j-  ♦ t2|f  '(z)|2 

we  see  that  r is  convex  from  below  in  t and  that  Hypothesis  II  of  Section 
1.4  holds.  Hence,  we  can  draw  the  conclusion: 

Theorem  4.  If  X^,(t)  are  the  eigenvalues  of  the  system  (3.4)  ordered 

N,  -i 

in  increasing  magnitude,  then  the  sum  Xv(t)  is  convex  from  below  as 

V 1 

a function  of  the  parameter  t. 

N 

The  curve  y„(t)  “ Xy(t)  has  at  every  point  a supporting  line 
-pTl 

with  the  slope 

(3.6)  aN  =*  2 1 f [Reff'(z)|*  t|f'(z)|2l  Xy(t)  1uydxdy 

L J V - 1 
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3.2.  In  the  special  case  which  we  consider  now.  we  can  make  a stronger 
statement  than  Theorem  4.  Let  «f\  ( v » 1,2,  . . . ,N)  be  the  first  N eigen- 
functions of  and  X-^/O)  the  corresponding  eigenvalues.  We  put 


(3.7) 
and  have 

(3.8) 
with 

(3.9) 


- - H -7=% 

v ■ i ^ Ap(o) 


N 

H 

R fu1  - 
1 ZZb 


X X 

ik  i k 


bik  ’ * tf/(z  'I2  ~T= 

D Vi 


§lK 


(o)  Ak(o) 


dx  dy 


Just  as  in  Section  1.4,  we  conclude  from  Poincare's  theorem 


(3.10) 


y-  l 


||l  ♦ tf '(z)l2  2Z  Xy(0)’1iydxdy 

D V-l 


W-?  have  by  Schwarz's  inequality 

N 


(3.11)  (Jj  Re{r(z)j-  yr  Av(0)"I$^2dxdy)2 

D V m 1 

A [f  ^y(0)'1^dx  dyjj|  f^z)]2  YUV0'1'1  &ydx  dy 

D V 9 1 D 2/" 1 

~ (YN(0)l2  fl  lf,(z)!2  t Xy(0)"1^dxdy  , 

4/  v - i 


N 


z^-1,1/2 


where  we  put 

(3.12)  Y„(t)  - i 7^  Av(t) 

We  derive  from  (3.10)  by  means  of  (3.1l)  the  weaker  inequality 

(3.13)  y„(t)>Y* 


,(0)  ♦ tYN(0)_1  |J Re ^f'(z)]'  ZZA^^dxdy 


Thus,  we  have  shown  the  existence  of  a supporting  line  to  the  curve  Y^(t) 
at  t ■ 0. 
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But  from  the  above  result  follows  the  existence  of  a supporting  line 
at  every  point.  In  fact,  let  zq  be  obtained  from  z by  a mapping  (3.l) 
with  the  parameter  t^.  We  can  express  z-F(zq)  by  mapping  D(tQ)  back 
into  D;  on  the  other  hand,  we  can  express  the  transformation  (3.l)  in  the 
equivalent  form 
(3.15) 

Thus,  using  t * t-t  as  a new  parameter,  we  can  repeat  the  preceding 
argument  and  establish  a supporting  line  at  the  point  t»tQ.  Henc^, 


z - z * (t - t )ffF(z  ) 1 
o oo 


we  have  proved : 

Theorem  5.  The  function  [ 
function  of  the  parameter  t. 


1 1 /2 

VJt)  "I  ' is  convex  from  below  as  a 


•V-  1 


3.3.  Let  D be  the  unit  circle  |z|<l;  then  all  eigenfunctions  of 
the  membrane  problem  (2.5)  have  the  form 


(3.16) 


, , . , .cos  no 

uix.y)  - A J (k  r)  , 

n,m  n n,m  sin  n$ 


where  J (x)  is  the  Bessel  function  of  order  n and  where  k is  the  m^  root 
n n ,m 

of  Jn(x).  The  corresponding  eigenvalues  are 

(3.17)  X * k^ 

n,m 

X has  only  one  eigenfunction  if  it  has  the  form  k and  is  otherwise  of 

o , m 

degeneracy  2. 

Consider  now  conformal  transformation  (3-l)  with 

(3. 18)  T{z)  m * .. . 

N .j 

The  curve  y^(t)  - ~)  ’ Xv(t)  ha3  for  t-0  supporting  lines  with  the  slopes 

^ n 

(3.19)  y(0)'1ujdldy 


If  \ i \ it  is  easily  seen  that 


6S  - 


(1.20) 


21 


, J (k 


n turn. 


•)“A 


n ,m 


2 

n*a 


and  does  not  depend  on  the  angular  variable  . From  the  fact  that  the 
harmonic  function 
(3.21) 


Re 


oo 

jf'(z) j - ^ ro<  ^ r'  cos  v <p  ♦ ^ 


r ^ sin  V<^  "* 


has  no  term  independent  of  Cp , follows  that  a^-0.  Since  the  curve 
7j.it)  is  convex  from  below  and  has  a horizontal  supporting  line  for  t-0, 
it  follows  that  y^(t)  has  its  minimum  for  t-0. 

Consider  next  the  case  that  Aj.  - A..+  ; now  we  can  choose  the 
eigenfunction  u L*i  various  ways,  namely  in  the  form 
(3.22)  u,,  - Jn^n  „r)^b<  cos  n<p*  £3 in  ncp)  , pC 2 * $ 2 - 1 

Now,  we  have  obviously 

(3  .23 ) />  * Uy  - A(r ) + B(  r ) ( cx  2 cos*-  n c^>  * (3 2 sin2  ncp>  + 2o(^cos  n<psin  ncp). 

Inserting  this  representation  into  (3.19)  and  denoting 


(3.24) 


2 

cos"  n cp  dx  dy 
2 

sin  n C^dx  dy 


sin  2n  dx  dy 


P 

q 


2s 


we  find  an  infinity  of  possible  0 lopes 

(3.25)  aN  - 2(  0(2p  * 32q  * 2o<(3s) 

Now  clearly  p*q  3 0;  put  cA  - cos  u>  , (3  * sin  so  and  obtain 

(3.26)  a^,  - i(p  cos  2 to  * a ain  2 00 1 

/~~2  21  12  21 

p+s  to  +2/yp+s  and  the  curve 

y^(t)  has,  in  general,  a corner  at  the  point  t-0.  Clearly,  both  branches 

of  the  curve  rise  from  that  point  and  since  the  curve  is  convex  from  below, 
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we  ae"  that  in  this  case  too  y^(0)  is  the  lowest  point  of  the  curve.  Thus, 
we  proved  for  every  integer  N,  under  th«  assumption  that  f(z)  has  the  form 
(3.18),  the  inequality 

(3.27)  yN(°)-yN(t) 

* 

Let  now  D be  an  arbitrary  domain  which  contains  the  origin  z » 0 
and  can  be  obtained  from  the  unit  circle  by  a univalent  mapping  function 

(3.28)  k(z)  - z*  ^ a zv  ; 

>-  2 V 

in  other  words,  we  assume  that  D*  ha3  the  mapping  radius  1 with  respect 

* 

to  the  point  z-0.  We  can  deform  the  unit  circle  continuously  into  D by 
the  transformation 


(3.29)  7.  -z*t(k(z)-z’’  , O^t^l 

By  virtue  of  (3.27),  we  can  then  assert  that 

N * -l  N . 

<3-30)  ZI  <X„>  7 II  X*1 

V - 1 v-l  v’ 

where  the  right-hand  sice  is  the  sum  of  the  reciprocal  eigenvalues  for  the 
unit  circle 

An  arbitrary  simply- connected  domain  D can  be  obtained  from  the  unit 
circle  by  conformal  mappings 
(3.31) 

such  that  the  center  z - 0 of  the  circle  goes  into  an  arbitrarily  prescribed 


k(z)  - a ♦ a,  z ♦ — ♦ a zn  ♦ . 
ol  n 


point  of  D.  If  we  consider  all  possible  choices  of  aQ.  there  exists  a 

maximum  value  for  a^  . This  value  r^  is  called  the  inner  radius  of  D. 

We  can  assume  that  the  mapping  connected  with  this  maximum  coefficient 

transforms  the  origin  into  the  origin;  for  otherwise,  we  achieve  this  by 

3 translation  of  D.  Next,  we  shrink  D similar  to  itself  in  the  ratio  l*.r; 

2 

each  eigenvalue  }\ v in  D goes  over  into  the  corresponding  }j  #-)\vr^  of 


the  shrunk  domain  D . We  map  the  unit  circle  into  D by  a map  function 
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0.2P).  Hence,  ne  conclude  from  (7.30): 
N N 

(3.32)  53  V-n  *•  EL 

v 1 »- 1 

Thus  we  have  proved 


N 


-1  -2 


Theorem  6 . For  N - 1,2, , the  expression  \ X r 
minimum  value  if  the  do~ain  D is  a circle. 


has  the 


It  is  easily  3een  that  the  circle  is  the  only  minimum  domain  for 
Theorem  6 and  the  inequality  ( 3 . 30 ) . For  we  deduce  from  inequality  (3.10) 
that  the  sum  doe3  actually  increase  with  t except  for  the  case  in  which 
f'(z)  = 0. 


This  fact  leads  to  the  following  consequence.  Let  D be  an  arbitrary 
simply-connected  domain  containing  the  origin  such  that 


(3.33)  \\ Re^f'tz)^^  X"1 

L Jv.i 


J v-  1 


u^dx  dy 


for  3ome  N and  all  functions  f(z)  analytic  in  D with  f (0)-0.  Then  D is 
necessarily  a circle.  In  fact,  from  (3.33)  end  the  convexity  theorem  4, 
we  can  shora  that  D leads  to  a minimum  in  Theorem  6 and  by  the  uniqueness 
statement  for  the  circle  it  follows  that  D must  be  a circle. 
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Chapter  TV 
•Symmetry 


l . Natation .and  .Results . 

Lx  l.  *e  consider  a bounded  domain  D in  a plane  referred  to  rectangular 
coordinates  x,y  and  several  functionals  of  D,  namely  A,  I,  P,  K,  k,  and  A. 

A is  the  area  of  D. 

I is  the  polar  moment  of  inertia  of  D with  respect  to  its  center  of 
gravity. 

P is  the  torsional  rigidity  of  D. 

K is  the  electrostatic  capacity  of  D (conceived  as  a plate,  or  a 

cylindrical  conductor  with  the  base  D and  infinitesimal  height). 

k is  only  defined  when  D is  doubly  connected  (r ingshaped ) , bounded  by 

an  interior  curve  C and  ar.  exterior  curve  C,  . Consider  the  two 

o 1 

infinite  cylinders  perpendicular  to  the  (x,y)-plane  that  intersect 
it  in  the  curves  C^  and  C^,  respectively.  The  capacity  of  the 
condenser  formed  by  these  two  cylinders,  per  unit  length,  is 
called  k. 

A denotes  here  the  first  eigenv.-!)ue  of  a membrane  stretched  over  D and 
fixed  along  the  boundary  ot  u.  Therefore,  A • A ^ , if  A^  is  taken 
in  the  meaning  given  to  it  in  several  preceding  sections;  see  III,  2.5). 

The  values  of  these  functionals  corresponding  to  another  domain  D# 
will  be  denoted  by  A#,  I,,  P#,  K# , k* , and  respectively. 

Tie  shall  assume  that  has  a certain  symmetry,  and  we  shall  distin- 
guish two  cases  in  this  respect. 

If  the  bisector  of  the  coordinate  axes,  with  equation  x*y,  is  a line 


of  symmetry  for  we  say  that  D#  has  simple  symmetry. 
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If  there  is  an  angle  2TT/n,  where  n is  an  integer  ^3,  rotated  through 
which  D,  coincides  with  itself,  we  sny  that  D,  has  Jiighgr  symmetry. 

Let  the  points  (x,v)  and(x,,y,,'  describe  the  domains  D and  D,, 
respectively.  If  these  points  are  connected  by  the  relations 

(1.1)  x 3 ax,  , y 3 bx, 

where  a and  b are  positive  constants,  ne  say  that  D,  is  transformed  into  D 
by  pure  dilatation.  If  the  points  are  connected  by  the  relations 

(1.2)  x - ax,  + by,  , y - cx, * dy, 

where  the  determinant  ad  - be  of  the  constants  a.  b,  c,  d is  positive, 
we  say  that  D#  is  transformed  into  D by  affinity. 

Now,  we  have  completed  the  necessary  preparations  and  can  state  our 
theorem. 


(1.3) 


(1.4) 


(1.5) 


(1.6) 


!Lfi  consider  tno  canes: 

(1)  D,  has  simple  symmetry  and  is  traas formed 

(2)  D#  ^as  hjgher  symmetry  and  is  iiansXarmed 

Fax  both  cas§s  wg  bays  tas  inequalities 


I.Ag 

PA'2 

ia’2 

ia"2 

P A-2 

1*^ 

i.*;2 

k2a_1 

^ ^"2 

ia’2 

k2*;1 

- I.‘I2 

-2 

y 

X. 

IA-2 

-2 

* — 

_ -2 

IA 

k 

-hi..  ^ _IA~.2 

v;2 

into  D pure  dilatation: 
into  d by  any  aCiinity- 


The  elementary  fact  that 


i.a;2^.ia-2 


enters  into  all  the  preceding 


inequalities.  We  can  rephrase  the  theorem  by  saying  that  there  are  two 


cases  in  which  all  the  3even  functionals 
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p'Vr1 

, PI'1 

k2ai_1 

, K'Vr1 

kA2!'1 

, k-Vr1 

> A3!'1 

attain  their  maximum  when  D becomes  D#:  first,  when  D*  has  3iraple  symmetry 

and  D varies  through  all  the  transforms  of  D#  by  pure  dilatation,  and 
second,  when  D#  has  higher  symmetry,  and  D varies  through  all  the  transforms 
of  D#  by  affinity.  In  comparison  with  the  first  case,  the  symmetry  is  more 
refined  and  the  class  of  ad-.issible  domains  is  wider  in  the  second  case. 

The  theorem  yields  both  upper  and  lower  estimates  for  the  ratios  P/P* » K/K* , 
k/k*,  and  an  upper  estimate  for  A/ A*-  We  shall  see  that,  in  fact,  these 
estimates  depend  only  on  the  transformation  that  changes  D#  into  D,  but 
are  otherwise  independent  on  the  shape  and  size  of  these  domains. 

1 .2 . We  consider  a given  domain  D in  space  referred  to  the  rectangular 

coordinates  x,  y,  z and  the  functionals  K,  V,  J,  and  H of  D. 

K is  the  electrostatic  capacity  of  D,  or  rather  of  the  body  B whose 

exterior  is  D. 

V is  the  volume  of  D . 

(1.7)  J - jjj(x2  * y2  * z2  )dx  dy  dz 

D 

provided  that  the  center  of  gravity  of  D is  taken  as  the  origin  of  the 
coordinate  systemj  J is  the  radial  moment  of  inertia  of  D. 

(1-8) 

[cf.  15,  F*  68-691.  The  double  integral  is  extended  over  the  whole  surfaces 
of  D an  element  of  which  is  denoted  by  dS , and  h is  the  distance  of  the 
origin  from  the  tangent  plant  at  dS.  We  assume  here  that  the  origin  is  at  the 


center  of  gravity  of  D and  that  D 13  convex.  (To  assume  that  D is  star- 
shaped with  respect  to  the  origin  would  also  be  sufficient.) 

The  quantities  K# , V# , J# , and  H#  are  connected  with  the  domain  D*  in 
the  3ane  way  as  K,  V.  J,  and  H are  with  D. 

We  assume  that  D*  is  symmetrical  and  distinguish  two  cases. 

If  the  three  bisecting  planes  of  the  first  octant  of  the  coordinate 
system  (with  equations  y-z,  z-x,  and  x-y,  respectively)  are  planes  of 
symmetry  for  D#,  we  say  that  L*  has  simple  symmetry. 

If  D has  the  rotational  symmetry  of  any  of  the  five  regular  solids 
(and,  possibly,  some  planes  of  symmetry  in  addition)  we  say  that  D*  has 
higher  ayiMStry-  (Therefore,  D*  has  higher  symmetry  if,  and  only  if,  it  is 
brought  to  coincidence  with  itself  by  all  the  operations  of  one  of  the 
following  seven  finite  groups:  T,  0,  I,  T^,  0^,  1^,  T(0;  for  the  notation, 

see  [8l.) 

Let  the  points  (x,y,z)  and  (x*,y»,z#)  describe  the  domains  D and  D*, 
respectively.  If  these  points  are  connected  by  the  relations 

(1.9)  x = ax*  , y - fcy#  , z - cz* 

where  a,  b,  and  c are  positive  constants,  we  say  that  D*  is  transformed 
into  D by  pure  411aia.ti2h  • If  the  points  are  connected  by  the  relations 
x - anX.  ♦ a12y*  ♦ a^z, 

(1.10/  y * a22^*  + “23^* 

z - a31x#*a32y**a33z. 

and  the  determinant  of  the  nine  constants  a^q>ai2’ * ’ ' ,a33  is  positive,  we 
say  that  D*  i3  transformed  into  D by  affinity. 

We  £Qn3.i4 SC  cases  : 

(1)  D*  baa  simple  symmetry  and  13  transformed  laic  D by  pure  dilatation: 

(2)  D#  has  higher  symmetry  and  is  transformed  into  D by  any  affinity. 
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In  b^t]}  cases 

(1.11) 


j v‘‘ 

^ JL  ^ H. 

jV- 


As  in  the  theorem  of  the  foregoing  Section  1.1,  both  sides  of  this 
inequality  depend,  ir.  fact,  only  on  the  transformation.  In  the  first  case, 
the  planes  of  symmetry  and  the  directions  of  dilatation  are  in  a definite 
geometric  relation  to  each  other,  yet  no  3uch  relation  is  necessary  in  the 
second  case;  this  also  is  analogous  to  the  situation  in  the  plane. 

LJ.  shall  drop  the  somewhat  cumbrous  notation  x# , y# , z#  and  we 
shall  instead  use  n,  ^ instead,  respectively.  3ometi-es,  when  convenient, 

we  shall  write  x^,  x.  , x^,  _ , and  ^ for  x,  v,  z,  tj,  and 

respectively.  For  example,  the  relation  ( 1 . 10 ) will  appear  in  matrix  form  as 


x ~ 

all  ®12  a13 

(1.12) 

y 

m 

a21  a22  a23 

z 

a31  ®32  a33 

L -J 

We  shall  need  sometimes  the  inverse  matrix  with  elements  (X^  and  use, 
instead  of  (1.12),  the  equivalent  relation 

T 1 r - ' - r 

*]  - I <X 


(1.13) 


'll 

*12 

*13 

(21 

*22 

*23 

31 

*32 

*33 

we 

can 

also 

k 

9 

J 1 

x 

y 

z 


( Y 
ik'Tc 


(1.14) 

The  first  partial  derivatives  of  the  function  f(x.y.z)  with  respect 


to  its  first,  second,  and  third  argument,  usually  denoted  by  f , f , f , 

x y z 

will  be  sometimes  denoted  by  f^(x,y,z),  {^{xfy,z),  f^(x,y,z),  respectively, 
ana  similar  notation  will  be  occasionally  used  for  functions  of  two  variables 


o - 


7,,>  Hown  -•  convert*  i or.  for  *.h--  p resent  chapter  T’-o  icn:;.lr.  of 
integration  is  D,  taken  with  2 or  3 dimensions,  if  the  element  of  integration  i3 
dxdy  or  dxdydz,  yet  it  is  D# , if  the  element  of  integration  is 
or  d<T  do  d t,  . T0  surface  integrals,  as  considered  in  Section  2.5,  the 
convention  does  not  apply.  of  course. 


2 . S/oa&trie  Figldii 

2tl.  We  say  that  the  function  f^(  r,  , ~Cf  ) defined  in  the  three 
dimensional  domain  D#  ad~its  the  symmetry  of  D#  if 
(^.i)  Q{  <^r  C..  <^3)  - <f( 

for  an  arbitrary  point  ( C , of  D#  provided  that  the  orthogonal 
transformation  (rotation,  or  rotatory  reflection) 

*$i  * bil  Vl  * fci2  ?2  * bi3  F3  (i-1,-,3) 

transforms  D#  into  itself.  Wp  say  that  the  vector  field  ^3^ 

(i—  1,2,3)  defined  in  D#  admits  the  symmetry  cf  D#  if,  for  an  arbitrary 

point  ( ^ of  . for  the  sane  group  of  transformations,  and  for 

i’  1,2.3, 


(2.2)  Kt(  cr1,  f2,  <>  3 


’ ?3^  * bi2*V?l’  ^2’  T3) 

K3  <?!•  ?2’$V 


For  example,  if  the  function  Cp  admits  the  symmetry  of  D#,  also  the  field 
of  its  gradient  Cp<~  CP,  C Py  admits  that  symmetry.  Similar  definitions 

» I v 

hold  for  scalar  and  vector  fields  in  a plane.  We  wish  to  prove  the  following: 

(a)  The  Inaction  qp(  ^ the  rector  tiib  casapoaenia 

k2^’  ace  defined  in,  and  admit  the  syaoetry  of,  Die  tao  dimensional 
domain  d#  if  d#  has  simple  ayametcy • then 

(2-3)  - SI • 

)lr-  ■ n f ^2C ? ’ 7 5 l2d<S  d*7 


(2.4) 
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9 .un>  nigner  symmetry.  we  have  in  addition 

^•5)  ■ 0 • 

(2.6)  J o J \(  f • ‘"l  ' ^o(  ^ ) d 9 d 7 - 0 

(b)  The  limstifi.0  (|)(  ^ » rj  w ) an3  the  v^sioj:  with  components 
<^(  9 • ri  > C) , K0(  v' . ') . y ) . K^(  C , r;,^  ) axg  defined  in,  and  admit 

the  symmetry  ol . the  three  diasn^ianai  domain  D. . If  Dw  w simple  symmetry. 


then 

(2.7) 

(2.8) 


■ b^v?^ 

H)  K2d\ d 7 d 4 - BKd?d7d^ 


II,  hawfiyer,  D,  has  higher  ayi&siry,  i»e  have  in  additiaa 


(2.9) 

(2.10) 


^V: d d rt 

l<2  r\-j  d *»  d n d ^ 


f f f K,  K 1 d ’T' d r(  dt; 

J J * ■*•  4 1 


S !i  J <?><?fPod?i  dn d 5 

. S' 


U Ki  K2^d7d7 


0 

0 


In  these  statements,  for  all  integrals  considered,  the  domain  of 
integration  is  D# , which  has  to  be  regarded  as  two-  or  three-dimensional, 
as  the  context  requires,  and  the  same  notation  is  adopted  till  the.  end  of 
Section  2.4,  in  accordance  with  what  has  been  said  at  the  end  of  Section  1.3. 

If  the  function  , r<  ,^)  has  the  symmetry  of  D#,  also  the  field  of 

it3  gradient  has  the  symmetry  of  D# . obviously.  Therefore,  of  the  eight 
relations  asserted  by  (a)  and  (b)  it  will  be  enough  to  prove  four,  ?or 
instance,  (2.9)  is  a particular  case  of  (2.10). 


Insofar  as  they  deal  with  simple  symmetry,  the  statements  (.■» ) 
and  (b)  of  the  foregoing  Section  2.1  are  obvious.  Therefore,  we  begin  by 
discussing  the  case  of  higher  symmetry  in  two  dimensions. 

He  assume  that  the  domain  D#  is  transformed  into  itself  by  the  rotation 
(2.1l)  <^>  - ^ cos  9-  1^  sin  9 , rl'  ~ f sin  9 * ^ cos  9 

where 
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(2.12)  0 ^ © - 2Tf/r.  -c  “TT 
and  that 

r|')  * K.,  ( - , r,  ) cos  ©-  K^(  ^ ,ir|  ) sin  6 , 

(2.13) 

K2(^'»  0 ' ’ * 1 ( \ ’"]  ) sin  © ♦ ^(9  ,<q)  cos  © 

ill  all  points  ( ^ , •'j  ) of  D#;  cf  course , ^ ',  rt  •’  in  (2.13)  are  expressed 
by  (2.1l).  We  introduce  the  abbreviations 

(2.H)  . P , jjKjH.d^dd  -9  . JJ^drj  R 

and  the  indeterminate:;  x,  y;  we  assume  that  x' , y7  are  so  expressed  in 
terns  of  x,  y as  vj7  are  m terms  of  ^ , by  (2.1l).  ^hen,  obviously 

(2.15)  j{[x'  K^',0')  * y'  K,(^',r|')]2d^dr]'  - 

Jj!x  Kl(?>'?)*5rK2(‘?-'i,,2dJd'?  • 

Introduce  r,  as  ner  variables  of  integration  of  the  left-hand  side, 
substitute  for  x7  and  y7  their  expressions  in  terms  of  x and  y,  compare 
the  coefficients  of  x^  and  xy  on  both  sides  of  (2.15)  and  use  (2.14). 

We  obtain  so 

(P-  R/sir^  S - 2Q  cos  9 sun  9 = 0 , 

(P  - R)cos  © sir.  © * 2Q  sin^  0*0  , 

that  is,  a system  of  two  homogeneous  linear  equations  for  the  unknowns  P-R 

2 

and  2Q  with  the  determinant  sin  0 which  is  positive  by  virtue  of  (2.12). 
Therefore,  P-R  and  2Q  vanish;  this  is  (sfs(2 .14/ ) precisely  the  desired 
conclusion  (2.4)  (2.6),  which  contains  (2.3).  (2.5)  as  a particular  case. 

The  theorem  (a)  of  -Section  2.1  is  proved. 

2.3.  The  proof  for  (b)  of  Section  2.1  is  similar.  Let  the  vector  field 
^ , fj  , Xf ),  K2(<^',  ^,9-  ),  *3(^,7,^)  kfiVe  syvme'try  °P  D*  and 

introduce  the  inde terminates  x,  y,  z.  The  quadratic  form  in  x,  y,  z 

(2.16)  JfftKp(^*7’^x*K2^,r/’^y+  f<3^’V?,^')Z]2d^dr/d4 


°9 


obviously  admits  all  the  orthogonal  substitution^  u.  the  group  that  brings 

D*  in  coincidence  with  itself.  In  other  words,  (2.16)  is  a quadratic 

2 2 2 

invariant  of  the  grcup.  It  is  obvious  that  x *y  * z is  an  invariant 
of  any  orthogonal  group  It  is  not  so  obvious,  but  it  can  be  shown  [81 

that  the  seven  groups  listed  in  lection  2.1  have  no  quadratic  invariant 

2 2 2/ 

linearly  independent  of  x *y  * r/' . (This  property  characterizes  the  seven 
groups.)  Therefore,  the  quadratic  form  (2.16)  has,  in  fact,  the  shape 

const . ( x * y * z ) 

Comparing  coefficients,  we  obtain  immediately  the  desired  conclusion  (2.8), 
(2.10)  from  which  (2  ' ' ),  (2.9)  follow. 

If  the  center  of  gravity  of  the  two  dimensional  domain  D#  is  at 
the  origin,  the  function  (x  * rp)/2  and  the  field  of  the  vector  9 , ^ 

(its  gradient)  obviously  admit  the  symmetry  of  . We  infer  from  (a)  that, 
if  D#  has  simple  symmetry 


(2.17) 


iJf' 


1 


d9  - ||  Ofcdv  d^ 


If,  however,  D#  has  higher  symmetry,  we  have  in  addition 
V2.18)  | j s f|d^  dr  - 0 

Similarly,  if  the  center  of  gravity  of  the  three  dimensional  domain 


E,  is  at  the  origin,  and  D#  has  simple  symmetry 

(2.19)  fj$2d^dnd£  - j j r-2d<^d  4 - J [J  $2d9  d rj  d $ 

If,  however,  has  higher  symmetry,  we  have  in  addition 

(2.20)  jj  Jrjjd^d^d^  - |f|^d9dr  d^  -JjJ^d^d^d^  - 0 . 

The  relations  (2.19)  and  (2.20)  express  the  fact  that  the  ellipsoid  of 
inertia  of  D#  around  its  center  of  gravity  reduces  to  a sphere,  We  inferred 
this  from  the  circumstance  that  D#  admits  one  of  the  seven  finite  groups 
listed  in  Section  2.1,  If  D,,  admitted  any  finite  group  different  from  these 
seven,  such  an  inference  would  not  be  valid. 
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2 . 5 . can  use  the  method  of  section  2.2  and  2.3  in  dealing  with  H#. 
In  fact,  we  may  assume  that  the  boundary  of  the  three-dimensional  domain  D# 
is  given  by  the  equation 


(2.21) 


( * » •'  * v ) = 1 


where  <£(  r ,»  , fy  ) is  a continuously  differentiable  homogeneous  function 
of  degree  1.  We  have,  th  -refore , 


(2.22) 


■*>  - - 
If  * r 


C > l 


and  so  the  equation  of  a tangent  plane  is 


(2.23)  x ($.  + y$*2.$’l 

» i v 

Let  cos  cs,  cos  p,  cos  $ stand  for  the  direction  cosines  of  the  normal  at 
the  point  f|,  X.  We  take  h to  have  the  meaning  defined  in  Section  2'.1, 
and  we  obtain  from  (2.23)  that 


(2.24) 


h » h * h *«' 


It  follows  from  (1.8)  and  (2.2 4)  that 


(2.25) 


rf(  = 

so  3- 


cos  </•  * $r  cos  jS  + ^ cos  ^ )dS 


This  double  integral,  as  all  those  that  follow  in  the  present  section,  is 
extended  over  the  whole  surface  of  . Consider  now  the  quadratic  form  in 
the  indeterminates  x,  y,  z 

(2.26)  j j*(x  cj)_*  y * 7,  ^>r ) ( x cos  * y cos  + z cos  ^f)dS 

It  obviously  admits  the  symmetry  of  D# . Therefore,  if  D#  has  higher  symmetry, 
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the  quadratic  form  (2.2 A)  must  be  a constant  multiple  of  the  form  x *j  *z 
and  so 

(2.27)  j , <^_cos  v.  dS  * J^i^cospdS  » ^ cos  ^ dS 

(2.28)  j j ($^  cos  y ♦ ^cos  |S)dS  - JJ($c  cos  c*  ♦ C^cos  y)dS 

■ J((£  ^.cos  6 + ^ co3  o^dS  » 0 

Obviously,  (2.27)  (but  not  (2.28))  remains  valid  under  the  different 
assumption  that  has  simple  symmetry. 
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We  shall  need  later  another  form  of  the  integral  (2.25): 


(2.29) 


) 


3.  Affine  Transplantation  of  Symmetric  Fields. 

3.1.  The  point  (x,y,z)  (or  x^,x2,x^)  describes  the  domain  D,  and  the 
point  ( ^ , *) , Xf  ) (or  ( v 1 , ^ ) ) describes  D#.  These  points  are  connected 
by  (1.12)  (or  (1.13),  or  (1.24)),  and  so  the  domains  D and  D#  themselves  are 
connected  by  the  affinity  ( 1 - 12 ) . 

A function  Cp(  > , ^ ) defined  in  D#  determines  a scalar  field  in  D*, 

and  three  functions  <i(^,r/,£),  i’ 1,2,3,  a vector  field.  We  transplant 
these  fields  by  affinity  into  D,  defining  there  f(x,y,z)  and  q^(x,y,z)  by 
the  equations 

(3.1)  f (x,y.z)  - <^>(^,rj,^)  , 

(3.2)  qi(x,y,z)  - ailK1(<^',rj,5)  ♦ ai2K2'(^,y,^)  ♦ ai3K^(^,  J)  , 

i*  1,2,3.  Equations  (3-l)  and  (3.2)  are  valid  at  all  points  of  D#;  of  course, 
x,  y,  z are  expressed  in  terms  of  ( ^»  *}»  $ ) according  to  (1.12)  and  so  f(x,y,z) 
and  q^(x,y,z)  are  defined  in  all  points  (x,y,z)  of  D. 

Taking  the  partial  derivative  of  (3.l)  with  respect  to  x^,  we  obtain 
from  (1.13)  that 


(3.3)  ^(x.y.z)  - {f1(^r),?) 

i- 1,2,3.  Comparing  (3.2)  and  (3.3),  we  see  that  to  take  the  gradient  of 
a transplanted  scalar  field  is  not  necessarily  the  same  as  to  transplant  the 
original  gradient  field.  (The  two  coincide,  however,  when  the  a^  form  an 
orthogonal  matrix,  which  was  the  case  considered  in  Section  2.1.) 

When  D and  D#,  described  respectively  by  (x,y)  and  (f,C|),  are  two- 
dimensional,  the  connecting  affinity  takes  the  shape  (1.2),  (with  9,n  for 


Oj  _ 

x*  ■ y# ' anC*  the  two--diner,siona3  f ields  and  their  transplantation  are 
correspond, ingly  defined. 

We  shall  study  the  trar.sj  Iantation  of  scalar  and  vector  fields  by 
assuming  more  arid  more  symmetry  of  the  original  fields  as  we  progress. 


3.2.  At  first,  we  assune  jic  symmetry  of  D,  • Still  we  can  assert, 
when  D ;u.d  D#  are  t*o-dimens icnal . that 


(3.4) 

*3.5) 

(3.6) 


A = (ad  - cc)A# 


a. 


•ix 


n 


- r<  ♦ r< 


in 


J (q^dy  - q„dx ' - (ad  - bc)J(K^dn-  r^dj) 


the  integrals  in  (3.6)  are  extended  along  corresponding  closed  curves  in 
D and  , respectively. 

When  D and  are  three-dimensional,  we  have 
(3.7)  V “ | a • ! V# 

X K 

°-8)  ql*‘ 

(3.9)  J j (q^dyd,  * q.dzdx  * q^dxdy)  . ;aik  1 , ( K^d^d?  * K^d^d  i^)  ; 

the  Integrals  in  (3  9)  are  extended  over  corresponding  closed  surfaces  in  D 
and  D#,  respectively  - 

We  can  express  (3-8)  by  saying  that  the  divergence  of  a vector  field 
is  an  affine  invariant.  Of  course,  (3.8)  contains  (3.5)  as  a special  case 
and  suggests  strongly  (3.9)  (but.  does  not  prove  it  immediately,  since  the 
interior  of  a closed  surface  that  lies  in  the  field  need  not  wholly  belong 
to  the  field.)  As  the  proofs  for  (3-4),  (3.5),  (3.7),  and  (3.8)  consist 
in  straightforward  calculations,  we  restrict  ourselves  to  a few  remarks  on 
the  proof  of  (3.9). 


We  let  the  parameters  p and  q determine  tin  position  of  a variable  point 


on  the  surface  lying  in  D.  Then,  or,  the  left-hand  side  of  (3.9) 
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(3.1°)  dydz 

and , by  ( 1 . 12 ) , 


dp  dq 


(3.11) 


p q 


a . 


C * 


:i  Ip  “22  r|p*  a23  vp  a21  ^ q * a22  ' q*  ®23 
31  ^P  * a32  r p * a33  ^P  a31  !q  + a32  ^q  * a33  "q 


a22  a23 

Ip  riq 

♦ 

a23  a21 

(?  q 

♦ 

a21  a22 

iP  U 

3 32  a33 

<p  £ q 

a33  a31 

ip  u 

a31  a32 

Ip  Iq 

Since  the  rf  , are  lements  of  the  matrix  inverse  to  that  of  the  a.,  , i»e 
ik  ik’ 

may  express  them  as  cofactors  of  the  latter,  and  nay  therefore  write  (3.1l), 
with  the  notation  (3.10)  as  the  first  of  the  following  three  equations 
dydz  - | aik|  (o(11d^d^*or21d^d^4'5t31d,rdr)  , 

(3.12)  dzdx  - , a^J  (cX12dr;d$  + (X^d^d  fj)  , 

dxdy  - |aik|(°(13dfj  df**23d^  **33dId7) 


Substituting  for  q^,...,  dydz,...  their  values  from  (3.2)  and  (3.12), 
respectively,  on  the  left-hand  side  of  (3.9),  we  obtain  the  right-hand 
side  by  observing  the  relations  between  the  elements  of  inverse  matrices. 


hi-  We  assume  now  that  D#  has  simple  symmetry,  but  we  assume  also  that  the 
affinity  connecting  3 with  D#  is  a pure  dilatation,  given  by  (l.l)  or  by  (1.9) 
(with  9,  rj , ^ for  x# , y# , z# ) according  as  D and  D#  are  two-  or  three- 
dimensional. 


If  D and  D#  are  two-dimensional 
(3.13)  I - 2 ab(a2  + b2)l# 

<3-U>  yf(r^>  ^)dxd,  - f <V  ^ij^2  “^‘V/ 

J'J (q^  + q2  )dx  dy  - 2 ab(a2  ♦ +K^)d^  drj 


(3.15) 
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If  E and  D#  are  three-dimensional 

( 3 .16)  r 1 - 2 


r . 1 

3 


o - - abc(a‘‘  * b*  * c )J, 


(3.17) 


H m Ski  l-U_  «.  -1_  4 -^-)U 

3 ' 2 ,2  2 * 

r.  be 


(3.18)  fJj(f‘-4-h)dxdyd,  . ais  (X  . X . 

• IS)  J * qydxdy  d2  - ^ abc  (a^  * b2  ♦ <3)j  j J (K^  ♦ + d^  d£  . 


f 

\ J 


The  proofs  for  (3-13),  (3.14),  (3.15),  (3 -16) , (3-17),  (3-18),  and  (3.19) 
are  based  or.  (2.17),  (2.3),  (2.4),  (2.19),  (2.27),  (2.7),  and  (2.8),  respec- 
tively. These  proofs  are  so  straightforward  that  we  can  leave  then  to  the 
reader,  who  may  find  it  advantageous  to  use  (2.29)  when  he  deals  with  H. 

3 .4 • Finally,  we  assume  that  D#  has  higher  symmetry  and  we  do  not  restrict 
the  affinity  that  changes  D#  into  D. 

If  D and  D.  are  two-dimensional 


(3.20) 

1*2  “ be  ) (a2  + b2  c2 

(3.21) 

SK-fJ.'dxdy 

(3.22) 

\ \(q?  + q^Jdxdy  - ~ (ad- 

* > 


If  D and  D#  are  three-dimensional 
(3.23) 


ze4j- 


(3.24) 


K ’ 3 * alkl 


(3-25)  |j[(f^rj.d)dxdyd2  • “?f2Kd7d? 

0-26)  XfJtqi  * q2  * q3)d>rdjr  d2-  ” 3 lalk  I2171aik  fXJ(Kl  " K2  * K3)dT dTd  ^ 

For  example,  the  proof  of  (3.25)  results  immediately  from  (3.3), 
(2.7),  and  (2.8),  and  the  proofs  of  the  other  relations  are  likewise 
straightforward,  except  perhaps  the  proof  for  (3.24).  In  order  to  derive 


('’.24)  let  F(x.y,z)  the  homogeneous  function  of  decree  1 that  arises  in 
the  equation 

h.27)  F(x,y.z)  - 1 

of  the  surface  of  D.  which  corresponds  to  ( 2 . 21 ) , the  equation  of  tho 

* 

surface  of  D , Obviously 

(3.28)  F(x,y,7.)  * $(<?  , , X(  ) 

for  nny  point  (^  • ri » 7 ' c*'  provided  that  x,  y,  ana  2 are  expressed  in 

terns  of  , and  y according  to  (1.12).  Then 

(3.29)  F^x.y.z)  - ^ 

for  i = 1,2, j and,  corresponding  to  (2.29). 

(3.30)  H ■ \ f (F  dy  dz  ♦ F dz  dx  * F dx  dy) 

o J x y z 

Fron  (3.30),  (3.29),  (3.12),  (2.27),  (2.23),  and  (2.29),  we  easily  obtain  (3.24). 

The  transfornation  formulas  of  this  section  go  over  into  those  of  the 
foregoing  Section  3.3  if  the  affinity  reduces  to  a pure  dilatation.  The 
difference  between  the  important  formulas  (3.25)  and  (3.26)  co  rresponds  to  the 
difference  between  (3.24;  and  (3.23)  and  i3  due  to  the  difference  between 
(3.3)  and  (3.2). 

The  Dirichlet  integral,  which  is  a conformal  invariant,  is  not  an  affine 
invariant.  The  simplicity  of  the  formulas  (3. IB)  and  (3.25)  is  due  to  the 
underlying  symmetry  assumptions  concerning  D#  and  Cp(^,r],i^). 


4.  Appiisaiioas . 

4.1.  We  apply  now  the  preceding  general  results  to  the  various  functionals 

enumerated  in  Section  1.1  We  start  with  the  torsional  rigidity  P which  was 

treated  in  the  general  case  in  decticn  1.1  of  Chapter  I.  It  was  found 

1 *2 

convenient  to  plot  tP(t)  against  the  variable  Z.  - t if  the  plane  domain  D 


war.  subjected  to  a one-sided  stretching  with  the  parameter  t.  If  f(x,y,z)  is 


the  stress  function  of  D.  ne  expressed  the  slope  of  the  above  curve  at  the 


abscissa  'L  by  (I.  1,  S)  os 


(4.1) 


a - t 3/2p-2  ) 


f*  dx  dv 
x J 


Let  us  suppose  nor  that  D#  is  of  simple  or  higher  symmetry;  since  the 
stress  function  of  a domain  is  uniquely  determined  and  since  D#  goes  into 
itself  under  a symmetry  transformation,  the  stress  function  ,^)  of  D# 

must  be  invariant  under  the  symmetry  transformations.  Thus , (2.3)  must 
hold  in  this  case  and,  by  (I,  14),  ne  have: 


(4.2) 


Thus,  in  this  special  cose,  the  slope  a can  b->  expressed  in  the  form 


(4.3) 


— ~ 3/2  , X 1 .-1/2  p-l 


_ 1 1 2.  - 1 — 

Since  we  plotted  the  ordinate  T ' P against  the  abscissa  TT  , we  may 

characterize  the  slope  of  this  curve  at  the  point  TT  belonging  to  the  symmetric 

* * 

domain  D as  follows.  The  tangent  at  the  point  belonging  to  D intercepts 

the  vertical  axis  at  or.o  half  of  the  ordinate  of  this  curve  point.  This 
fart  leads  to  an  easy  construction  of  this  tangent  if  this  curve  point  is 
known.  ^ince  r,e  ;<r.ow  oir^ndy  that  the  curve  is  convex  everywhere  the 
knowledge  of  such  a point  leads  to  a very  'useful  piece  of  information  on 
the  location  of  the  whole  curve.  This  fact  will  be  utilized  in  the  Appendix 
in  order  to  obtain  numerical  estimates  for  the  torsional  rigidity  of  triangles. 

The  symmetry  of  a domain  can  be  applied  in  the  same  way  in  various  other 
problems  in  order  to  determine  the  slope  of  curves  for  given  functionals  under 
stretching  of  their  domain  of  definition.  Consider,  for  example,  the  electro- 
static capacity  of  a symmetric  domain  D#  in  space.  By  reason  of  symmetry,  we 
have  the  equations  (2.7)  for  the  conductor  potential  cp(^,ry ,^)  of  D#  since  this 
uniquely  determined  domain  function  must  admit  the  symmetry  of  D# . Hence: 
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u-4>  *. 

j * • 

_ i _ 2 

On  the  other  hand,  if  wo  plot  the  expression  t K(t)  against  V ■ t 
in  the  case  of  one-sided  stretching,  we  find  that  the  slope  of  the  curve 
is  given  by  (I.  4.6)  which  leads  now  ta  the  result: 

(4-5)  a = 3 tK#  ’ v"  * 3 ^ ^K*) 

This  can  again  be  interpreted  geometrically:  The  tangent  to  the  curve  at 

the  point  corresponding  to  D,  intercepts  the  vertical  axis  at  2/3  of  the 
ordinate . 

finally,  1 it  D#  be  a plane  symmetric  domain  and  <p(^,Kj)  denote  the 
membrane  eigenfunction  belonging  to  the  lowest  eigenvalue  X#  • 3ince  the 
lowest  eigenvalue  is  non-degenerate  the  eigenfunction  is  uniquely  determined 
and  has  necessarily  the  symmetry  of  the  domain.  Thus,  by  (2.3), 

(4.6)  A*  - •<¥?)d?dr  -2jj<p2d^drj 

We  showed,  on  the  other  hand,  in  Chapter  III  that  if  we  plot  X(t)  against 
_2 

TI  ■ t in  the  case  of  one-sided  stretching,  we  obtain  a curve  with  the 
slope 

(4.7)  a ''■JJ  u^dxdy 

D(t)  X 

* 

Thus , if  TT  is  Just  the  abscissa  of  D#  . we  have  by  (4.6) 

(4.8)  a » p X* 

Z * 

Again,  the  tangent  at  the  point  (T  ,A  ) has  the  intercept  on  the 

vertical  axis. 

It  should  be  observed  that  no  similar  conclusion  is  valid  for  the 
case  of  the  eigenvalue  |0,  of  the  free  boundary  value  problem  (III. 2.14). 

This  eigenvalue  may  be  degenerate  and  admit  more  than  one  eigenfunction. 
Hence,  it  is  not  sure  that  the  eigenfunctions  of  D#  will  have  the  symmetry 
of  the  domain  and  formula  (2.3)  cannot  be  U3ed. 
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4.2.  Wo  consider  next  a dona  in  D#  in  space  with  simple  symmetry.  Let 
^(^>^>4)  denote  its  conductor  potential,  K*  is  capacity  and  K - V^p  the 
corresponding  gradient  vector  field.  Let  D be  obtained  from  D#  by  a pure 
dilatation  (1.9).  We  define  th'*n  a test  function 

(4.9)  r(x,y,z)  - J , *) 

which  is  well-defined  in  D,  has  the  boundary  value  1 and  is  admissible  in 
the  Dirichlet  inequality  (I,  4.4).  Thus,  we  have 

(4.10)  K 


4ir 


Using  now  (3.18)  which  holds  because  of  the  symmetry  of  cp . we  obtain 

(4.11)  K-'  ^ (A,* *^)K# 

a b c 

Another  estimate  for  K may  be  deduced  from  Thomson's  principle  (I.  4.38) 
We  transplant  the  vector  field  by  (3.2)  into  the  domain  D and  use  it 
there  as  a test  field  in  the  Thomson  inequality  which  we  formulate  as 
follows  (see  (I.  4.35)): 


(4.12) 


4TTK 


(Jfq.n  dCT 


dy  dz 

n2 


Hence,  we  find  by  (3.9)  and  (3-19) 

2 + 2 2 

, 4U  | | /vq,  * ♦ c 

(4.13) 


1 ^ 4TriJ/(ql  *q2^1)dxdydz  _ a2  * b2 1 C2  X 

[ jj  (q^  dy  dz  ♦ dx  dz  ♦ q^  dx  dy"]2  3abc 


Thus , we  have ; 

(4.14) 


iafes ^ JL  ^ afes 

2.2  2 — K*  3 2 ,2  2 

a + b * c ^ a b c 


This  may  also  be  formulated  as 


a3  V,3  r3 

JL_+  + 
aK#  oh#  cK# 


y i a 4 „ bK*  c^j. 

<4  JL  ^ 1 (—If f ) 

abc  ""  3 ' 3 .3  3 ' 

a b c 


(4.15) 


- 89  - 


K 

Here  the  lower  and  upper  bounds  for  the  ratio  appear  as  the  harmonic 

C aK,  bK#  cK, 

and  arithmetic  mean,  respectively,  of  the  ratios  ~ , — , ~ ~ . One 

b' 

calls  sometime:;  the  ratio  K.'V  the  specific  capacity  of  a body  if  V denotes 
its  volume.  With  this  notation  (4.15)  may  be  interpreted  as  fellows: 

Under  a pure  dilatation  with  parameters  a.  b.  c the  specific  capacity  lies 
between  the  harmonic  and  ar i thne t : c means  of  the  specific  capacities  which 
are  obtained  under  similarities  with  the  coefficients  a,  b,  and  c, 
respectively. 

If  the  domain  is  the  exterior  of  a sphere  of  radius  1,  the  domain 
D is  the  exterior  of  an  ellipse-id  with  the  axes  a,  b,  c and  (4.15)  leads 
to  the  interesting  result: 

The  specific  capacity  of  an  ellipsoid  with  axes  a,  b,  c lies  between 
the  harmonic  and  arithmetic  means  of  the  specific  capacities  of  three 
spheres  whose  radii  are  a,  b,  c i 14,  Problem  95,  p.  5Tl. 

If  we  compare  (4.14)  with  the  transformation  laws  for  J and  H 
under  pure  dilatation,  we  recognize  that  (4.14)  implies  the  estimate  ( 1 .11) 
which  was  announced  in  Section  4.1. 

In  the  particular  case  that  D#  is  the  exterior  of  the  unit  sphere, 
the  transformed  domains  are  bounded  by  ellipsoids  with  ax»s  a,  b,  and  c. 

In  this  case,  the  capacity  K can  be  expressed  explicitly  in  terms  of 
elliptic  integrals  and  better  estimates  than  (4.15)  csr.  be  obtained,  namely 
[16,  p.  4,  formula  ( 1 . 7 ) t 

(4.15  ' ) (bc)^^  * (ca)^^  * (ab)^^  1^-K^-^[a  + b+  cl 

If  D*  13  the  exterior  of  the  unit  cube,  the  domains  D will  be  the  out- 
sides of  rectangular  boxes  with  3ides  a,  b,  and  c.  Here  symmetri zation 
leads  to  the  result  [15,  p.  1591: 

(4 . 15*' ) K*(abc)l/3^K 
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"he  result  (4.15/  seer:.;  non,  however,  for  the  general  class  of  domains 
with  higher  symmetry.  i?e  mention  as  nn  example  the  tetrahedron  defined  by 
mhe  inequal ties 

(4.15'"'  x * y • z *r-  1 , x ^ " , y "7  0 , z ^ 0 

This  "trirecta.ngular  isosceles  tetrahedron"  has  the  required  symmetry  and 
its  exterior  may  serve  as  0 domain  D# . 


4.3.  In  the  case  that  D#  has  a higher  symmetry  we  can  extend  our  results 

to  the  case  of  ~e.ne:al  affinities.  If  and  ^ are  the  correct 

conductor  potential  and  gradient  field  of  D#  and  if  D has  been  derived  from  D 

by  an  affinity  (l.l')  we  transplant  q>  and  K by  means  of  (3-l)  and  (3.2).  We 

use  the  new  fields  f(x,y,z)  and  q, (x,y,z)  as  test  fields  in  the  ^irichlet  and 

± 

Thomson  principle  for  the  capacity  K.  We  obtain  by  (3.25)  and  the  Dirichlet 


principle 

(4.16) 


4tt 


3 ' aik 


iSE«(kK. 


rhile  Thomson's  principle  (4. If)  and  (3. 3),  (3.26)  yield 


(4-17) 


X51aik  ’ rri 


4TTK  3 lvr’ik'  4-.4l-.-ik  4-n-K* 

Combining  these  inequalities,  we  arrive  at 


(4.18) 


ZI< 


3 1 or 


ik 


ik 


If  we  observe  the  transformations  laws  (3.'7),  (3.23),  and  (3.24)  for  the 
functionals  V,  J,  and  H under  affinities,  we  see  that  (4. 18)  is  equivalent 
to  the  result  ( 1 . 1 1 ) announced  in  dection  4.1. 

Inequality  (4.18)  contains  the  inequality 


(4.19) 


of. 


which  is  just  the  Cauchy  inequality'.  In  view  of  (3.7),  (3-23).  and  (3.24) 


thi3  obvious  inequality  lead.;  to  the  geometric  formulation 
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(...20)  HJV-*" 

Thus,  this  particular  functional  has  a minimum  for  the  domain  with  highest 
symmetry  within  the  class  of  all  affinely  related  domains.  In  the  case 
that  the  dor.ain  D#  i3  a sphere,  we  can  xake  a more  general  statement: 

Among  nil  solids  starshaped  with  respect  to  their  center  of  gravity 
the  sphere  has  the  mininun  value  for  HJV-^. 

In  order  to  prov^  this  statement  we  introduce  the  distance  r from 
the  center  of  gravity  and  the  element  of  solid  angle  dw . fie  have  then 

(4.21)  J - ^ r5dw 

(4.22)  ^ 3 S 1 r^T| 

(4.23)  ^ ^ ^ ] r^dw 

Here  ail  integrations  arc  extended  over  the  boundary  of  the  solid  and  J 

n 

denotes  its  surface  area  . By  Holder ' 3 inequality  with  p*l/3,  q»2/3,  we 
have 

(4.24)  3V  -^(r5)1/3(r2)2/3d.  6(S5r5d.)1/3(Jfr2d.)2/3 
which  means  by  (4.2l)  and  (4.23) 

(4-5)  (3V)3^5JS2 

On  the  other  hand,  it  is  known  [15,  p.  69,  (7)1  that 

(4.26)  32^  3vh 
T hus , (4.25)  leads  to 

(4.27)  HJV"2  £ 2 , {4rrr)(ZJLr!)(^Kr2)  2 

This  proves  the  asserted  minimum  property  of  the  sphere  and  the  above 
derivation  shows  also  that  the  sphere  is  the  only  solid  with  this  minimum 
property . 

4.4.  In  a similar  way  we  can  prove  the  inequalities  ( 1 . 3 ) — (1.6).  fie 
shall  indicate  the  methods  only  very  briefly.  In  the  case  of  the  torsional 
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rigidity  P re  utilize  the  L'irichl-',t  principle  (I,  1.5)  and  the  Thomson 
principle  (I,  1.1S)  and  combine  thee,  with  the  methods  of  transplanting  a 
scalar  and  vector  field.  If  D#  has  simply  symmetry  and  D is  obtained  from 
it  by  a pure  dilatation  (l.l),  we  use  the  identities  (3.5' * (3.14),  and 
0.15)  and  obtain 


(4.2S) 


— 4 — ^ ^ ab  (a2  * b2 ) 

X * _1_  P„  2 aD  D ' 


If  D#  has  a higher  symmetry  and  D arises  through  a general  affinity  (1.2), 
we  apply  (3 .21)  and  (3.22)  and  arrive  at  the  estimates 

(4.29)  ^ (ad  - bc)(a2  * b2  * c2  * d2) 

h * * c"  * d * 


If  we  compare  these  formulas  with  the  transformation  laws  (3.4)  and 
(3.20),  we  recognize  the  validity  of  (1.3).  We  observe  that  all  inequalities 
(l.3)-  (1.5)  contain 

(4.30)  j#a;2^ja"2 


In  the  case  of  a general  affinity  this  is  by  (3-4)  and  (3.20)  equivalent 

with  the  trivial  inequality 

(4.3l)  2(ad  - be ) ^ a2  * b2  * c2  * d2 

_2 

It  is,  moreover,  well  known  that  the  circle  minimizes  the  expression  JA 
among  all  domains . 

He  come  next  to  the  electrostatic  capacity  K for  which  we  use  the 
Dirichiet  principle  (4.10)  and  the  Thomson  principle  (4.12).  We  have  to 
extend  the  linear  transformations  (l.l)  and  (1.2)  into  the  entire  (x,y,z)- 
space;  we  do  this  by  adding  the  equation  z-  with  arbitrary  real  £. 
Suppose  that  D#  has  a higher  symmetry  in  the  (x,y)-plane;  let  f ,7^) 
be  its  conductor  potential  and  K - Vcp  its  gradient  vector  field.  We  use 


the  transformation 
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v4.32;  x = o^*be  , v - c>*  dr  , z -»  i $ 

j j 1 

and  the  transplantation  formulas  (3.l)  and  (3.2). 

From  Dirichlet's  principle  and  (3.2l),  we  derive  easily 

(4.33) 


2 2 0 

K 41  V^- 


2 (ed - bcT 


if  M r?|  fJI?c2dfd,df 


We  put  this  inequality  into  the  abbreviated  form 


( , -1  . N 


? O *>  9 

, „ A ♦ h*' ■ ♦ r 1 . 

^ -c  - / . , \ f ♦ ? vad  - be  1W 

^vaa  - be;  *■ 


and  observe  that 

(4.35)  k#-p*q 

We  utilize  the  freedom  in  the  choice  of  £ in  order  to  minimize  the  right- 
hand  side  of  our  estimate  and  find 

(4.36)  K2  £ ^ (a‘  * b2  * c^  * d^)4Pi 

Since  4PQ  ^ (P  * Q)*1  * , we  obtain  finally 


(4.37) 


^ ^ 1 (a*  * b^ 
v2  - 2 U D 


c2*d2) 


which  leads  by  virtue  of  (3.20)  and  (3.4)  to  the  right-hand  inequality 
of  (1.4).  We  can  derive  the  left-hand  inequality  in  exactly  the  same  way. 
We  prove  also  by  the  same  reasoning  the  validity  of  (1.4)  in  the  case  of 
simply  symmetry. 

In  order  to  prove  (1.5),  we  have  to  characterize  the  quantity  k by  a 
Dirichlet  and  a Thomson  principle.  It  is  easily  seen  that 


(4.38) 


k » min 


_1. 

2T r 


? U <7“>: 


dx  dy 


where  u(x,y)  is  any  continuously  differentiable  function  in  the  ring  domain 
D which  has  the  value  0 on  one  boundary  curve  and  the  value  1 on  the 
other  C^.  The  minimum  is  attained  for  the  conductor  potential  U which  is 
harmonic  in  D.  We  hove  by  Green's  formula 


(4.39) 


k - - 


2TT 


_2J 

an 


ds 


2tt 


jj'  ^ + ^y  ^ 


- 94  - 


Consider  new  any  continuously  differentiable  vector  field  a with 
(4.49)  V - q - 0 

and 


Then  we  have 

1 . 1 rr  2 . . 

k - '“*n  2TT  JJ  q UXOy 

D 

and  the  minimum  is  attained  for  the  particular  vector  field  q - ^ 7U.  The 
proof  for  these  statements  i3  the  same  as  in  the  case  of  the  three-dimen- 
sional capacity  and  is,  therefore,  omitted. 

The  inequalities  (l.S/  are  an  immediate  consequence  of  the  extremum 
principles  ( 4 . 3 S ) and  (4.42)  ard  of  the  transformation  formulas  (3.13),  (3.14), 
(3.15),  (3.4),  (3.20),  (3.21),  and  (3.22). 

Inequality  (1.6)  is  obtained  in  the  3ame  way  by  use  of  the  minimum 
definition  of  Section  III,  2.1,  for  the  lowest  eigenvalue. 


4.5 . The  symmetry  of  a domain  D can  also  be  used  in  order  to  obtain 
estimates  for  functionals  with  respect  to  domains  D(t)  which  arise  from 
D-D(l)  by  the  partial  one-sided  stretching  discussed  in  Section  I,  4.3. 

In  order  to  illustrate  the  method,  we  consider  the  following  example. 

Tie  use  the  notation  of  Section  I,  4.3,  and  assume  that  the  body  B considered 
is  a sphere  of  radiu3  R around  the  origin.  The  body  B(t)  will  consist  of 
a hemisphere  of  radius  R in  the  half-space  x 0 and  of  one  half  of  the 
ellipsoid  with  axes  tR,  R,  and  R in  the  half-space  x ^0.  Let  Kq  be  the 
capacity  of  B and  K(t)  denote  th^  capacity  of  3(t).  If  V(x,y,z)  is  the 
conductor  potential  in  the  exterior  of  B,  we  can  estimate  K(t)  by  means  of 
(I,  4.26)  with  tQ*l.  On  the  other  hand,  we  have  for  reasons  of  symmetry 
(4.43)  V(x,y,z)  ’ V(-x,y ,z) 
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Consider  now  any  continuously  differentiable  vector  field  q with 


(4.40) 
and 

(4.41) 

Then  we  have 

1 , 


y.  q - 0 


27 r 


q-n  da  - -1 


1 . 1 17  2 . 

k “ rain  2~ir  JJ  q aX  Gy 

D 

and  the  minimum  is  attained  for  the  particular  vector  field  VU.  The 

proof  for  these  statements  is  the  same  as  in  the  case  of  the  three-dimen- 
sional capacity  and  is,  therefore,  omitted. 

The  inequalities  (1.5)  are  an  immediate  consequence  of  the  extremum 
principles  (4.3S)  and  (4.42)  and  of  the  transformation  formulas  (3.13),  (3-14), 
(3.15),  (3.4),  (3.20),  (3.21),  and  (3.22). 

Inequality  (l.6)  is  obtained  in  the  same  way  by  use  of  the  ninimura 
definition  of  Section  III,  2.1,  for  the  lowest  eigenvalue. 


4.5.  The  symmetry  of  a domain  D can  also  be  used  in  order  to  obtain 
estimates  for  functionals  with  respect  to  domains  D(t)  which  arise  from 
D->D(l)  by  the  partial  one-sided  stretching  discussed  in  Section  I,  4.3. 

In  order  to  illustrate  the  method,  we  consider  the  following  example. 

17e  use  the  notation  of  Section  I,  4.3,  and  assume  that  the  body  B considered 
is  a sphere  of  radius  R around  the  origin.  The  body  B(t)  will  consist  of 
a hemisphere  of  radius  R in  the  half-space  x 0 and  of  one  half  of  the 
ellipsoid  with  axes  tR,  R,  and  R in  the  half-space  x ^-O.  Let  Kq  be  the 
capacity  of  3 and  K(t)  denote  th^  capacity  of  3(t).  If  V(x,y,z)  is  ths 
conductor  potential  in  the  exterior  of  B,  we  can  estimate  K(t)  by  means  of 
(I,  4.26)  with  t -1.  On  the  other  hand,  we  have  for  reasons  of  symmetry 
(4.43)  V ( x , y , z ) V(-x,y , z) 
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Hence,  if  V is  anv  of  the  partial  derivatives  V V or  V we  have 
i r x y z 


U.44)  V^dx  dy  dz  - JJ  dx  dy 


dz 


From  (I,  4.2")  and  the  equations  (2.?),  we  conclude 


(4 . ) JJj  V^dxdydz  " Jff  ^ dxdydz 


2f  K 

3 o 


D L 

Hence,  we  derive  fron  (I,  4.26)  the  e3ti~ate 

(4.46)  K(t)6  Kj**  £•*  * -.) 

Thus,  we  obtained  a simple  estimate  for  the  capacity  of  the  rather  complicated 
body  B(t)  in  terms  of  the  parameter  of  the  partial  stretching.  Various 
similar  results  can  be  obtained  in  the  same  way  for  analogous  cases  of 


symmetry . 
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Chapter  V 

Trans  • Ir.ntation  of  Harmonic  Functions 

1 • iQESiQAttA  Riciditv  and  the  'reen's  Function. 

IjJ  . In  this  chapter  we  shall  discuss  the  dependence  of  various  domain 
functionals  on  the  parameter  t in  the  transformation 

( 1 . 1 ) Z = Z * t f ( Z ) 

This  transformation  carries  a given  domain  D"»D(0)  into  a sequence  of  domains 
D(t);  it  is  conformal  if  me  assume  that  f(z)  is  regular  analytic  in  D(0)  and 
will  be  one-to-cne  for  sufficiently  small  valuer,  of  t,  if  we  assume  f(z'  to 
be  regular  in  the  closure  of  D. 

Let  H(t)  be  the  class  of  all  functions  h(x.y)  which  are  regular  harmonic 
ir.  D(t);  we  shall  write  for  sak«  of  shortness  h"h(z)  but  observe  that  h does 
not  depend  analytically  or.  the  com}  lex  variable  can  map  the  class  H(t) 

in  a one-to-one  manner  into  the  class  H(C)  of  functions  harmonic  in  D(0)  by 
the  relation 

(1.2)  h (z)  •»  h(z  ♦ tf  (z  ) ) 

o 

In  fact,  the  right  side  function  i3  well-defined  in  D(0)  and  is  harmonic 
there,  since  a conformal  mapping  preserves  the  harmonicity  of  functions. 

We  may  say  that  we  have  transplanted  the  class  of  all  harmonic  functions 
in  L(0)  into  the  class  of  harmonic  functions  in  D(t.)  by  means  of  the 
correspondence  (1.2). 

Various  important  domain  functionals  are  closely  related  to  the  theory 
of  harmonic  functions  in  the  domain  considered  and  their  study  is  facilitated 
by  the  fact  that  harmonic  functions  esn  be  transplanted  in  the  simple  manner 
described . 

1.2 . Let  g(z,£)  be  the  Green's  function  of  a simply-connected  finite 
domain  D;  that  is,  we  suppose  that 
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a)  g(z,  £ ) is  harmonic  in  the  variable  7.,  except  at  the  point  £ 

b)  g(z.  Cr  ' * r^T  log  | rf I is  harmonic  at  rf 

c)  g(z,  ^)  vanishes  if  th*’  variable  z tends  to  the  boundary  C of  D. 

fle  consider  now  the  function 

(-.9)  f(^  ' - 2 11  <:(z.  ^ 'ri^dn  , 

D 

This  function  is  indefinitely  differentiable  in  D gnd  satisfies  by  virtue 
of  b'  the  differential  equation 

(1.4)  Af  - 

By  the  property  c)  of  the  Green’s  function  we  recognize  that  f(z)  vanishes 
on  the  boundary  C and  hence,  we  may  identify  f(z)  with  the  stress  function 
of  D,  defined  in  Section  I,  1.1.  By  means  of  (I,  1.4),  we  compute  the 
torsional  rigidity  of  D: 

(1.5)  P - 4 

This  formula  is  often  a convenient  tool  for  dealing  with  torsional 
rigidity.  It  permits  computation  of  P by  quadratures  if  the  conformal 
mapping  of  D onto  the  unit  circle  is  kncwr. . In  fact,  let 

(1.6)  z » F(w) 

be  the  function  which  is  univalent  and  analytic  in  |w|<l  and  maps  this 
domain  onto  B.  If  £ -F(^->),  we  have  clearly  (transplantation  of  the 
Green’s  function) 

(1.7)  g(z .r/)  ■ icg 

and  we  can  bring  the  formula  0 .5)  into  the  form 

(1.8)  P - f-  ^ ^ log  | |F'(w)|  2|  F'(co)|2du  dvdpdtf" 

\W|  < 1 i sJi  < 1 

if  w»u+iv,  u~>  • p+i(T.  From  ( 1 . 8 ) , we  could  derive  a formula  for  F in 
terms  of  the  coefficients  of  the  Taylor  series  for  F(w)  around  the  origin 

[9,  151. 


gU.  )dxd y d> d 
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1-3  As  another  application,  vie  deal  with  the  following  situation. 

Let  L lie  in  a larger  domain  whoso  Green's  function  and  torsional 
rigidity  vie  denote  by  g,  ( r. . ''  ■ nd  P1  . respectively.  From  the  maximum 

definition  for  the  torsional  rigidity  (I_,  1.5),  we  deduce 

(1.9)  Px  ^ F 

and,  from  the  maximum  principle  for  harmonic  functions,  we  derive  easily 

(1.10)  g1(z,  Xf  ) ~>y  g(",h£  ) 

Obviously,  ( 1 . 9 ) is  a consenuence  of  ( 1 . 10 ) and  the  identity  (1.5). 

Tie  can  obtain  a better  upper  bound  for  P than  ( 1 . 9 ) in  the  form 

(1.11)  P—4  [J  g,  (z,^  )dxdy  drj 

°D  D 

This  formula  is  often  useful  if  the  Green's  function  of  the  domain  D is 

not  available  or  is 'very  complicated  while  the  Green's  function  of  a slightly 

larger  domain  is  easier  to  handle. 

1.4.  Tle  consider  now  the  sequence  of  domains  D(t)'  which  are  obtained 
from  the  given  domain  D»D(C)  by  the  conformal  mappings  (l.l).  Let 
g (z,^)  be  the  Green's  function  of  D(t)  and  P(t)  its  torsional  rigidity. 

Tie  have 

(1.12)  g (s  + tf(z),  tf  * tf(?))  - g(z,^) 

where  g(z,^)  is  the  Green's  function  of  D.  In  fact,  both  si  tj  of  (1.12) 
are  harmonic  in  D and  satisfy  the  requirements  a)-  c)  of  Section  1.2,  which 
determine  the  Green's  function  in  a unique  way. 

We  apply  the  identity  (l.5)  for  the  domain  D(t)  and  its  corresponding 
functionals;  we  then  refer  back  to  D(0)  by  means  c f (l.l)  and  (1.12).  In 
this  way,  we  obtain 

(1.13)  P(t)  - 4jfJ  JJg(z,'4)|l  + tf/(z)|2<|l  + tf/(^)|2dxdyd^dr| 

D D 

This  formula  shows  that  P(t)  is  a polynomial  of  the  fourth  degree  in  t. 


Re  calculate  fro-  ( ':  ]}): 


(l.H)  -f-z  P(t)  - ojj  \\  g(z.^)(lf/(7.)|^|l  * tf'(^)  T 


‘ ' f /(^:(-|l  * t.f ' { r. ) | ^ )d x dy  d r dq  ♦ 32  f f f g(z  ^ ) (Re  £ ' (z  )| 


* 1 1 f ' iz)  | “ ) (Re  ’(£  • t|  f 'iff)  | ^ )dx  dy  d^  d rj 

It  is  a we  1 1- kncwr.  fact  tm*.  thrt  dreer. '.■■  fa  r-naitive -definite ; 


that  is. 

(1.15) 


r. U . y)  £ ( r.  ■ j?  ( dx  dy  d c d r > 3 


- - - > 

for  every  continuous . not  identically  vanishing,  real  function  0(z). 

Hence,  we  have 

(1.16)  P‘(t)>0 

except  for  the  case  that  f '(s'  = '.  Thus,  we  have  proved: 

Theorem  1.  If  the  domains  D(t)  are  obtained  from  a given  domain  D(0) 
by  the  conformal  transfor-.ations  (l.l),  their  torsional  rigidity  P(t)  is 
ccr.vex  from  below  as  a function  of  the  parameter  t. 

This  convexity  statement  could  have  been  obtained  also  from  the 
maximum  definition  (I,  1.5/  by  transplantation  of  the  extremal  function 
f(x,y)  and  by  using  the  invariance  of  the  Dirichlet  integral  under  conformal 
mapping.  The  above  derivation  is  perhaps  more  illuminating  since  it  shows 
the  explicit  dependence  of  P(t)  upon  the  parameter. 

l«i.  ife  maKe  the  following  application  of  Theorem  1.  Let  D be  the 
unit  circle  |zj<l  and  consider  conformal  transformations  (l.l)  with 
f(0)  - f/(0)  - 0.  Re  derive  fiom  ( 1 .13)  for  all  such  transformations 

(1.17)  P'(C)  - 8 


g(z,y)(Re-(f'(z)^+  Re  jf 7 (£)  jldx  dy  d 7 d g 

16  [J  JJ  g(z,^  )Re|f '(^)|dxdy  d7  dd 


D 

(J 

D'  U 


* 


b;. 

virtu 

e of  • 

,he  symr.et 

-1 

o 

(1 

.18) 

Jf“ 

' r ^ 

J 

D 

in 

fact, 

by  t: 

ip  Laplace 

-Foi 

have  the 

a me 

La  pine  inn 

- Re 

)r  ri'f  fir,  - - y Re< zf (z) - ^ f ( z) 
**  . z 


,o 


the  unit  circle,  both  do  vanish  there.  Hut  these  two  facts  clearly 
imply  the  identity  of  both  s id-ms. 

'"e  ob.erve  next  that  ir.  view  of  our  restrictions  regarding  f(z),  we 


have 


(1 .1?) 

i ^ co 

HeJ  r.f(r. ) - ^ f(z  ® *)  ^(a_  (r ) cos  ^ ^ * 

Hence , 

clearly  we  obtain  from  (1.17) 

(1.23) 

p'(o)  - o 

Because  of 

the  convexity  of  F ( t ' from  below,  we  have 

(1.21) 

P'(t)  >3  for  t > 0 

and 

(:  .22) 

P(t)  > F (0 ) 

bv(r)  sin-^C"^ 


Using  the  sane  cons  ideratio.ns  as  in  section  III,  3.3,  we  obtain 
Theorem  2.  Let  r.  be  the  Lnner  radius  of  a simply- connected  finite 
domain  and  Pis  torsional  rigidity.  Then  Pr^  has  the  minimum  value  if 
U is  a circle 


Theorem  2 expresses  the  ineoyuality 

(1.23)  P>/Tri 

valid  for  ail  sinj  iy-connected  domains  E [15"'. 

2.  Furiher  Applications . 

2,1.  Various  other  functionals  are  connected  with  the  Green's  function 
g(z,^)  of  a domain  D.  S7e  mention,  .as  an  example,  the  iterated  Green's  function 

gV^U,^)  ’ JJ  g(z,w)g(w,y)dudv 

D 


(2.1) 


W -*  u + iv 


tir 


t icr. 


the  ireer.'s  function  for  the  biharmonic  differential 


equation 

(...2)  ^ A-.i  ^ " 

with  the  boundary  conditions 

(*. . 3 ) u - A u ’ 7 on  3 . 

It  plays  a role  in  the  listicimy  theory  of  plane  plates  with  the  boundary 

f „ \ 

curve  C.  It  is  easily  seen  that  g'-  \z,Cf)  is  positive  and  continuous  in  D. 

dupiose  now  that  we  transform  D»^(0)  by  conformal  mappings  (l.l)  into 
a sequence  of  domains  D[x.' . Let  g '(z,^)  be  the  corres ponding  iterated 
Green's  functions.  Using  the  definition  C 2 . 1 ) with  respect  to  the  domain 
D(t)  and  referring  bach  to  the  original  domain  D by  means  of  (l.l)  and 


(1.12). 

(..4) 


we 


nr. : 


(z 


g(z  ,w)g(w,!£  )\  1 * tf  '(w  ) | 2 du  dv 
D 


Let  us  choose  an  arbitrary  but  fixed  pair  of  points  z,^  in  D.  The 
functional 

(2.5-  d«(t)  - g^2  ; (z  * tf  (s' , y * tf  U) ) 

(2 ) * * . 
is  the  value  of  g ^ at  the  transplanted  points  z • We  see  from  (2.4) 

3(t)  is  a second  degree  polynomial  in  t.  We  calculate  immediately: 

(2.6)  G"(t)  ~ G"(0)  ■*  2.  jj  g(  z ,w)g(w,"£)(  f '(w  ■ du  dv  >0 

D 


(2.?)  G'(t)  " 2 g^z,w)g(w,£)Re-^f '(w^dudv  + tG"(0) 


We  see  that  G(t)  is  convex  from  below  in  dependence  on  t. 


2.2.  In  order  tc  apply  our  preceding  results,  we  choose  again  D to 
be  the  unit  circle  |z|<l  and  restrict  ourselves  to  functions  f(z)  vhich 

are  analyic  in  the  -unit  circle  and  satisfy  f(0)  -*  f'(0)-*0.  We  specialize 

* * 

further  the  pair  of  points  z,£  by  putting  z*£  = 0.  Observe  that  z - X - 0 


J • '* 


for  nil  m sf.  rings  vl.l)  with  the  oho  .e  functions  f(z  . 
• "i-i  have  in  the  case  of  the  unit  circle 


and  hence 

(2.9) 


, . 


TT 


iog  I n| 


/■  / 


(0) 


A r 

\\  ( lop:  Re(f ' (w)}dudv 

TT*"  JJ  ^ 


In  view  of  the  radial  symmetry  of  the  Green's  function  g(w,3)  and  of  the 
normalization  of  f(w  . we  recognize  that 

(n .10)  :-'(o)  - 3 

Hence,  the  convexity  of  G(t:  implies  that 
(2.11)  :-(t)  5*3(0) 

and  equality  is  cnly  p<jss  ible  if  f'(z)  = 0,  i.e.,  if  D(t)-D(0). 

In  the  usual  way,  we  may  deduce  from  (2.1i)  the  following 
Theorem  3.  If  r.  is  the  inner  radius  of  D and  g^(z,Z.)  is  the 

I 7 

- 2 ( 2 ) \ 

iterated  Green's  function  of  the  ions  in . the  expression  r^  SK  ;(z,z)  has 

a minimum  if  D is  a circle  with  z as  center. 

(2)  , 

It  is  easy  to  calculate  g <0,0 ) for  the  unit  circle;  we  find  for 
its  value  (8T)  * and  may,  therefore,  express  Theorem  3 by  the  inequality 


(2.12) 


J2),  -1-  r2 

g U , - / V-  g-p-  r . 
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Heinz,  Heifer..; tein 


FUdi/Uty  _q£  an  isTdSslQd  Triangle- 

1.  Let  n denote  the  H ire , h the  height,  and  A the  area  of  an 

isosceles  triangle  whose  torsional  r:gi:ity  is  P.  If  ne  keep  a ° 2 fixed, 

the  shape  v f th--  triangle  i.;  deterr.  hie  : by  the  variable  quantity  h.  We 
car.  determine  the  value  of  P .r.  tr.e  following  special  cases : 
infinitely  acute  trian  *!•>:  h --»o>  , P ^ ^ • h ; 

equilateral  triangle:  1.  - V 3 , P ~ ‘/TV  5 ; 


•igr.t  isosceles  triangle: 


n 


, F - 0.104364  ; 


infinitely  flat  triangle:  h— =>C  , P 's-/ h V*6 

The  result.;  for  3 and  h-1  are  well  known  15,  p.  256,  25"’'.  The 

result  for  h — >a;  easily  follows  from  the  value  of  P for  an  infinitely 
narrow  sector  .15,  p.  2**6’.  The  result  for  h — * v can  be  obtained  by  two 
symnetriostions  . In  fact,  for  a right  triangle  with  legs  2 and  h,  P can 
estimated  from  the  P of  a narrow  sector  when  h — >C.  Hymmetr ization  with 
respect  to  a line  parallel  to  th  ’ short  leg  h changes  this  right  triangle 
into  a (flat'  isosceles  triangle,  which  is  changed  b‘  a second  synir.etr ir.ntion, 
with  respect  to  a line  parallel  to  the  base,  into  a rhombus,  the  torsional 
rigidity  of  which  can  bo  estimated  from  above  by  use  of  an  inequality  due 
to  E,  Nicolai  15,  p.  112(l)\  In  view  of  the  fact  that  the  symmetrica tion 
increases  P,  the  torsional  rigidities  of  the  three  figures  considered  turn 
out  to  be  equal  as  h — >). 

vi . We  shall  '.use  four  inequalities  for  the  P of  our  isosceles  triangle 
f 12  \ If  ;ve  introduce  the  notation 

p - 2 ( 3 * hT  , q - {1  • h‘:)(9*  h2)  , r-4(l  + h2)2 


1 -w 

- i - >U  ~ 


tr.er. 


(A.l) 


5-fT 


?(43p3-  ?3pq  ‘ dl6r) 


P > _ ,v2a_3U;p^  - lfrgL 
45 (5p  * l6pq * 64r) 


The  fir^t  (simpler)  upper  bound  for  P ex; reuses  the  fact  tnah  of  all 
triangles  with  a given  area  A th  a -.uilateral  triangle  has  the  largest 
torsional  rigidity  "15,  ; 155'.  *he  second  upper  bound  for  P is  based 

cn  Thor.. on’s  principle;  detail*,  of  the  derivation  ■sill  be  published  else- 
where. The  first  loner  bound  for  P is  a consequence  cf  the  convexity  of 
the  quantity  oJ  ■ h. /P,  considered  as  a function  of  T • l/bX  (see  section  I, 
l.l)  and  of  the  slotw  cf  the  curve  <*>(/C  ) at  51=3  " given  in  (IV,  4.3). 

The  second  iovj->r  bound  is  iiseussod  in  a previous  report  7 1 3 , p.  2i  , formula 
(6.18)  \ 

3.  He  can  derive  closer  bounds  for  P if  we  combine  the  above 
inequalities  with  the  fact  thnt  00  » h/P  is  a convex  and  increasing 
function  of  'C  ■»  l/rT . 

The  four  triangles  of  which  the  torsional  rigidity  can  be  precisely 
computed,  yield  four  points: 

X ” 0 , ^ ■ 3/2  ; 

t - 1/3,  W - 5 ; 

^ = 1 , w = 9.5321  ; 

lin  nr-  •»  6 
1—*  oc 

The  last  relation  yields  the  slope  of  the  asymptote  of  the  curve  to('C) . 

The  inequalities  (A.l)  can  easily  be  translated  into  the  (oo,T)-not3 1 inr. . 


He  write 


\ *Z  * l) 


■>  . 


CaX 


al  iii'c  ♦ iH:  c*  3) 
u 13:  t/  * i%r  * is 


_*  £ :_15  iilC  iJj&lS—l 

- ir..  r 


'n  t*) 


Then,  and  uj  are  ujj'or  bounds,  (*>„  and  u->  lower  bounds  for  uATr)  . In 

*3  .:  4 

a (l\  ui-diagram,  <*->,  represent;  the  conn.cn  tangent  of  ul  , uJ  , and  oj  at  the 

1 ^ 4 

point  ( 1 / 3 , 5 - • The  exj.ression  lox  is  in  other  points  always  less  than  c*j 

(*his  follows  also  from.  the  derivation  of  th°se  quantities)  and  so,  in 

looking  for  an  u.por  bound,  we  always  prefer  to  60^  . By  numerical 

confutation  it  is  found  tint  and  cj,  intersect  (besides  at  (1/3,5)  for 

T * 3.12145,  and  we  have  u>  > c*j  for  0-£  T •<  X , :Ah  ^ sj  for  "T  >/Tl  . 
o ’ 4 2 o’  u 4 o 

The  curve  oj  passes  through  ti'.e  correct  joints  (3,3/2)  and  (l/3,5).  From 

A 

the  convexity  we  conclude  that  the  straight  line  that  joins  these  points 
and  may  be  colled  cur  : another  lower  bound  between  them,  but  ah  and  ux 
are  better  ones.  The  same  is  true  in  comparison  with  the  tangent 
<*>*12.5  Z * 1.5  from  the  point  ( ',3/2)  to  uh. , since  the  tangent  of  has 
a greater  slope  in  this  point  (viz.  12.70).  U)  and  uj  have  a common 
tangent  which  is,  by  convexity,  the  best  lower  bound  for  o>  in  the  neighbor- 
hood of  their  intersection . However,  since  these  curves  intersect  under 
a very  small  angle  the  improvement  is  very  unimportant  and  valid  for  a 


very  small  range  of  values  of  f only. 

The  line  joining  the  precisely  known  points  ( l/3 , 5 ) and  (l,  9.5321) 
is  again  a lower  bound  between  these  points  from  convexity.  But  it  car. 
be  improved  by  drawing  the  tangent 

0>n  - 6.3413 T ♦ 2.7409 
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from  ('  . :.:i ' 

bes  t lower  bound 
dince  do 
little  hi. v.er  (l 
point,  improve  t 
parallel  to  the 
onus tior.  car.  be 


t.o  w which  .ouches  the  latter  curve  at  0./,G06.  The 
in  this  interval  are  therefore  (O  and  <<,. 

•u  not  pe ,u  through  the  correct  joint  (l,  9 . 5B2 1 ) but  a 
, n .64-19 ! we  can,  in  the  left-hand  neighborhood  of  this 
r.e  upper  bound  by  drawing  through  (l,  °.532l)  the  line 
nsyonrtote  of  u(T)  • Since  we  know  its  slope  to  be  6.  its 
found 


^ " 6 T*  * 3 . 1 32 1 

n 

olj  and  ^ intersect  for  T * 1.9133.  By  convexity  is  the  best  upper 

bound  of  .v1  on  the  right-hand  . lie  of  the  point  (l,  9.5821),  as  long  as  it 

is  less  than  IA^  and  u)^  intersect  for  V - 1.3884.  The  best  lower 

bound  for  co  on  the  right-hand  side  of  1 is  ul.  since  ana  cA.  are  both 

<■'  <.  u 

less  than  uj,  . For  large  values  of  f"  the  curve  A can  be  replaced  by  its 
b 3 

asymptote  .<1  - 6 . 54"6  X * 3.3113- 

Sunning  up  we  have  found  the  following  bounds: 


0 

Z ^ 0.12145 

: (A,  ^2 

a) 

c. 12145 

X 4z  0.40064 

0.40064 

0.91327 

40  ^ 

0.91327 

^7' 

1 

£r  Z 1 .33340 

1.38840 

^ X 

Figure  II  represents  the  points  (0,  3/2),  (l/3,  5^  and  (l,  9.5821); 
the  shape  of  the  corresponding  triangle  is  indicated  near  to  the  baseline. 
Essential  portions  of  the  straight  lines  and  ana  of  the  parabola 

ul  are  also  given.  Yet  of  uXy  u)^,  and  40  only  the  points  of  intersection 
with  the  left  or  right  border  cf  the  figure  are  indicated. 


I'*’  - 


iTroro  the  •■luntities  :j  ■*  h/a  - ~ * — and  5 

2 .<fc  ~4 

and  translate  the  above  inequal r * i which  yield 

1 O C - 


•;  UJ-  i • '0 

i.43‘3  > a ^O.voi 
0.791  u>  q ^0.^46 

0.326^  q ^ .••00 


0.500  ^ c ^ .4*4 


0.424 


lx] 


3 ^ 


tj 


12 ,3a  ^ ^ jc  1-  .5a 


to 

2 


lu  * 3 q / ^ s 12  * 3 q 

^ ' a> 

3 7 


6 


/ q / li-xla. 


a; 


2^-0-  zt  Q ^14^-4- 

- cJ 


uJ 


b 
12  r 


^ 3 ^ ^ ■-■  J---4- 

cu 

0 


100  P _ 25 

a-  4 <ff  us 


3y  means  of  these  inequalities  the  Table  I has  been  prepared.  It  is 
remarkable  that  for  ail  values  cf  q the  relative  error  never  surpasses  12.5j6, 
and  that  for  q <1  2 it  is  even  less,  than  6.J.  In  the  neighborhood  of  the 
precisely  knonn  cases  it  drops  iow.n  to  much  smaller  values.  In  about  two 
thirds  of  the  cases  listed  i.n  Table  I tne  relative  error  is  under 

TABLE  I.  Torsional  rigidity  of  an  isosceles  triangle. 

P is  the  torsional  rigidity  cf  ar.  isosceles  triangle  with,  base  a and 
height  h.  The  columns  lists 

(1)  IT  - a^/4h4 

(2)  h/a 

(3)  lower  be. and  H for  1X>  P/a^* 

/ 

(4)  upper  bound  u for  190  P/V 

(5-  approximate  value  2lu/\l  * u)  for  13C  ?/a* 

(6)  extreme  possible  relative  error  10D(u- / )/(u  *1) , computed  in 


percentages 


of  trie  approximation  given  in  column  (5). 


r\.j  _ 


vi : 

U } 

v 3; 

(4) 

(5) 

r 

h/n 

lower 

wrr.er 

100  p/a 

' . '0 

QO 

^ 6.4 '15  h 4o 

'4  3 3 h/a 

~ 7.2917  h 

0.05 

2 ../V>0 

1 0 . Oh 1 5 

13.0053 

11.9171 

o . 10 

1, 5;'il 

6 . ,7. 17 

7.1503 

6.7775 

0.15 

1 .291 0 

/ r ~r;  3 

4.8113 

4 -68«7 

0.20 

1.1180 

3 . 5 32"’ 

3.6081 

3.5702 

0.25 

1 .0030 

2.8671 

2 . 9368 

2 . 3769 

0.30 

0/3 120 

2 .4  ''32 

- .<*  '56 

2.4045 

0.35 

'.-1452 

2 . 061c 

2 . 0620 

2.0619 

0.40 

1 ' .j 

’..'Ol- 

1.9042 

1 .3016 

0 . /,  ^ 

*"7  £ ; 

’ 01/ 

**  O 

t r ~ ~ 

J 

•N  - . > 

^ . >4 

■>  *■7*0  ■» 

1 . 2.2  3 0 

1 .4345 

1.4287 

3.55 

0.6742 

1 ;< 

1 .2958 

1.2897 

v . oo 

0.6/.5-': 

1 .167} 

1 . l'”-'~ 

1.1728 

•J  .65 

0.6202 

1 .C6~7 

1.0785 

1 .0731 

0.70 

0.5-74 

■ . ■>('■>•"  3 

0.9921 

0.9785 

: . 75 

0 . 5 7”4 

0 9080 

0.9160 

0.9124 

0.80 

0.55°: 

0.9430 

*•  or,  ^r' 

- * - ; ~ / 

0.9468 

0.85 

0 • 542  3 

'.'’956 

•' . 3 

7 r»7'2Q 

0.90 

^ r - 

• 7346 

0/2404 

.7375 

0 .9f- 

0.513  ■ 

.6908 

0.6940 

0.69^4 

1.00 

0 . 5000 

0.6523 

0.6523 

0.6523 

1 ,0C 

0.4980 

"*.6146 

0.6172 

0.6159 

1.10 

0.4^67 

0 . 5S0r 

'.6-53 

0.5828 

1 1 c 

J.  • * s 

0 . 4462 

0.5494 

0.5560 

3 . f c 2 7 

' .20 

1 . 4 

'-5210 

0.52-2 

0.5251 

1.25 

O / / '7" 

' .4950 

0.5044 

0.4997 

1.30 

0 .49-15 

' .45*12 

0 .4816 

0.4763 

1.35 

0.4303 

9. 44 91 

'.4605 

0.454 7 

1.40 

' .4226 

0.4289 

0.4409 

0.4348 

1.50 

0.4092 

/nii 

://.056 

0.3993 

1.60 

0.3953 

' . oo-. : 

0. 3’  48 

0 . 3634 

1.70 

"•3  '35 

0.3351 

0 . 34  "3 

0.3413 

oo  o aj  XH  Hah  ^ 122  aiil 

4 275  4 12 

a a 


(6) 

error  i? 


12 . 50 

8.72 

5.21 

2.55 
1.06 

0.34 
0.05 
3 . 31 
0.01 

1.30 

0.40 

0.46 

0.50 

0.51 

1.37 

0.48 

0.46 

0.43 

0.30 

0.23 

0.00 

0.21 

0.41 

0.60 

0.77 

0.94 

1.10 

1.24 

1.38 

1.56 
1 .72 
1.87 


4.36 
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1 1 . Qa  Vy  lar^i-anl Jlisisiity  _£  C_a  Jk  ; tangle  • 

1.  Let.  n,  1)  be  the  jj.de.;  of  a rectangle,  q-a  b their  ratio,  and  P 

its  torsional  rigidity.  '.Ve  keep  b - 1 fixed  and  let  a vary.  P is  known  in 

the  case  of  the  3quare  \n  -*  l),  nomeiv  P - u.1406. 

" o 

Moreover,  the  following  bounds  will  be  used: 


p^‘  a 


1 ■> 

I _a_b_ 

0 J 


which  goes  over  in  an  equation  both  for  a/b  — *0  and  for  a/b — >oo, 


p > p ^0--^ 

r V 1 o ~ > 

0 1 . V*- 

a * b 

T^=  P n2b2 
o 

The  last  inequality  follows  from.  the  fact  that  of  all  quadrilaterals 
with  a given  area  the  square  ha.,  the  highest  torsional  rigidity  C 15 » p.  153 1 . 
The  first  inequality  (an  upper  bound  for  P)  is  a particular  cose  of  an 
ineq>uility  due  to  K.  Nicolai  [15,  p.  112  ( 1 ) 9 . The  middle  inequality  (a 
lower  oound  for  P)  is  a particular  case  of  (IV,  4.23). 


2.  Similarly  as  in  Section  A.  1 . 3 the  cuar.tity  a ^P  * kJ  or  Cl m P 

o 


ap 

T 


r . 


is  a convex  function  of  the  variable  0 

a 

We  define  the  following  f 'unctions: 

Ql  - 3Po(l^)  , n2  * 2 (1  + Z)  ’ °3  = ^ 

Then  the  above  inequalities  take  the  form 

, Q*  n2  , o3 

and  the  precisely  known  cases  are  expressed  by  the  equations 

.Cl  - ?P  for  T •»  0 , and  Cl  = 1 for  T - 1 

c 

In  a T,0  diagram  Cl^  is  the  tangent  of  the  parabola  at  the 
point  (l,l).  From  the  convexity  of  Q( V ) we  conclude  that  the  tangent 
from  the  correct  point,  (0,  3P^)  t.o  Cl^  is  a better  lower  bound  for  .O. 


113 


than  CX^ . The  e;  nation  of  thin  tang^r.t  in 

O - T-V- a'  * lp  - - .5927 X + 3.421? 

!•-  C 

This  tangent  i.;  ever,  a b-'ttor  bound  than  the  lire  joining  the  two 

known  points  (',  '?  un.i  (l.l),  since  the  slope  of  the  latter  is  only 

1-3P  The  point  of  contact  of  ri  and  71  is  X - 0.71 17, 

o • 3 4 

Cl  = 0.?4  36. 

Likewise  the  tangent  to  Cl . parallel  to  w.,  in  a better  lower  bound 

for  O than  . Its  equation  is 

Or  - J?  -Z  * - 1.421? X * 0.5927 

o 

and  its  point  of  contact H • 1 .4052 , Q-*  1.1854. 

‘-‘sing  the  convexity,  the  best  available  bounds  for  Cl  are  therefore 


found  to  be: 

0 ^ X ^ 9.7117 
'■.7117  ^ X i ./.0-2 
1.4152  4:  c 


a4 .4.  o - a2  , 
q1^ci^o2 
o5  ^ o. 


If  in  these  inequalities  we  introduce  the  quantities  q 


a/b,  Q - P/ab\ 


instead  of  T,  Q.,  they  take  the  form.: 


1 . i.  , > /* 


1.1854  ^ q 0.8436 


0.8436 


IP  o' 
c ~ 


2P  q' 

O 


4 Q4 


12pV 

Q 

36P4q2  ♦ 1 
0^ 


4 Q 4 P q 


9 90 

2P  cT  12P  q" 

__C_L  ^ q ^ 2__ 


2 2 
36P  + oT 
o 


Column  (4)  of  Table  II  contains  the  best  approximation  to  the  torsional 
rigidity  P of  a rectangle  that  we  can  compute  or.  the  basis  cf  the  inequalitie 
listed.  Column  (6)  of  the  same  table  shows  the  extreme  possible  error  (in 


orcentages ) cf  the  approximate  value  1 luted  in  column  ( 4 J this  error  is 
estimated  a rr icri . without  reference  fc  t'ne  exact,  value,  only  on  the  basis 


of  the  lc..er  and  upper  estimate.;  listed  in  columns  (•-)  and  (3),  respectively. 
1 he  actual  deviation  of  the  a:;reximato  value  in  colunr.  (4)  from  the  actual 
value  ir.  column  ( f ■ is  -r.uch  less  ir.  most  cases  than  the  extreme  error  given 
ir.  column  (6\ 


more  is  a 


the*,  may  considerably  enhance  the  interest  of 


Table  II. 

Let  P/  denote  the  t 
and  b.  Then  the  foil  on ; 


P 


rsional  rigidity  of  a rhombus  with 
g inenua li tics  hold 


diagonals  a 


which  goes  over  in  an  equation  both  for  a/b  — >0  and  for  a/b  —*■  co, 

P .13 

a*"  * b" 

P aV/7.  ; 

o 

P » 0.1406  is  the  torsional  ririditv  of  the  unit  ssuare.  These  three 
o * 

inequalities  are  derived  from  the  some  sources  as  the  cor  res  pond  Lmg  three 
inequalities  for  P listed  in  Section  1,  and  we  can  summarize  them  by 
saying  that  the  sa”.e  three  inequalities  .mold  for  4 P7  as  for  P.  Therefore, 
the  columns  (l',  (3  . (4).  and  (6)  of  Table  II  yield  1 ower  and  upper  bounds, 
on  apiroximote  value,  and  an  estimate  of  the  extreme  possible  error  of  this 
approximate  value  for  4P/  just  as  well  as  for  P. 

Mow,  4P'  precisely  coincie.es  with  P for  a"b,  a/b— >0,  a/b-^oo.  yet 
\tc  do  not  know  and  can  scarcely  ex;>ecl  that  a precise  coincidence  takes  place 
for  any  other  value  cf  th<->  rati-  a /b . However,  it  seem  reasonable  to  take 
the  (rather  well  known)' value  P/4  as  an  approximation  for  the  (unknown)  P ' . 
Table  II  could  oe  used  tc  estimate  the  error  of  this  approximation  numerical! 
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Ta^ie  : : . Totai.cnai  .Rigidxv** _of  ja JiectaagJie . 

F i.;  the  torsionn'  rigidity  of  a re'tangle  with  aides  a and  b,  a i^.  b. 
The  co limn.;  list: 


(1) 

q ~ 

a/b  ratio  of  Ion 

•or  side 

a to  shorter  side  b. 

(2) 

lowe 

r bound 

JL  for  P/ 

ab5. 

(3/ 

uppe 

r bound 

u for  P, ab  ' . 

(4' 

approximate 

value  2 It 

u '(/♦•-) 

for  P/ab3  , 

(5) 

exac 

t value 

for  F ab 

3 

» 

(6) 

extreme  pcs. 

>ible  relative  er 

ror  100(u-£)/(u‘£) , computed 

(«rc 

entases , 

. of  the 

approximation  given 

in  column  (4). 

(1) 

(2) 

(3) 

(4) 

(5) 

\ 6) 

a/b 

lower 

u;  :er 

Pa_1b"3 

Fa"V3 

error  i 

prox. 

ex. 

1 

0 . 1 406 

0 . 1406 

0.1406 

0.1406 

0.00 

1.1 

0.1540 

0 . 1 54° 

0.1543 

0.1540 

0.23 

1.2 

0.1660 

O.lbS'7 

0.1673 

0.166 

0.32 

1.25 

0.1715 

0.1755 

0.1735 

0.1717 

1.16 

1.5 

0.1947 

0.2152 

0.1998 

0.1953 

2.63 

1.75 

0.2120 

0.2235 

0.2199 

0.2143 

3.75 

2 

0.2250 

0.2467 

0.2353 

0 .228/ 

4.60 

2.5 

0.2424 

0 2722 

0.2564 

0.2494 

5.78 

3 

0.2531 

1.2883 

0 .2696 

0.2633 

6. 51 

4 

0.2647 

0 . 3064 

0.284 0 

0.2808 

7.31 

5 

0.2704 

0.3156 

0.2912 

0.2913 

7.71 

6 

0.2736 

0.3208 

0.2953 

0.2983 

7.94 

7 

0.2756 

0.3240 

0.2078 

0.3033 

8.03 

3 

0.2769 

0.3262 

0.2995 

0.3071 

3.17 

10 

0.27S4 

0.3287 

0.3015 

0.3123 

8.28 

12 

0.2793 

0 3301 

0.3026 

0.3158 

8.34 

100 

0.2812 

0.3333 

0.3050 

0.3312 

8.48 

00 

0.2812 

0.3333 

O.3051 

0.3333 

8.48 

113  - 


Bibliography 

[ 1 ] Bochor,  M.,  Introduction  Higher  Algebra.  New  York,  1912. 

[21  Courant,  R.,  and  iiilbert,  D.,  Methoden  dex  Math.  Phvs ik . Vol.  I, 

Berlin,  1931 . 

H 

C3l  Fekete,  M...  "Uber  die  Verteilung  der  Wurzeln  bei  gewissen  algebraischan 
Gleichungen  mil  gar.zzahligen  Koeffisienten,"  Math.  Zeitschrift.  Vol.  17 
(1923),  pp.  22S-2A9. 

n 

[4"1  Fischer,  E.,  "Uber  quadratische  Fornen  mit  reellen  Koef f izienten,n 
Monatshefte  fur  Math . una  Phvsik.  V0i,  16  (1905),  pp.  234-249. 

[51  Friedrichs,  K.,  "Ein  Verfahren  cer  Variations rechnung,  das  Minimum 

eines  Integrals  als  das  Maximum  eines  anderen  Ausdruckes  darzustellen,n 
Nachr . der  Ges.  der  Missensch.  zu  Gottingen.  1929,  pp.  13-20. 

[6]  Garabedian  P.  R.,  and  ochiffer,  M.,  "Convexity  of  Domain  Functionals," 
Journal  d1 Analyse  Mathematique.  Vol.  2 (1953),  pp.  281-368. 

C7]  Hardy,  G.  H.,  Littlewood,  J.  E«,  and  Polya,  G.,  Inequalities . 

Cambridge,  1934- 

[8]  Meyer,  B.,  and  Polya,  G„,  "Sur  les  symmetries  des  fonctions  spheriques 
de  Laplace,"  Comptes  Rendus . Paris,  Vol.  228  (1949 ),  pp.  23-30. 

[91  Morris,  R.  M.  , "The  internal  problems  of  two-dimensional  potential 
theory,"  Math . Annalen.  Vol.  116  (1939),  pp.  374-400. 

[101  Rooney,  Grove  C.,  "On  the  vibrations  of  triangular  membranes. " 

Technical  Report  No.  35,  Office  of  Naval  Research  Contract  N6ori 
106  Task  Order  5. 

[11"  Poincar4,  H.,  "Sur  les  equations  aux  derivees  partielles  de  la  physique 
mathematique,"  Amer.  Journal  of  Math. . Vol.  12  (lS90),  pp.  211-294- 

[12 1 Polya,  G.,  "Sur  le  role  des  domaines  symm4triques  dans  le  calcul.  de 

certaines  grandeurs  physiques,"  Comptes  Rendus . Paris,  VDp>  235  (1952), 
pp.  1079-1081. 

[131  , "A  simple  method  jf  approximation  applicable  to  some 

problems  depending  on  a minimum  principle,"  Technical  Report  No.  20, 
Office  of  Naval  Research  Contract  N6ori-106  Task  Order  V,  July  7,  1952. 


114 


and  Cl.,  Aufgation  und  Lehrsat/.e  ana  der  Ana_ly3is . 


Vo]  . I..  Near  ;'cr‘:,  1 u.'.  5 • 


. . n r:d 


. , lacpcri-ctric  Inequalities  qT 


Fhvjic-- . Princeton,  1951. 


r i6 " 


_ , and 


, " Inor.ua  2 Lti oj  for  the  capacity  of  a 


ccr.den-or ,"  kacr . J.  af  Matr.c:  m tics  . Vo1..  6"  (’.945),  p p.  1-3?-. 

n 

[I"'  dchiffer,  M.,  and  dzopo,  G..  "Virtual  mans  and  polarization." 
Irar.i.  Ar.er . WfltG-  4qc.,  Voi . 67  (1949),  pp.  133-205. 


STANFORD  UNIVERSITY 


Technical  Reports  Distribution  List 
Contract  Nonr  225(ll) 

(NR- 041- 086) 


Chief  of  Naval  Research 
Office  of  Naval  Research 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  Code  i,} 2 5 

Scientific  Section 
Office  of  Naval  Research 
Department  of  the  Navy 
1000  Geary  Street 

San  Francisco  ?,  California  2 

Attn:  Dr.  J.  D.  Wilke?  1 

Director,  Naval  Research  Lab. 
Washington  25,  D.  C. 

Attn:  Technical  Information 

Officer  9 

Dr.  Edward  Paulson 
Head,  Statistics  Branch 
Office  of  Naval  Research 
Department  of  the  Navy 
Washington  25,  D.  C. 

A l In:  Cone  433  1 

Flsnning  Research  Division 
Deputy  Chief  of  Staff 
Comptroller,  U.  3 A.  F. 

The  Pentagon 

Washington  25,  D.  C.  1 

Headquarters,  U.  3,  A.  F, 

Director  of  Research  and 
Development 

Washington  25,  D.  C.  I 

Chairman 

Research  and  Development  Board 
The  Pentagon 

Washington  25,  D.  C.  - 1 

Asst  Chief  of  Staff,  G-4 
for  Research  and  Development 
U.  S.  Army 

Washington  25,  D.  C.  1 

Chief  of  Naval  Operations 
Operations  Evaluation  Group- 

OP342E 

The  Pentagon 
Washington  25,  D.  C. 


Office  of  Naval  Research 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  Code  438  2 

Chief,  Bureau  of  Ordnance 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn : Re3d  1 

Re6a  1 

Chief,  Bureau  of  Aeronautics 
Department  of  the  Navy 
Washington  25,  D.  C.  1 

Chief,  Bureau  of  Ships 
Asst.  Chief  for  Research  & 
Development 

Department  of  the  Navy 
Washington  25,  D.  C.  1 

Director,  ONR 

Branch  Office 

495  Summer  Street 

Boston  10,  Mass.  1 

Director,  ONR 

Branch  Office 

844  N.  Rush  Street 

Chicago  11,  Illinois  1 

Director,  01® 

Branch  Office 
346  Broadway 

New  York,  13,  N.  Y.  1 

Director,  Of® 

Branch  Office 
1030  E,  Green  St. 

Pasadena  1,  California  1 

Office  of  the  Asst.  Naval 
Attach^  for  Research 
Naval  Attach^ 

American  Embassy 

Navy  No.  100 

Fleet  Fost  Office 

New  York,  N.  Y.  10 


1 


ii 


Commander,  USNOTS 
Pasadena  Annex' 

3202  E.  Foothill  Boulevard 
Pasadena  S,  California 

Attn:  Technical  Library  1 

Office  of  Ordnance  Research 

Duke  University 

2127  Myrtle  Drive 

Durham,  North  Carolina  1 

Commanding  Officer 
Ballistic  Research  Laboratory 


Aberdeen  Proving  Grounds 
Aberdeen,  Maryland 

Attn:  Mr.  R.  H.  Kent  1 

Director,  David  Taylor  Model  Basin 
Washington  25,  D.  C. 

Attn:  L Landweber  1 

P.  Lisenberg  1 

Hydromechanics  Lab . 1 

Technical  Library  . 1 


Director,  National  Bureau  of 
Standards 
Washington,  D.  C. 

Attn:  Natl.  Hydraulics  Lab.  1 

Commanding  General 

U.  S.  Proving  Grounds 

Aberdeen,  Maryland  1 

Commander 

U.  S.  Naval  Ordnance  Test  Station 
Inyokern,  China  Lake,  Calif.  1 

N.  A.  C.  A.. 

1724  F St.,  N.W. 

Washington  25,  D.  C. 

Attn:  Chief,  Office  of  Aero- 

nautical Intelligence  1 

Hydrodynamics  Laboratory 
National  Research  Laboratory 
Ottawa,  Canada  1 

Director 

Penn.  State  School  of  Engineering 
Ordnance  Research  Laboratory 
State  College,  Pa.  1 


Statistical  Laboratory 
Dept,  of  Mathematics 
University  of  California 
Berkeley  4>  Galifornia  1 

Hydrodynamics  Laboratory 
California  Inst,  of  Technology 
1201  East  California  St. 

Pasadena  4,  California 

Attn:  Executive  Committee  1 

Commanding  Officer 
Naval  Ordnance  Laboratory 
White  Oak,  Silver  Spring  19,  Md. 

Att-p:  Technical  Library  1 

Ames  Aeronautical  Laboratory 

Moffett  Field 

Mountain  View,  California 

Attn:  Technical  Librarian  1 

Mr,  Samuel  I.  Plotnick 
Asst,  to  the  Director  of  Research 
The  George  Washington  University 
Research  Laboratory  - 

Area  B,  Camp  Detrick 
Frederick,  Maryland  1 

Mathematics  Library 
Syracuse  University 
Syracuse  10,  N.  Y.  1 

University  of  Southern  California 
University  Library 
3518  University  Avenue 
Los  Angeles  7,  California  1 

Library 

California  Institute  of  Technology 
Pasadena  4,  California  1 

Engineering  Societies  Library 

29  W.  39th  Street 

New  York,  N.  Y.  1 

John  Crerar  Library 

Chicago  1,  Illinois  1 

National  Bureau  of  Standards 
Library 

3rd  Floor  North  West  Building 
Washington  25,  D.  C. 


1 


iii  - 


Library 

Hass.  Inst,  of  Technology 
Cambridge  39,  Hass.  1 

Louisiana  dtate  University  Library 

University  Station 

Baton  Rouge  3 , La . 1 

Library 

Fisk  University 

Nashville,  Tennessee  1 

Mrs.  J.  Henley  Crosland 
Director  of  Libraries 
Georgia  Institute  of  Technology 
Atlanta,  Georgia  1 

Librarian 

G-16,  Littauer  Center 

Harvard  University 

Cambridge  38,  Mass.  1 


Dr.  E.  P.  Cooper 
U.  S.  Naval  Shipyard 
U.  S.  Navy  Radiological  Defense 
Laboratory 

San  Francisco,  California  1 

Prof.  R.  Courant 

Inst,  for  Mathematics  & Mechanics 

New  York  University 

New  York  3,  N.  Y.  1 

Dr.  A.  Craya 

Dept,  of  Aeronautical  Engineering 

Columbia  University 

New  York.  27,  N.  Y.  1 

Dr.  K.  S.  M.  Davidson 
Experimental  Towing  Tank 
Stevens  Inst,  of  Technology 
711  Hudson  Street' 

Hoboken,  N.  J.  1 


Prof.  L.  V.  Anlfors 

Dept . of-  Mathematics 

Harvard  University 

Cambridge  38,  Mass.  1' 

Prof.  Peter  G.  Bergmann 

Dept,  of  Physics 

Syracuse  University 

Syracuse  10,  N.  Y.  1 

Prof.  G.  Birkhoff 

Dept,  of  Mathematics 

Harvard  University 

Cambridge  38,  Mass.  1 

Prof.  H.  Busemann 
Dept,  of  Mathematics 
University  of  So.  California 
Los  Angeles  7,  California  1 

Dr.  F.  H.  Clauser 
Dept,  of  Aeronautical  Engineering 
Johns  Hopkins  University 
Baltimore  18,  Maryland  1 

Dr.  Milton  Clauser 

Dept,  of  Aeronautical  Engineering 

Purdue  University 

Lafayette,  Indiana  1 


I — 


Prof.  R.  J.  Duff in 
Dept,  of  Mathematics 
Carnegie  Inst,  of  Technology 
Pittsburgh  13,  Pa.  1 

Dr.  Carl  Eckart 

Scrippsc Inst,  of  Oceanography 

La  Jolla,  California  1 

Prof.  A.  Erd^lyi 
Dept,  of  Mathematics 
California  Institute  of  Technology 
Pasadena  4,  California  1 

Prof.  K.  0.  Friedrichs 

Inst,  for  Mathematics  & Mechanics 

New  York  University 

New  York  3,  N.  Y.  1 

Prof.  D,  Gi]barg 

Grad.  Inst,  for  Applied  Mathematics 

Indiana  University 

Bloomington,  Indiana  1 

Prof.  A.  T.  Ippen 
Dept,  of  Civil  and  Sanitary 
Engineering 

Mass.  Inst,  of  Technology 
Cambridge  39,  Mass.  1 


F: 

n . • . • • i r . * i . 

Prof.  C.  Debar i 

L 

I s!  *„  i v*.o 

Dept.  of  Mathematic:; 

• f Tec  no 

Washington  university 

■ ■ Cf«  ; f ' ' I'  ".  i 'l 

1 

. Lg- is,  M . 

■ 

*'  • K;  i :k- 

Dr.  ’uL.  o.  Plesr.et 

D y‘  . 

*.  wia  • : •*  ■ v i wi 

California  Inst,  of  Technology 

V:'. . 

. :*.■  f Pennsylvania 

Pasadena  A,  CaLLfor-.ia 

i 

r*.: : :: 

' y:. . n u . Po  . 

* 

Prof.  W . Proger 

Lr  . K. 

T.  r.nap; 

->e  it.  of  Mi  them  tic  a 

. • ,‘r 

ics  laboratory 

Br  -nn  Drivers i ty 

z\  i : 

r::  i>*  lost.  of  Technology 

Providence.  12 , H I. 

1 

Fnsa.le 

no  .. . Cr.  1 for.-,  la 

1 

Prof.  T.  C.  Hcser.bloor 

r , - •> 

P A . Lo  r * r.vr. 

De {.■  t . of  V.a the na t i c s 

Ac rcrn 

D • : • . 

University  of  Minnesota 

C;  i f 

ruin  "rut.  ' f Technology 

.V.inr  e s } s i i .»  1 A . I’inr. . 

a 

Parade 

r.n  ruin 

") 

Dr.  :i.  Pc  use 

Dr  . M, 

rt.  m Lesser. 

Btate  Inst,  of  Hydraulic  Rosea 

rch 

Aer  or.fl 

utiral  Engir.-.ers.:;  y Do  t. 

University  or  Io*a 

F<nr 

state  C- 

Iowa  City,  Icnn 

i. 

L ♦ * «p 

n „ > • . r 

^ . .C,/.  *•:  . 

n 

Prof.  Charles  B^iltzer 

Prof. 

:i . G . Lew 

Case  Institute  of  Technology 

1>  -• ; . 

of  Aeronautical  Engineer 

:ng 

Cleveland,  Ohio 

1 

r^r.r. . 

c - : : o <’o 

* r f.  A 2.  Bchaeffer 

^ ‘a  vr> 

Coil eg'-.  Pa. 

1 

De:t.  of  ‘.'a the na tics 

T V, 

: t - i r. . : L. v;y 

University  of  Wisconsin 

D-:  pt . 

of  ttathe-vtt  :c:> 

Ma'ison  fc.  Wisconsin 

i 

Vr.ive: 

j;iy  of  Cali  for::  La 

LVr  ke] 

try,  Co  1 i for. -i  in 

* 

Prof.  L.  I . ochiff 
Be. t . of  Physios 

I r -'f . 

v w Err: 

Bx-anford  University 

Lc 

of  V..-:  •, he  -at:  or 

Ctanforc,  California 

1 

* ♦ 

Ins  3.  of  Te-imology 

Cn-i  : i 

ore  3?,  Mis... 

1 

Pro'*.  i7 . Bears 

Crad . School  of  Aeronautical 

P:  jf 

W . T . M,rt:r. 

Engineering 

L <■ ; *■  . 

of  M'»the~.:it. 

Cornell  University 

K:,l  . 

Insr. . of  Technology 

Ithaca,  e.  :• 

1 

Ca~Vr ; 

ago  3v.  ^3.-:.;. 

1 

l 

Trof.  Herbert  oclonon 

D an  Paul  E.  Mohr. 

Teachers  College 

3r.  hnol 

of  Er.  g ir.ee  rir.o 

Coluitibia  University 

fr.e  ’J- 

-ver  ' i?y  of  Buffalo 

Mer  Ycr’-  17,  N . Y.' 

1 

Buffo: 

*•  \r 

c,  .. 

1 

Prof.  D.  sponcer 

Tr.  f . 

j Ls , Lc  r r e v 

Pine  Hall , Box  °DS 

D e y .. 

o!’  tnerratici 

Princeton.  !i  J. 

■> 

i. 

\’n  i ver 

. r.v  of  Cal : for n in 

Bert  j 

e /..  California 

V 


I r<  f.  J J . . >toker 

t.  for  Mathematics  '•  Mechanic.* 
Ne  ' V :<  rk  ’Jn  : v*t  i tv 


Dr  . V . L . Streeter 
Fur.dame.ntii  1 Mechanics  Research 
lilin  -is  Inst.  of  T och.no  logy 
Chicago  16,  Illinois 

Frof.  3.  Trues  del  l 
Sr  a 1 Inst  . f.  r A;-;  1 ; ou 

Indian?  l-niver^ity 

3 i o or;  1 i n gtc n . I n<;  \ nr. 1 

Prof.  Joserh  - 1 1 nun 
D p t . of  Ma  t he  Katies 
University  f Michi ■•an 
Ann  Arbcr . Michigan  1 

Prc  f . J.  K.  Venn a rd 

Dey:  of  3 : v i 1 Engineering 

Sta  : 3c  r d Vr.  i ve  rs  i ty 

Stanford.  California  1 

Prof.  Martin  Vitnusek 

University  of  .'Dim i i 

Honolulu  11,  Hawaii  1 

Prof.  3.  E Wars  char  ski 
D ; i of  Mathura  tics 
"niver.  ity  of  Minnesota 
Minneapolis  14.  Minr..  1 

Frof.  A.  Weinstein 

Inst,  for  Fluid  Dynamics  & 

Applied  Mathematics 
University  of  Maryland 
College  Park,  Maryland  1 

Prof.  A.  Zygmu.nd 
Dept,  of  Mathematics 

Jnp  ’’n  1 yr  r j a *y  O f*  lC3  JO 

Chicago  17,  Illinois  1 

Mathematics  Dept. 

University  of  Colorado 
Boulder,  Colorado  1 

Dos  Angeles  Engineering  Field  Office 
Air  Research  and  Development  Command 
5504  Hollywood  Bc'ulevard 
Los  Angeles  28,  California 
Attn:  Cnpt.  f. . E.  Nelson 


Navy  Department 
Naval  Ordnance  Test  Station 
'Underwater  Orlnance  Dept. 
Pasadena,  California 

Attr.:  Dr.  0.  V.  Schliestett 

Co-e  P8101  1 

A3TIA,  Western  Regional  Office 
f504  Hollywood  Blvd . 

Hollywood  28.  California  1 

Armed  Services  Technical 
Infcrmutio.n  Agency 
Documents  Service  Center 
Knott  Building 
Dayton  2,  Ohio 

Mr . R . T . Jones 
Ames  Aeronautical  Laboratory 
Moffett  Field 
Mountain  Vice,  California  1 

Mr.  J.  D.  lierson 
•'l°nr.  D.  Martin  Company 
Middle  River 

Baltimore.  Maryland  1 

Mr.  E.  G.  Straut 
Consoli  iated-Vultee  Aircraft 
Corporation 

Hydrodynamics  Research  Lab. 

San  Diego.  California  1 

Dr.  George  E.  F0rsyt.he 

National  Bureau  of  Standards 

Institute  for  Numerical  Analysis 

University  cf  California 

435  Hilgard  Avenue 

Los  Angeles  24,  California  I 

Additional  copies  for  project 
leader  and  assistants  and  re- 
serve for  future  requirements  50 


M a * t.  i c.; 

i 


1 


\.n 


